\jd^?    &ex>G\2^j&      &o£>t^\ 


f^O;  VBCU 


FWS/OBS-80/40.17 
October  1983 


Air  Pollution  and  Acid  Rain 
Report  No.  17 


ROCKY  MOUNTAIN 
ACIDIFICATION  STUDY 


iA/A^!9!iAL  PARK  SERVICE 
W^TS?  RESOURCES  DIVISION 
RFORTCOLUNS  COLORADO 
RESOURCE  ROOM  PROPERTY 


Office  of  Research  and  Development 
U.S.  Environmental  Protection  Agency    M 

Fish  and  Wildlife  Service 

U.S.  Department  of  the  Interior 


REPORTS   ISSUED 


FWS/OBS-80/40.1 
FWS/OBS-80/40.2 
FWS/OBS-80/40.3 


FWS/OBS- 
FWS/OBS- 
FWS/OBS- 
FWS/OBS- 
FWS/OBS- 
FWS/OBS- 
FWS/OBS- 
FWS/OBS- 


80/40.4 

80/40.5 

80/40.6 

80/40.7 

80/40.8 

80/40.9 

80/40.10 

80/40.11 


FWS/OBS-80/40.1 2 
FWS/OBS-80/40.1 3 
FWS/OBS-80/40.1 4 
FWS/OBS-80/40.1 5 
FWS/0BS-80/40.16 
FWS/0BS-80/40.17 


Effects  of  Air  Emissions  on  Wildlife  Resources 

Potential    Impacts  of  Low  pH  on  Fish  and  Fish  Populations 

The  Effects  of  Air  Pollution  and  Acid  Rain  on  Fish, 
Wildlife,  and  Their  Habitats:     Introduction 

:  Lakes 

:  Rivers  and  Streams 

:  Forests 

:  Grasslands 

:  Tundra  and  Alpine  Meadows 

:  Deserts  and  Steppes 

_______^__ :  Urban  Ecosystems 


Critical  Habitats  of 


Threatened  and  Endangered  Species 

Effects  of  Acid  Precipitation  on  Aquatic  Resources: 
Results  of  Modeling  Workshops 

Liming  of  Acidified  Waters:  A  Review  of  Methods  and 
Effects  on  Aquatic  Ecosystems 

The  Liming  of  Acidified  Waters:  Issues  and  Research  - 
A  Report  of  the  International  Liming  Workshop 

A  Regional  Survey  of  Chemistry  of  Headwater  Lakes  and 
Streams  in  New  England:  Vulnerability  to  Acidification 

Comparative  Analyses  of  Fish  Populations  in  Naturally 
Acidic  and  Circumneutral  Lakes  in  Northern  Wisconsin 

Rocky  Mountain  Acidification  Study 


For  sale  by  tho  Superintendent  of  Documents,  U.S.  Government  Printing  Office.  Washington,  D.C.  20402 


UNITED  STATES 

DEPARTMENT  OF  THE   INTERIOR 

FISH  AND  WILDLIFE  SERVICE 


Dear  Colleague: 

The  Eastern  Enerqy  and  Land  Use  Team  (EELUT)  is  pleased  to  provide  you 
this  report  on  the  evaluation  of  the  sensitivity  and  potential  effects 
of  acidic  deposition  in  watersheds  characteristic  of  the  Rocky  Mountain 
Region.  This  report  is  part  of  the  series  of  technical  reports  on  air 
pollution  and  acid  rain  developed  at  EELUT.  Previous  reports  are   listed 
on  the  inside  front  cover. 

Areas  within  Rocky  Mountain  National  Park  and  Yellowstone  National  Park 
were  selected  as  representative  of  geologic  types  in  a  large  portion  of 
that  in  the  total  Rocky  Mountain  Region.  In  addition  to  determining  the 
sensitivity  of  characteristic  watersheds,  the  study  also  evaluates  the 
extent  of  current  acidification,  the  impacts  on  fish  populations,  and 
recommendations  for  assessment  of  future  trends  in  both  changing  water 
chemistry  and  fishery  impacts. 

Please  feel  free  to  send  suggestions  or  comments  on  this  report  to  EELUT, 


Sincerely, 

R.  Kent  Schreiber 

Acting  Team  Leader,  EELUT 


Digitized  by  the  Internet  Archive 

in  2012  with  funding  from 

LYRASIS  Members  and  Sloan  Foundation 


http://archive.org/details/rockymountainaciOOgibs 


FWS/OBS-80/40.17  Air  Pollution  and  Acid  Rain 

October  1983  Report  No.   17 


ROCKY  MOUNTAIN  ACIDIFICATION  STUDY 
by 


1. 

J.  H.  Gibson 

2. 

James  N.  Galloway 

3. 

Carl  Schofield 

4. 

William  McFee 

5. 

Robert  Johnson 

6. 

Sandy  McCarley 

7. 

Nancy  Dise 

8. 

David  Herzog 

(See  addresses  on  page  ii) 


Project  Officer 

R.   Kent  Schreiber 

Eastern  Energy  and  Land  Use  Team 

U.S.    Fish  and  Wildlife  Service 

Route  3,   Box  44 

Kearneysville,  WV  25430 

Cooperatively  Produced  For: 

Research  and  Development 

Eastern  Energy  and  Land  Use  Team 

and 

Region  6,  Habitat  Resources 

U.S.  Fish  and  Wildlife  Service 

U.S.  Department  of  the  Interior 

Washington,  D.C.  20240 


NATIONAL  PARK  SERVICE 

WATER  RESOURCES  DIVISION 

FORT  COLLINS,  COLORADO 

RESOURCE  ROOM  PROPERTY 


Fish  and  Wildlife  Service 


U.S.   Department  of  the  Interior 


DISCLAIMER 


Although  the  research  described  in  this  report  has  been  funded  wholly  or  in 


part  by  the  U.S.  Environmental  Protection  Agency 

No.  EPA-81-D-X0581  to  the  U.S.  Fish  and  Wildlife 

subjected  to  the  Agency's  peer  and  policy  review 

sarily  reflect  the  views  of  the  Agency.  Mention 

products  does  not  constitute  endorsement  or  recommendation  for  use  by  the 

Federal  Government. 


through  Interagency  Agreement 
Service,  it  has  not  been 
and  therefore  does  not  neces- 
of  trade  names  or  commercial 


1.  J.  H.  Gibson 

Natural  Resource  Ecology  Laboratory 
Colorado  State  University 
Fort  Collins,  Colorado  80523 

2.  James  N.  Galloway 

Dept.  of  Environmental  Science 
University  of  Virginia 
Charlottesville,  Virginia  22903 

3.  Carl  Schofield 

Department  of  Natural  Resources 
Cornell  University 
Ithaca,  New  York  14853 

4.  William  McFee 
Agronomy  Department 
Purdue  University 

West  Lafayette,  Indiana  47907 


5.  Robert  Johnson 
Earth  Resources 
Colorado  State  University 
Fort  Collins,  Colorado  80523 

6.  Sandy  McCarley 
Versar  Inc. 
P.O.  Box  1549 
Springfield,  Virginia  22151 

7.  Nancy  Dise 

Dept.  of  Environmental  Sciences 
University  of  Virginia 
Charlottesville,  Virginia  22903 

8.  David  Herzog 
1618  12th  Avenue 

San  Francisco,  California  94122 


This  report  should  be  cited  as 


Gibson,  J.H.,  J.N.  Galloway,  C.  Schofield,  W.  McFee,  R.  Johnson,  S.  McCarley, 
N.  Dise,  and  D.  Herzog.  1983.  Rocky  Mountain  Acidification  Study.  U.S.  Fish 
and  Wildlife  Service,  Division  of  Biological  Services,  Eastern  Energy  and  Land 
Use  Team,  FWS/OBS-80/40.17,  137  pp. 


EXECUTIVE  SUMMARY 


As  a  result  of  the  growing  concern  for  the  potential  effects  of  acid 
deposition  in  western  mountain  watersheds,  this  project  was  undertaken  to: 

1.  Determine  the  sensitivity  of  watersheds  characteristic  of  the  Rocky 
Mountain  Region  and  the  relationship  of  watershed  sensitivity  to 
geology  and  soils. 

2.  Evaluate  the  extent  of  current  acidification  and  the  potential  for 
increasing  acidification  with  increasing  deposition  of  nitrate  and 
sulfate. 

3.  Evaluate  the  results  of  the  above  in  terms  of  impacts  on  fish  popu- 
lations. 

4.  Develop  recommendations  for  assessment  of  future  trends  in  both 
changing  water  chemistry  and  impacts  on  fish  populations. 

Areas  were  selected  for  study  which  had  minimal  human  impact  and  for 
which  the  maximum  amount  of  data  on  soils,  geology  and  water  chemistry  already 
existed.  The  Rocky  Mountain  National  Park  (RMNP)  and  Yellowstone  National 
Park  (YNP)  areas  selected  exemplified  two  different  geologic  types.  The 
geology  of  these  areas  is  representative  of  a  large  portion  of  that  in  the 
total  Rocky  Mountain  region.  In  addition,  data  on  precipitation  chemistry 
were  available  from  the  National  Atmospheric  Deposition  Program  (NADP)  moni- 
toring studies  in  the  two  parks.  In  Yellowstone  National  Park  up  to  30  years 
of  water  chemistry  data  were  available,  along  with  information  on  geology  and 
soils.  In  Rocky  Mountain  National  Park  a  lake  and  stream  sampling  program  was 
conducted  to  collect  the  water  chemistry  information.  Throughout  the  project 
the  relationship  between  water  quality  data  and  fish  responses  determined  in 
research  studies  in  Scandinavia  and  eastern  North  America  were  relied  upon  to 
assess  potential  impacts  on  fish  populations. 

RESULTS 

The  two  parks  represent  a  contrast  in  geologic  materials;  Rocky  Mountain 
National  Park  being  primarily  underlain  by  granite  and  Yellowstone  National 
Park  by  volcanic  materials,  although  the  geologic  material  in  one  watershed  in 
Rocky  Mountain  National  Park  (the  Upper  Colorado  River  Basin)  is  of  volcanic 
origin.  The  examination  of  the  geochemistry  of  Rocky  Mountain  National  Park 
has  shown  that  many  areas  in  RMNP  are  sensitive  to  acidic  deposition  and  that 
this  sensitivity  is  primarily  determined  by  bedrock  geology.  In  addition, 
sensitivity  varies  inversely  with  elevation.  The  analyses  show  that  water- 
sheds underlain  by  granite  and  biotite  gneiss  and  schist  are  equally 
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sensitive  to  acidic  deposition.  The  lakes  and  streams  in  these  watersheds  had 
alkalinities  <200  jjeq/1 .  Waters  at  the  higher  elevations  (>3300  m)  were  very 
sensitive  (alkalinity  <  100  ueq/1).  The  Upper  Colorado  River  Basin  and  the 
Upper  Fall  River  Basin  contains  tertiary  intrusive  rocks  in  their  drainage, 
resulting  in  low  sensitivity  (alkalinity  >200  ueq/1).  Waters  in  Yellowstone 
National  Park  had  alkalinities  generally  above  200  ug/1 ,  with  a  few  <  200  ug/1 
on  rhyolite  or  basaltic  flows. 

The  analysis  of  sensitive  aquatic  systems  have  been  extrapolated  to  the 
Central  Rocky  Mountain  Region  by  delineating  areas  underlain  by  granite 
biotite  gneiss  and  schist  and  similar  gneisses  and  schists.  Areas  underlain 
by  these  formations  are  classified  as  sensitive  (alkalinity  <200  peq/1),  lakes 
and  streams  located  at  higher  elevations  (>3300  m)  can  be  classified  as  very 
sensitive  (alkalinity  <  100  ueq/1).  Areas  underlain  by  tertiary  intrusive 
rocks  were  classified  as  nonsensitive  (alkalinity  >200  ueq/1). 

The  Central  Rocky  Mountain  Region  is  currently  receiving  precipitation 
that  is  somewhat  acidic  (ave.  annual  pH  =  5.0,  NADP  1982).  The  analysis  of 
the  data  collected  in  RMNP  shows  that  little,  if  any,  acidification  of  lakes 
or  streams  has  occurred';  however,  areas  that  are  subject  to  periodic 
deposition  of  pollutants  during  upslope  air  movements  from  population  centers 
such  as  Denver  may  be  experiencing  some  acidification. 

As  stated,  high-elevation  lakes  and  streams  in  the  Central  Rocky  Mountain 
Region  are  very  sensitive  to  acidic  deposition.  Much  of  the  region  is 
underlain  by  rock  with  low  buffering  capacities  that  is  covered  by  highly 
permeable  soils  with  low  ion-exchange  capacities.  As  in  RMNP,  high-elevation 
lakes  and  streams  in  this  poorly  buffered  region  probably  will  become 
acidified  if  acidic  deposition  increases  to  the  level  currently  experienced  in 
the  northeastern  United  States.  An  evaluation  of  the  potential  impact  of 
increased  acidic  deposition  on  specific  lakes  and  streams  would  require  a 
drainage-by-drainage  assessment,  as  local  variability  in  bedrock,  hydrologic 
flow  path,  and  soil  development  may  have  an  overwhelming  influence  on 
sensitivity. 

With  respect  to  fish  population  there  is  currently  no  evidence  of  chronic 
acidification  and  thus  no  apparent  impact  on  fisheries.  However,  the  very  low 
base  cation  concentration  observed  in  the  headwater  drainages  of  Rocky 
Mountain  National  Park  suggests  extreme  sensitivity  to  acidification.  Fish 
populations  present  in  these  low  calcium  waters  may  be  particularly 
susceptible  to  osmoregulatory  stress  from  episodic  acidification.  The  few 
remaining  native  trout  (Salmo  clarki)  located  in  the  interior  regions  of  the 
Rocky  Mountains  persist  only  in  small  isolated  headwater  drainages.  The 
displacement  of  these  rare  and  endangered  genotypes  to  headwater  drainages 
also  makes  them  most  susceptible  to  potential  acidification  in  these  sensitive 
habitats.  Waters  in  volcanic  areas  such  as  Yellowstone  National  Park  are 
generally  of  high  alkalinity  and  thus  do  not  represent  potentially  sensitive 
habitats. 


RECOMMENDATIONS 

It  is  recommended  that  long-term  watershed  studies  be  established  at 
several  points  in  the  Rocky  Mountain  region.  Water  samples  should  be  analyzed 
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for  major  inorganic  constituents,  DOC,  aluminum,  alkalinity  and  total  acidity. 
In  addition,  surveys  of  water  chemistry  should  be  performed  in  other 
mountainous  areas,  not  only  in  the  Rocky  Mountain  region  but  also  in  other 
areas  of  the  western  United  States.  In  addition,  an  effort  should  be  made  to 
determine  the  rate  of  dry  deposition  of  neutral  salts,  i.e.,  calcium  sulfate. 

Behavioral  responses  and  immigration  tendencies  of  S.  clarki  populations 
exposed  to  episodes  of  acidification  should  be  determined  experimentally. 
Potential  for  aluminum  mobilization  in  the  Rocky  Mountain  watersheds  exposed 
to  increased  acid  deposition  should  be  studied.  Studies  of  winter-spring 
water  chemistry  in  headwater  catchments  should  also  include  evaluation  of  the 
movement  of  trout  populations  in  response  to  chemical  change.  The  potential 
sensitivity  to  acidification  of  watersheds  currently  occupied  by  endangered  or 
threatened  populations  of  S.  clarki  or  of  watersheds  considered  as  potential 
candidate  sites  for  reintroduction  should  be  determined  and  given  special 
consideration. 


CONTENTS 


Page 

EXECUTIVE  SUMMARY  iii 

FIGURES ix 

TABLES xi 

ACKNOWLEDGMENTS  xiii 

INTRODUCTION  1 

ROCKY  MOUNTAIN  NATIONAL  PARK  6 

Introduction  6 

Description  of  Study  Watersheds  6 

Lithologic  Units  9 

Sampling  Scheme  and  Methods  11 

Chemistry  of  Surface  Waters  11 

Soils  and  Surficial  Materials 16 

Results 18 

Soils  and  Surficial  Materials 18 

Chemistry  of  Surface  Waters  27 

Discussion 27 

Relationship  of  Surface  Water  Chemistry  to 
Atmospheric  Deposition,  Elevation  and  Soil 

and  Geology   27 

Alkalinity  (Sensitivity)  of  RMNP  Lakes  and  Streams  48 

Current  Acidification  Status  of  RMNP  48 

Future  Impacts  of  Acidic  Deposition   54 

YELLOWSTONE  NATIONAL  PARK  57 

Introduction  57 

Description  of  Geology   57 

Determination  of  Surface  Water  Chemistry  59 

Results  and  Discussion  61 

Lake  Alkalinities  (Sensitivities)  62 

Factors  Influencing  Alkalinity  (Sensitivity)  68 

Historical  Changes  in  Surface  Water  Chemistry  70 

Current  Status  and  Future  Trends  in  Surface  Water 

Chemistry 78 

FISH  POPULATIONS 82 

Introduction  82 

Assessment  Approach  82 

Results  and  Discussion  83 

Fish  Resources  in  Potentially  Sensitive  Areas  of  the 

Rocky  Mountains 83 

Fisheries  Management  Policies  in  Potentially  Sensitive 

Areas   86 


vn 


CONTENTS  (Concluded). 

Page 

Potential  Impacts  of  Acidification  on  Fish 

Populations 89 

CONCLUSIONS 95 

Introduction  95 

Rocky  Mountain  National  Park 95 

Soils  and  Geology 95 

Surface  Water  Chemistry  95 

Yellowstone  National  Park  97 

Effects  of  Current  and  Future  Acidification  Levels  in  Fish 

Populations  in  Yellowstone  and  Rocky  Mountain  National  Parks  ...  97 

Sensitivity  Evaluation  of  the  Central  Rocky  Mountain  Region  ...  98 

RECOMMENDATIONS  FOR  RESEARCH  AND  MANAGEMEMENT  102 

Introduction  102 

Rocky  Mountain  National  Park  and  Rocky  Mountain  Region  102 

Yellowstone  National  Park  103 

Fishery  Research  and  Management  103 

REFERENCES 105 

APPENDICES 

A.  Soil  and  Water  Chemistry  Sampling  Site  Maps,  Rocky 

Mountain  National  Park  113 

B.  Soil  Laboratory  Procedures,  Profile  Descriptions  and 

Chemical  Analyses  118 

C.  Water  Chemistry  Data,  Rocky  Mountain  National  Park  124 

D.  Water  Chemistry  Data,  Yellowstone  National  Park  133 


vm 


FIGURES 


Number  Page 

1    Location  of  study  areas,  Rocky  Mountain  National  Park, 

Colorado  7 


2 


Conceptual  behavior  of  CR,  CA  and  alkalinity 
concentrations  as  a  function  of  altitude  .  . 


12 


3  Distribution  of  pH  values  in  the  mineral  soil  and  the 

surface  organic  layer  of  all  four  watersheds  21 

4  "Excess"  base  cations  vs.  silica  for  two  sensitive 

watersheds  in  RMNP 31 

5  Graphs  of  CI  vs.  elevation  for  subbasins  in  RMNP 33 

6  Graphs  of  SO,  vs.  elevations  for  subbasins  in  RMNP  34 

7  Graphs  of  NO.,  vs.  elevation  for  subbasins  in  RMNP 35 

8  Graphs  of  alkalinity  vs.  elevation  for  subbasins  in 

RMNP 37 

9  Graphs  of  excess  base  cations  vs.  elevation  for 

subbasins  in  RMNP 38 

10  Graphs  of  silicate  vs.  elevation  for  subbasins 

in  RMNP 39 

11  Spatial  distribution  of  sensitivity  in  RMNP 49 

12  Status  of  23  lakes  in  RMNP 54 

13  Predictor  nomographs  projecting  the  best  and  worst-case 
estimates  of  the  effects  of  increased  acidic  deposition 

on  lakes  in  RMNP 56 

14  Regional-geological  map  of  Yellowstone  National  Park  ....  58 

15  Alkalinity  map,  Yellowstone  National  Park  63 

16  pH  map,  Yellowstone  National  Park 65 

17  Sulfate  map,  Yellowstone  National  Park  66 

ix 


FIGURES  (Concluded). 

Number  Page 

18  Location  of  hot  springs  and  geysers,  Yellowstone 

National  Park 67 

19  Alkalinity  versus  elevation,  lakes  draining  rhyolite 

bedrock,  Yellowstone  National  Park  70 

20  Historical  chemistry  comparisons,  High  and  Crescent 

Lakes 72 

21  Historical  chemistry  comparisons,  Grebe  and  Wolf  Lakes  ...     73 

22  Historical  chemistry  comparisons,  Ice  and  Crag  Lakes  ....     74 

23  Historical  chemistry  comparisons,  Cascade  Lake  75 

24  Chloride  versus  alkalinity,  Yellowstone  National  Park 

lakes 77 

25  Sulfate  versus  alkalinity,  Yellowstone  National  Park 

lakes 78 

26  Precipitation  pH  and  sulfate  concentration,  Yellowstone 

National  Park,  WY,  1980-1981  79 

27  Cross-hatching  shows  area  found  barren  of  fishes  by 
Jordan  in  1889,  with  the  exception  of  sculpins  in  the 

Gibbon  River  above  Gibbon  Falls  85 

28  Indigenous  distribution  of  Salmo  clarki  pleuriticus  91 

29  Indigenous  distribution  of  Salmo  clarki  s torn i as  92 

30  Spatial  distribution  of  sensitivity  in  the  central 

Rocky  Mountain  region:  Colorado  99 

31  Spatial  distribution  of  sensitivity  in  the  central 

Rocky  Mountain  region:  Wyoming  100 


TABLES 

Number  Page 

1  Major  bedrock  units:  Composition  and  relative 

susceptibility  to  chemical  weathering  10 

2  Analytical  techniques  15 

3  Representative  landform  types  (surface  features 

and/or  materials)  observed  in  research  area  17 

4  Seismic  refraction  data  summary  19 

5  Summary  of  soil  properties  at  four  locations  in  the 

East  Inlet  Valley 22 

6  Summary  of  soil  properties  at  nine  locations  in  the 

Upper  Colorado  drainage  23 

7  Summary  of  soil  properties  at  sample  locations  in 

the  Fall  River  drainage  basin 24 

8  Summary  of  soil  properties  at  sample  locations  in 

Glacier  Gorge  Basin  25 

9  Mean  ionic  concentrations  in  Rocky  Mountain  National 

Park  watersheds 28 

10  Comparison  of  mean  of  uncorrected  (U)  and  "excess"  (E) 

base  cations 30 

11  R-mode  varimax  factor  matrix  of  chemical  data  for  88 
lake  and  stream  samples  underlain  by  granite  and  biotite 

gneiss  and  schist  in  Rocky  Mountain  National  Park  41 

12  R-mode  varimax  factor  matrix  of  stream  chemistry,  %  soil 
organic  matter,  and  %  unreactive  rock  for  40  lake  and 
stream  samples  underlain  by  granite  and  biotite  gneiss 

and  schist  in  Rocky  Mountain  National  Park  43 

13  R-mode  varimax  factor  matrix  of  stream  chemistry 
for  33  lake  and  stream  samples  underlain  by  tertiary 

intrusive  bedrock  in  Rocky  Mountain  National  Park  44 


XI 


TABLES  (Concluded). 

Number  Page 

14  Discriminant  analysis  of  surface  water  chemistry 
from  117  lake  and  stream  samples  from  Rocky  Mountain 

National  Park 46 

15  Reclassification  results  46 

16  Comparison  of  Rocky  Mountain  National  Park  and 

Hubbard  Brook  precipitation  chemistry  51 

17  Calculation  of  acidification  status  for  Rocky  Mountain 
National  Park  watersheds  according  to  ion  balance 
considerations  52 

18  Analytical  methods,  Yellowstone  National  Park  61 

19  Regional  alkalinity  of  Yellowstone  National  Park 

lakes 69 

20  Relative  concentration  of  alkalinity  and  sulfate, 

lakes  with  historical  data 76 

21  Sulfate  levels  in  selected  acidified  waters  80 

22  Current  and  historical  fish  population  status  of 

lakes  and  streams  in  Yellowstone  National  Park 86 

23  Total  alkalinity  and  fish  population  status  for 

sensitive  lakes  in  Yellowstone  National  Park  87 

24  Endangered  and  threatened  Salmo  genotypes  in  the 

Rocky  Mountain  Region 89 


XI  i 


ACKNOWLEDGMENTS 


As  in  any  project  as  extensive  as  the  Rocky  Mountain  Acidification  Study, 
there  are  a  large  number  of  individuals  who  make  invaluable  contributions  who 
are  not  credited  in  the  author  list.  Such  contributions  were  critical  to  the 
success  of  this  study. 

In  particular,  we  wish  to  express  our  appreciation  to  Jill  Baron  of  the 
National  Park  Service  Water  Resources  Field  Support  Laboratory  for  her 
generous  assistance  in  the  field  sampling  program  at  Rocky  Mountain  National 
Park  (RMNP),  and  Robert  E:  Gresswell,  fishery  biologist  with  the  U.S.  Fish  and 
Wildlife  Service,  Yellowstone  Fisheries  Assistance  Office  at  Yellowstone 
National  Park  (YNP),  for  his  assistance  in  data  gathering  and  advice 
throughout  the  project.  Without  his  contributions  we  would  not  have  been  able 
to  evaluate  the  YNP  data. 

Another  important  contribution  to  the  data  collection  was  made  by  William 
Locke  of  Montana  State  University,  who  provided  information  on  soils  on  RMNP 
and  aided  in  the  selection  of  sampling  sites.  We  also  wish  to  thank  Dave 
Stevens  of  RMNP  and  Wayne  Hamilton  of  YNP,  who  facilitated  our  access  to  the 
parks  and  to  the  available  information  on  soils,  geology  and  water  chemistry. 

We  also  wish  to  acknowledge  the  support  for  the  project  from  the  U.S. 
Fish  and  Wildlife  Service  and  Kent  Schreiber,  who  provided  the  major  funding, 
and  the  National  Park  Service  Water  Resources  Field  Support  Laboratory,  and 
Raymond  Herrmann,  who  provided  support  for  one  of  the  graduate  students. 

Lastly,  we  wish  to  express  our  sincere  appreciation  to  Janice  Hill  of  the 
Natural  Resource  Ecology  Laboratory  Publications  Section  and  her  staff,  Katie 
Curry,  Kay  McElwain  and  Susan  Taylor  for  their  efforts  "beyond  the  call  of 
duty  in  organizing,  typing  and  editing  this  report. 


XI  11 


INTRODUCTION 

Acid  precipitation  and  its  subsequent  impact  on  aquatic  and  terrestrial 
ecosystems  is  a  matter  of  extensive  research  and  debate  in  much  of  the 
industrialized  world  in  the  northern  hemisphere.  Northern  Europe,  United 
States,  Canada,  and  Japan,  have  all  experienced  an  increase  in  the  acidity  of 
precipitation  and  in  many  of  these  areas  this  increase  has  been  associated 
with  deterioration  in  aquatic  and  terrestrial  ecosystems.  The  extent  of  this 
deterioration  and  its  relationship  to  acidic  precipitation  has  been  the 
subject  of  considerable  controversy.  However,  it  appears  clear  at  this  time 
that  in  the  Scandanavian  countries,  and  limited  areas  of  northeastern  U.S.  and 
eastern  Canada,  lake  and  stream  acidification  as  a  result  of  acidic  precipi- 
tation has  been  demonstrated  beyond  reasonable  doubt.  In  addition,  in  many  of 
these  areas  there  have  been  documented  declines  in  fish  populations  associated 
with  the  increase  in  acidity.  Effects  on  terrestrial  systems  and  their 
association  with  acidic  precipitation  are  less  well  established,  but  there  do 
appear  to  be  declines  in  forest  productivity  in  areas  of  Europe  and  north- 
eastern United  States  that  are  correlated  with  atmospheric  deposition.  These 
observations  have  spurred  initiation  of  extensive  research  programs  to 
determine  the  relationship  between  acidic  deposition  and  declines  in  the 
productivity  of  terrestrial  and  aquatic  systems  (National  Research  Council 
1981). 

The  emphasis  of  research  to  date  has  been  in  eastern  North  America,  areas 
of  which  now  receive  rainfall  with  an  average  pH  of  4.1  or  below  along  with 
elevated  levels  of  sulfate  and  nitrate  (NADP  1980).  In  addition,  areas  in 
which  the  most  significant  impacts  of  acidic  deposition  have  been  demonstrated 
in  aquatic  systems  are  underlain  by  geologic  material  of  low  buffering 
capacity.  Such  areas  are  found  in  eastern  Canada  and  northeastern  United 
States.  The  combination  of  limited  buffering  capacity  of  soils  and  parent 
materials  and  elevated  levels  of  acid  deposition  have  led  to  lake  and  stream 
acidification,  with  subsequent  effects  on  the  biotic  community.   The  effects 


on  aquatic  systems,  particularly  the  decline  in  fish  populations,  has  been  the 
most  significant  factor  in  encouraging  governments  and  industry  that  we  are 
dealing  with  a  serious  problem  which  requires  not  only  increased  levels  of 
research,  but  also  the  development  of  policy  to  bring  about  reductions  in 
anthropogenic  emissions  of  sulfur  and  nitrogen  oxides. 

While  the  emphasis  in  the  eastern  U.S.  and  Canada  is  understandable 
because  of  the  documented  problems,  in  terms  of  both  levels  of  acidic 
deposition  and  impacts,  it  has  been  noted  that  the  same  conditions  of 
sensitivity  exist  in  areas  of  the  western  United  States  including  the  Rocky 
Mountains,  the  Sierras,  Cascades,  and  other  western  ranges.  These  areas  are 
primarily  of  granitic  bedrock,  and  the  limited  data  available  has  shown  that 
the  high  elevation  lakes  and  streams  have  sensitivities  comparable  to  those  in 
the  eastern  United  States  and  Canada.  On  the  other  hand,  data  from  the 
National  Atmospheric  Deposition  Program  (NADP)  and  other  research  has  not 
demonstrated  that  these  areas  are  receiving  precipitation  with  acidity 
comparable  to  that  of  eastern  North  America.  Preliminary  data  from  NADP  would 
indicate  that  the  pH  of  rainfall  in  the  Rocky  Mountain  West  is  averaging  4.8 
to  5.0,  with  a  possibility  that  in  some  areas  at  higher  elevations  the  pH  may 
be  somewhat  below  these  values.  Precipitation  with  pH  values  below  5  would 
generally  be  considered  to  be  more  acidic  than  normal.  This  raises  the 
question  whether  the  most  sensitive  lakes  and  streams  at  the  higher  elevations 
are  being  acidified.  Because  of  the  importance  of  these  lakes  as  trout 
fisheries  and  their  aesthetic  value  to  millions  of  visitors,  it  is  important 
to  gain  a  better  understanding  of  the  current  status  of  these  lakes  with 
respect  to  buffering  capacity  and  to  assess  the  potential  for  acidification 
and  subsequent  impact  on  fish  populations.  A  survey  of  the  literature 
indicates  however,  that  little  water  chemistry  data  are  available  for  high 
lakes  (>2500  m)  in  the  central  Rocky  Mountains. 

The  most  extensive  watershed  study  in  this  region  has  been  at  the  Como 
Creek  watershed  in  the  northern  front  range  of  the  Colorado  Rockies  (Lewis  and 
Grant  1979).  At  least  150  weeks  of  consecutive  stream  and  precipitation 
chemistry  data  are  available  (Lewis  and  Grant  1979;  Lewis  1982),  and  the  first 
report  of  acidic  deposition  in  the  Rockies  was  made  at  this  site  (Lewis  and 
Grant  1980).  The  alkalinity  in  the  waters  of  Como  Creek  averaged  192  ueq/1  as 
determined  by  potentiometric  titration  to  pH  4  (Lewis  and  Grant  1979).  Other 
studies  of  high  mountain  lakes  include  those  by  Dodson  (1981)  in  the  Elk 


Mountains  of  westcentral  Colorado.  This  area,  located  on  the  west  slope  of 
the  Continental  Divide,  is  underlain  by  quartzite,  siltstone,  sandstone,  and 
shale  (Harte  et  al.  1983).  Dodson  (1982)  found  these  lakes  to  have  low 
alkalinity,  8  to  350  (jeq/1.  A  recent  study  by  Harte  et  al.  (submitted)  in  the 
same  area  reveals  alkalinities  ranging  from  8  to  250  ueq/1  with  pH  values  of 
5.9  to  7.88.  Turk  and  Adams  (1983)  have  carried  out  a  study  of  the  chemistry 
of  high-elevation  lakes  in  the  Flat  Tops  Wilderness  area  of  western  Colorado. 
The  bedrock  in  this  area  is  predominantly  basalt  with  some  granitic  outcrops. 
Alkalinities  in  this  region  range  from  70  to  1400  ueq/1.  In  addition,  data 
from  a  lake  survey  of  the  Mount  Zirkle  Wilderness  area  by  Turk  (unpublished) 
shows  alkalinities  ranging  from  12  to  315  ueq/1.  The  bedrock  in  this  area  is 
primarily  granite.  In  his  study  of  the  South  St.  Vrain  Watershed  in  the  Front 
Range  of  Colorado,  Thurmon  (in  press)  found  that  alkalinities  in  the  head- 
waters averaged  82  ueq/1.  This  area  is  underlain  by  silverplume  granite  and 
biotite,  gneiss,  and  schist. 

These  studies  of  water  chemistry  of  headwater  lakes  and  streams  in  the 
Central  Rocky  Mountain  Region  have  only  recently  been  conducted.  The  lack  of 
historical  data  prevents  the  evaluation  of  any  long-term  trend  in  acidity. 
Lewis  and  Grant  (1979)  reported  an  increase  in  stream  hydrogen  ion  concen- 
trations in  the  Como  Creek  watershed  over  a  150  week  period.  This  record, 
however,  is  too  short  to  be  considered  evidence  of  increasing  acidification. 
Lewis  (1982)  reported  an  average  decrease  in  alkalinity  of  180  ueq/1  for  104 
lower  elevation  lakes  and  streams  surveyed  by  Pennak  in  1938-1942,  and 
resurveyed  by  Lewis  in  1979.  Although  Lewis  attributed  this  loss  to  acidic 
deposition,  the  concurrent  increase  in  total  dissolved  solids  (TDS)  suggests 
that  hydrologic  variability  is  responsible  for  this  decrease  in  alkalinity. 
In  addition,  a  180  ueq/1  loss  of  alkalinity  is  much  greater  than  any 
alkalinity  loss  observed  in  strongly  acidified  lakes  in  the  northeastern 
United  States  (Hendrey  et  al.  1980).  Thus,  no  convincing  evidence  of 
acidification  in  the  Central  Rockies  has  yet  been  presented,  although  several 
areas  have  been  shown  to  be  potentially  sensitive  (alkalinities  less  than  200 
ueq/1)  to  acidic  deposition.  Also,  none  of  the  studies  in  the  Rocky  Mountain 
Region  to  date  have  associated  the  current  alkalinity  and  pH  of  high-elevation 
lakes  and  streams  with  potential  changes  in  the  future  and  subsequent  effects 
on  fish  populations. 


In  general,  the  objective  of  this  project  was  to  provide  an  assessment  of 
lake  and  stream  sensitivity  of  selected  areas  in  the  Rocky  Mountain  region  and 
to  relate  this  to  potential  effects  on  fish  populations.  It  was  also  felt 
that  by  coupling  water  chemistry  data  with  data  on  the  geology  and  soils  in 
the  sampled  watersheds,  that  a  relationship  could  be  established  between  water 
chemistry  and  geology  and  soils  types  and  that  this  relationship  could  be 
extrapolated  to  other  areas  of  the  Rocky  Mountain  Region.  More  specifically 
the  objectives  were  to: 

1.  Determine  the  sensitivity  of  watersheds  characteristic  of  the  Rocky 
Mountain  region  and  the  relationship  of  watershed  sensitivity  to 
geology  and  soils. 

2.  Evaluate  the  extent  of  current  acidification  and  the  potential  for 
increasing  acidification  with  increasing  deposition  levels  of 
nitrate  and  sulfate. 

3.  Evaluate  the  results  of  the  above  in  terms  of  effects  on  fish  popu- 
lations. 

4.  Develop  recommendations  for  assessment  of  future  trends  in  both 
changing  water  chemistry  and  effects  on  fish  populations. 

The  plan  called  for  selecting  areas  for  the  study  which  had  minimal  human 
impact  and  for  which  the  maximum  amount  of  data  on  soils,  geology,  and  water 
chemistry  might  already  exist.  The  plan  relied  heavily  on  the  established 
relationship  between  water  quality  data  and  fish  responses  as  determined  in 
research  studies  in  Scandinavia  and  eastern  North  America.  After  considering 
potential  study  sites  available,  two  areas  of  different  geologic  type  were 
selected:  Yellowstone  National  Park  and  Rocky  Mountain  National  Park.  The 
geology  of  these  areas  is  representative  of  a  large  portion  of  that  in  the 
total  Rocky  Mountain  Region.  Another  factor  considered  was  that  each  of  these 
parks  participated  in  the  National  Atmospheric  Deposition  Program  (NADP) 
precipitation  chemistry  monitoring  network,  and  therefore  had  available  the 
necessary  precipitation  chemistry  data  to  be  used  in  the  analysis  of  the 
extent  of  acidification  of  the  surface  waters.  The  existence  of  information 
on  geology  and  soils  in  both  parks,  while  incomplete,  proved  of  considerable 
value.  In  Yellowstone  National  Park,  water  chemistry  data  was  available  for 
up  to  thirty  years  in  some  lakes.  A  single  season  of  lake  and  stream  sampling 
was  conducted  in  Rocky  Mountain  National  Park  to  determine  water  chemistry. 


The  first  sections  of  the  report  begins  with  an  evaluation  of  the 
sensitivity  of  lakes  and  streams  in  Rocky  Mountain  National  Park  and 
Yellowstone  National  Park.  Sensitivity  (alkalinity)  is  related  to  both  soil 
and  bedrock  influences,  as  well  as  elevational  gradients.  The  following 
section  discusses  the  relationship  of  the  water  chemistry  data  to  fish 
populations  based  on  currently  available  fisheries  data.  Conclusions  and 
recommendations  are  presented  in  the  final  two  sections  of  the  report. 


ROCKY  MOUNTAIN  NATIONAL  PARK 

INTRODUCTION 

A  discussion  of  the  studies  performed  in  Rocky  Mountain  National  Park 
(RMNP)  is  presented  in  this  section  of  the  Rocky  Mountain  Acidification  Study 
(RMAS)  report.  Rocky  Mountain  National  Park  is  located  in  northcentral 
Colorado  along  the  Continental  Divide.  It  has  an  area  of  106,700  ha  and 
ranges  in  elevation  from  2,329  m  to  4,345  m.  Approximately  one-third  of  the 
area  and  the  107  named  mountains  are  over  3,353  m  high.  The  mountain  building 
was  accomplished  by  a  regional  uplift  of  Precambrian  igneous  and  metamorphic 
rocks  during  Late  Mesozoic  and  Early  Cenozoic  time  (Richmond  1974).  Volcanic 
activity  has  more  recently  altered  the  western  side  of  the  park.  Alpine 
glaciation  during  the  Pleistocene  formed  characteristic  U-shaped  valleys, 
steep-sided  cliffs,  and  lateral  moraines.  Cirque  lakes  are  common.  These 
geologic  characteristics  made  this  park  an  ideal  site  in  which  to  evaluate 
lake  and  stream  sensitivity. 

DESCRIPTION  OF  STUDY  WATERSHEDS 

Four  watersheds  were  selected  within  the  park  boundaries  in  which  to 
conduct  the  studies.  The  geologic  criteria  used  in  selecting  watersheds  were 
geologic  control  (bedrock  and  glacial  erosion/deposition)  as  determined  from 
existing  information,  with  attention  being  given  to  the  representativeness  of 
selected  watersheds  to  the  Rocky  Mountain  region  in  general  and  RMNP  in 
particular.  Climatic  factors  suggested  the  selection  of  watersheds  on  both 
the  east  and  west  side  of  the  Continental  Divide.  Access  to  these  watersheds 
by  foot  trail  was  also  a  consideration. 

The  specific  watersheds  selected  for  study  were  Fall  River  Basin,  Glacier 
Gorge  Basin,  East  Inlet  Basin  and  Upper  Colorado  River  Basin  (Figure  1). 
Information  on  the  bedrock  geologic  control  of  these  watersheds  is  available 
with  the  bulk  of  the  data  compiled  on  a  1:48000  scale  map  (unpublished)  by  W. 


UPPER   COLORADO 
RIVER    BASIN 


Figure  1.  Location  of  study  areas,  Rocky  Mountain  National  Park,  Colorado. 
Subbasins  are  (A)  Roaring  River,  (B)  Ypsilon  Creek,  (C)  Upper  Fall  River,  (D) 
Tyndall  Gorge,  (E)  Loch  Vale,  (F)  Glacier  Creek. 


A.  Braddock,  University  of  Colorado,  Boulder.   There  are  also  two  pertinent, 
unpublished  theses  (Abbott  1974;  Shroba  1977). 

The  geology  in  RMNP  is  similar  to  about  70%  of  the  mountainous  areas 
within  the  Rocky  Mountain  region.  The  most  extensive  geologic  material 
present  in  the  Colorado  Rocky  Mountain  region  is  Precambrian  granites  and 
metamorphic  rocks.  These  rock  types  are  present  in  three  of  the  four  basins 
selected  for  this  study:  Fall  River,  Glacier  Gorge,  and  East  Inlet  Basins. 
The  fourth  watershed,  the  Upper  Colorado  River  Basin  has  as  its  major  rock 


types  Tertiary  vol  cam' cs  with  intrusives  of  andesitic  to  rhyolitic  composi- 
tion. Glacial  till  also  constitutes  a  major  geologic  material  in  all  four  of 
the  watersheds,  but  especially  in  the  Fall  River  Basin. 

Fall  River,  Glacier  Gorge,  and  East  Inlet  Basins  have  tributaries 
beginning  at  the  Continental  Divide  at  elevations  above  timberline.  The 
tributaries  start  in  material  that  is  predominantly  rock  talus  with  some 
tundra  vegetation  present.  Soil  development  above  timberline  is  generally 
weak  thus  making  geologic  control  the  most  significant  factor  controlling 
water  chemistry.  Most  of  the  tributaries  have  as  their  source  some  type  of 
snow  field,  glacier  or  cirque  lake. 

The  Upper  Colorado  River  Basin  watershed,  with  its  bedrock  of  inter- 
mediate to  rhyolitic  composition,  intrusive  and  extrusive  igneous  rocks,  has 
approximately  analagous  geology  to  that  present  in  the  San  Juan  mountains  of 
southwest  Colorado  and  the  northwest  section  of  Yellowstone  National  Park. 
This  area  is  the  only  section  of  RMNP  where  these  geologic  conditions  exist. 
These  rock  types  contribute  more  in  terms  of  buffering  capacity  to  the  soils 
because  of  mineralogy  and  more  rapid  weathering. 

The  third  broad  type  of  geologic  material,  glacial  deposits,  is  present 
in  all  four  basins  with  Fall  River  having  the  largest  percentage.  Glacial 
deposits  range  from  late  moraines  of  up  to  200  feet  thick  to  thin  veneers  in 
upstream  regions.  These  deposits  have  a  major  bearing  on  the  soils  encoun- 
tered as  well  as  the  general  water  quality.  The  morainal  material  in  general 
is  coarse  grained  with  abundant  boulders,  gravels  and  sands.  This  material  is 
well  drained  and  thus  water  interacts  not  only  with  the  soils  developed  on  the 
glacial  deposits  but  also  with  the  deposits  themselves. 

Another  factor,  climatic  regimes,  has  profound  influences  on  watershed 
characteristics.  It  is  noted  that  two  watersheds,  Fall  River  and  Glacier 
Gorge  basins,  are  located  on  the  east  side  of  the  continental  divide,  while 
the  other  two  watersheds,  East  Inlet  and  Upper  Colorado  River  Basins  are  on 
the  west  side  of  the  divide.  This  provides  a  contrast  in  precipitation  levels 
and  vegetation  types,  both  major  factors  in  soil  development. 

A  bedrock  and  surficial  geology  map  using  a  U.S.  Geological  Survey  7.5 
minute  quadrant  topographic  map  was  prepared  for  each  watershed  using 
primarily  existing  data  (Braddock,  unpublished;  Cole  1977;  Abbott  1974). 
Also,  a  map  using  the  same  base  has  been  prepared  showing  the  geomorphic 
provinces  in  each  watershed.   The  map  of  the  geomorphic  provinces  (talus 
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slopes,  scoured  bedrock,  moraine  veneer,  etc.),  shows  the  prime  source  of  the 
soils  (parent  material)  that  have  developed  in  those  provinces  along  with  the 
primary  vegetation  present  (Locke,  unpublished). 

The  four  watersheds  were  subdivided  on  the  basis  of  tributary  drainage  as 
follows  (Figure  1): 

Glacier  Gorge:        Tyndall  Gorge,  Loch  Vale,  Glacier  Creek 

subbasins 

Fall  River:  Roaring  River,  Ypsilon  Lake,  Upper  Fall 

River  subbasins 

East  Inlet:  Not  subdivided 

Upper  Colorado  River:  Box  Canyon  subbasin 

Lithologic  Units 

Igneous  and  metamorphic  rocks  are  present  in  the  four  watersheds  and 
range  in  age  from  Precambrian  to  Tertiary;  the  metamorphic  rocks  are  confined 
to  the  Precambrian.  Unconsolidated  materials  of  Quaternary  age  are  abundant 
in  all  basins.  The  major  bedrock  units  exposed  in  the  four  watersheds  are 
(Abbott  1974;  Cole  1977;  O'Neill  1981): 

Xqs:  Biotite  gneiss  and  schist 

Xam:  Amphibolite 

Xgg:  Quartz  diorite  gneiss 

be:  Boulder  Creek  granodiorite 

Ysp:  Silver  Plume  granite 

PEa,  PEab,  PEap:  Andesite  flows,  breccia  flows,  prophyry 

Ngd:  Granodiorite  of  Mt.  Richtofen  stock 

Ngr:  Granite  of  Mt.  Cumulus  stock 

Nvr:  Ash  flow  tuff 

The  mineralogical  compositions  of  bedrock  geologic  units  are  summarized 
in  Table  1.  Using  the  Goldich  (1938)  weathering  sequence  and  the  modal 
percentages  of  the  lithologic  units,  a  suggested  weathering  stability 
classification  is  shown.  The  susceptibility  to  chemical  weathering  increases 
to  the  right.  The  classification  is  essentially  based  on  the  relative  modal 


<D 

> 

•i — 

-!-> 

03 

r— 

<D 

S- 

T3 

CZ 

03 

c 

o 

•r— 

+-> 

•( — 

CO 

o 

en 

Q. 

cz 

E 

•i— 

O 

S- 

CD 

cu 

-C 

+-> 

.  . 

03 

in 

CD 

+-> 

5 

•  r— 

cz 

l 

Z3 

fO 

a 

j*. 

•i — 

u 

£ 

o 

CD 

S- 

x: 

-o 

u 

a; 

jd 

o 

+-> 

s- 

o 

>> 

03  .i- 


r-l  +-> 

a. 

CD  QJ 

r—  U 

_Q  CO 

03  Z3 

I—  CO 


E 

JD 

i 

TO 

03 

I*-. 

X 

0 

s- 
<v 

-CZ 

+-> 

03 

Q. 

a» 

03 

? 

LU 

o_ 

CO 

»> 

o 

JD 

s: 

a> 

03 

c 

03 

LU 


CO 

cr 
x 


CO 

-l-> 

•  1 — 

CD 

cz 

CD 

Z3 

X 

u 

•  ^ 

CD 

o 

1 — 

T3 

o 

CD 

sz 

z 

+-> 

ra 


5- 
> 


03 
CD 


CD 


CD  O  O  CD 

^"  S-  CNJ  CD  ro 

I  -Q  I  I  I 

CD  03  1^.  O  CD 

CSJ  —I  CSJ 


o 


O         -i- 


s_ 

o 

c 


CD 


■a       cr> 

CZ  -i- 


m 


CD 

CSJ  T- 

CS1  r— 

o 


00 


en 


LO 

cr> 

CD 

o 
m 

LO 

m 

1 

CD 
CO 

1 

1 

o 

1 

LO 
I— 1 

CD 


CM 
CO 


CD 


CD 


S- 

o 


-a 

cz 


00 


S- 
o 


u 

i— 1 

C£> 

LO 

■o 

LO 

_Q 

CD 

-Q 
03 

S- 

co 

cn 

CSJ 

CSJ 

cz 

I— 1 

Q. 

sz 

CD 

CO 

+-> 

r^ 

00 

"3" 

•< — 

r^. 

>- 

03 
CO 

+-> 

CO 

CSJ 

OO 

CSI 

o 

i- 

03 

CD 

CO 

LO 

CD 

-3" 

>— i 

-z. 

1 

cn 

CJD 

CSJ 


CSJ 


-O. 

a> 

a. 

>> 

+-> 

ai 

CO 

QJ 

ai 

CO 

CO 

CO 

CD 

CZ 

CO 

03 

03 

■a 

CD 

-!-> 

>> 

■  i — 

03 

1 — 

i — 

i 

cz 

+-> 

•  i — 

S- 

i — 

1 — 

u 

u 

CD 

OJ 

•r— 

s_ 

o 

M 

u 

<_) 

o 

o 

-!-> 

1 — 

-t-> 

C0 

CO 

-t-> 

o 

o 

•  1 — 

•  ) — 

•  i — 

-Q 

OJ 

•) — 

CO 

S- 

£_ 

-C 

CD 

CD 

+-> 

CZ 

cz 

T3 

CD 

03 

u 

+-> 

03 

03 

O 

S- 

CD 

s_ 

U 

Z3 

•  1 — 

S- 

r— 

i — 

•  i — 

o 

03 

o 

U 

cr 

s: 

o 

Q. 

Q_ 

CCt 

:n 

s: 

CD 

<C 

CO 

a. 
+-> 

U 

o 

S- 


5 
o 

cr 
o 

-o 

CD 

CO 

03 

JD 

CD 


-t-> 
03 
OH 


cu 


CO 

CO 

CO 

o 

-t-> 

Q. 

•r— 

S- 

-(-> 

o 

O 

■a 

cz 

03 

s- 

CO 

JD 

03 

03 

5 

1 

CO 

II 

cz 

03 

. 

CD 

S- 

E  -C 

03 

i — 

1 

03 

U 

~ 

•  i — 

CD 

-t-J 

CZ 

CL 

•  f— 

O 

CO 

CD 

>>-a 

.a 

CZ 

«=C 

CD 

CD  II 

03 

+-> 

CZ 

■o 

CD 

CZ 

u  <c 

t- 

C0 

~ 

a.  ai 

CO 

C(- 

03 

o 

( — 

u 

c 

o 

o 

CD 

•1 — 

•  i — 

•*-> 

i — 

03  O 

CZ 

•  1— 

II 

E 

S- 

CD 

CD 

+-> 

t — 

CD 

CD 

O 

03  JD 

10 


percentages  of  mafic  minerals,  hornblende  and  biotite,  and  plagioclase,  as 
these  minerals  are  more  susceptible  to  chemical  weathering  (Goldich  1938). 
Some  attention  has  been  given  to  the  physical  stability  of  the  rock,  i.e., 
schistosity  in  the  biotite  gneiss  and  schist  (Xqs).  The  location  of  the  ash 
flow  tuff  (Nvr)  in  the  classification  is  questionable  because  of  possible 
influence  of  the  moderately  high  temperature  of  formation  and  the  glassy 
matrix  on  weathering. 

Of  the  four  basins  selected  for  study  (two  on  each  side  of  the 
Continental  Divide),  three  are  dominated  by  metamorphic  rocks  and  soils 
derived  from  those  rocks.  The  third  is  in  a  region  of  volcanic  rocks  typical 
of  a  significant  section  of  the  Rocky  Mountain  region.  The  three  basins 
representing  metamorphic  materials  differ  in  the  percentage  of  glacial 
deposits  in  the  basin  vs.  the  steep  slope  areas  with  either  thin  or  no  soil. 

SAMPLING  SCHEME  AND  METHODS 

Chemistry  of  Surface  Waters 

Conceptual  basis.  The  general  objectives  of  the  RMAS  were  to  determine 
the  sensitivity  of  waters  in  the  Rocky  Mountain  region  to  acidification  by 
acid  precipitation  and  to  determine  whether  this  will  have  an  impact  on  fish 
populations.  An  earlier  section  of  this  report  outlined  the  strategy  for 
selecting  the  sampling  sites  and  described  how  existing  geologic  and  soils 
information  had  been  used  to  select  systems  and  subsystems  in  the  Park  that 
are  representative  of  the  Rocky  Mountain  region.  Lakes  and  streams  in  nine 
large  watersheds  (four  basins)  were  sampled  on  an  elevational  gradient.  For 
each  subsystem,  it  was  expected  there  would  be  an  elevational  trend  in  the 
concentrations  of  the  base  cations  (calcium,  magnesium,  sodium,  and  potassium) 
because  these  ions  are  derived  from  terrestrial  ecosystems  by  the  process  of 
soil  or  parent  material  weathering.  An  elevational  gradient  in  the  concen- 
trations of  the  acid  anions  (sulfate,  chloride,  and  nitrate)  was  not  expected. 
This  concept  is  illustrated  in  Figure  2,  a  theoretical  illustration  in  which 
concentration  of  the  sum  of  base  cations  (CR)  and  the  sum  of  the  acid  anions 
(C.)  is  plotted  as  a  function  of  elevation.  At  the  point  where  they 
intersect,  the  system  has  lost  all  alkalinity  and  is  about  to  develop  strong 
acidity.  In  the  lower  elevation  areas  where  CR  is  larger  than  C*,  there  has 
to  be  another  anion  for  charge-balance  considerations.  In  these  systems,  the 
anion  is  bicarbonate,  which  is  equated  to  alkalinity  because  of  the  absence  of 
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Figure  2.  Conceptual  behavior  of  Cn,  CA  and  alkalinity  concentrations  as  a 
function  of  altitude. 

Al  and  organic  compounds  that  could  contribute  to  the  alkalinity.  For  the 
higher  elevations  where  C.  are  larger  than  Cn,  there  have  to  be  other  cations, 
again  for  reasons  of  charge  balance.   The  cations  are  hydrogen  and  aluminum. 

By  definition  this  region  has  already  developed  strong  acidity.  The  effect  of 

-2    -1      -1 
increased  atmospheric  deposition  of  the  acid  anions  SO,  ,  NCL  and  CI   will 

cause  the  horizontal  line  of  C.  to  move  upward  resulting  in  a  larger  number  of 
lakes  at  the  lower  elevations  developing  strong  acidity.  This  suggests  that 
the  most  sensitive  (lowest  alkalinity)  lakes  will  be  at  the  watershed  head- 
waters. 

The  dotted  line  in  Figure  2  is  worthy  of  mention.  At  lower  elevations 
where  soil  development  is  more  extensive,  the  process  of  sulfate  adsorption 
may  remove  atmospherically  derived  sulfur  from  the  waters  of  the  watershed. 
One  of  the  goals  of  this  project  was  to  determine  the  degree  of  sulfate 
adsorption  as  a  function  of  elevation.  Therefore,  if  the  elevational  gradient 
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for  the  acid  anions,  specifically  sulfate,  follows  the  dotted  line,  we  will  be 
able  to  field  test  this  by  the  determination  of  sulfate  adsorption  capacity. 

This  concept,  relating  the  concentrations  of  various  species  to  eleva- 
tion, was  the  basis  for  devising  the  general  water  chemistry  sampling  plan. 

Sampling  plan.  Water  chemistry  sampling  within  each  sub-basin  was 
carried  out  in  a  one  to  four  day  period  to  reduce  variation  in  hydrological 
conditions.  Samples  were  collected  under  base  flow  conditions,  i.e.,  sampling 
did  not  occur  within  24  hours  after  rainstorms.  Lake  samples  were  collected 
at  each  lake  inlet,  outlet,  and  center  location.  Stream  samples  were  taken  25 
meters  below  each  confluence  and  at  approximately  150  meter  elevation 
intervals.  (Sample  location  maps  are  presented  in  Appendix  A.)  Stream 
samples  were  collected  in  the  middle  of  the  stream  under  falling  water,  while 
lake  samples  were  taken  0.3  meters  below  the  surface  at  the  center  of  the 
lake. 

Two  250  ml  samples  were  collected  at  each  site  in  clean  polyethylene 
bottles.  Each  bottle  had  been  washed  with  hot  water  and  detergent,  rinsed 
five  times  with  hot  tap  water  and  five  times  with  deionized  water.  The 
analysis  of  blank  samples  indicated  no  contamination  from  this  cleaning 
procedure.  Each  bottle  was  rinsed  three  times  in  situ  with  sample  water 
before  filling.  One  sample  was  immediately  preserved  with  reagent  grade 
chloroform  for  later  base  cation  and  acid  anion  analysis.  The  sample 
identification  number,  location,  date,  and  time  of  collection  were  printed  on 
the  bottle  and  recorded  in  the  field  log  book. 

Conductivity  and  pH  were  measured  in  a  separate  aliquot  of  sample  at  the 
site.  Unpreserved  250  ml  samples  for  alkalinity  analysis  were  stored  at  5°C 
for  one  to  two  weeks  before  analysis.  The  preserved  samples  were  analyzed  for 
chloride,  nitrate  +  nitrite,  sulfate,  phosphate,  silicate,  calcium,  magnesium, 
sodium,  potassium,  and  ammonium  within  two  to  three  months  after  sample 
collection.  The  pH  of  the  stream  waters  was  too  high  to  include  aluminum 
analysis. 

Field  analyses.  Conductivities  were  measured  in  the  field  using  a 
Beckman  model  RC-16C  meter  and  a  Yellow  Spring  conductivity  cell  (YSC  model 
3404).  Corrections  to  25°C  were  calculated  according  to  Standard  Methods 
(A.P.H.A.  1976).  Field  pH  measurements  were  taken  with  a  Corning  digital  pH 
meter  and  a  Corning  model  476182  pH  electrode.  The  meter  was  periodically 
calibrated  with  cold  dilute,  strong  acids  in  the  laboratory.  Each  measurement 
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was  preceded  by  a  two  point  calibration  using  pH  7.00  and  pH  4.00  buffer 
solutions.  Corrections  for  the  temperature  difference  between  the  sample  and 
buffer  were  made  using  meter  adjustments.  Stream  samples  for  the  field  pH 
analysis  were  collected  in  polyethylene  beakers  from  the  center  of  the  stream. 
The  pH  was  allowed  to  stabilize  before  a  reading  was  taken  in  a  quiecent 
sample.  Lake  samples  for  the  field  analysis  were  collected  in  polyethylene 
bottles  and  brought  immediately  to  shore  where  pH  and  conductivity  measure- 
ments were  taken.  pH  was  also  measured  in  the  laboratory  as  part  of  the 
alkalinity  titration. 

Laboratory  analyses.  Calcium,  magnesium,  sodium,  and  potassium  were 
analyzed  by  atomic  adsorption  spectroscopy  using  an  instrumentation  labora- 
tories model  751  AA/AE  spectrophotometer.  Samples  were  spiked  with  a  solution 
of  lithium  and  lanthanum  to  suppress  ionization  in  the  magnesium  and  calcium 
analysis  (E.M.S.L.  1978).  A  standard  curve  consisting  of  one  blank  and  five 
standards  was  prepared  between  every  34  samples  to  check  for  drift.  If 
significant  (>10%)  drift  had  occurred,  the  samples  were  reanalyzed.  The 
instrument  was  recalibrated  after  every  standard  curve  (Emmel  1977). 

Sulfate,  phosphate,  nitrate,  chloride,  silicate  and  ammonia  were  analyzed 
by  automated  wet  chemistry  using  a  Technicon  II  Auto-Analyzer.  Sulfate  was 
measured  using  a  modification  of  the  Thorin  technique  developed  by  the 
Norwegian  Institute  for  Air  Research.  Phosphate  and  silicate  were  analyzed  by 
standard  molybdenum  techniques  (T.I.S.  1973b;  T.I.S.  1976b).  Nitrate  plus 
nitrite  was  measured  by  the  standard  cadmium  reduction  technique  (T.I.S. 
1972),  while  ammonia  was  analyzed  using  an  Indophenol  Blue  method  (T.I.S. 
1973a).  Chloride  was  analysed  using  the  standard  ferricyanide  method  modified 
for  low  levels  (T.I.S.  1976a).  In  all  of  the  wet  chemical  analyses,  one  blank 
and  six  standards  were  run  between  every  33  samples. 

Alkalinities  were  measured  using  a  potentiometric  method  developed  by 
Gran  (1952).  A  50  ml  sample  was  titrated  with  0.010  N  HC1  from  pH  4.0  to  pH 
3.3.  The  Gran's  function  of  this  titration  curve  was  extrapolated  to  an 
endpoint  to  determine  the  sample  alkalinity.  A  Fisher  Accumet  420  pH  meter 
and  a  Corning  model  476182  pH  electrode  were  used  for  this  titration.  The 
meter  was  calibrated  with  pH  4.0  and  pH  7.0  buffers  at  the  beginning  and  end 
of  each  set  of  titrations.  Both  samples  and  buffers  used  were  at  room 
temperature  and  air  equilibrium.  An  initial  air-equilibrium  pH  was  recorded 
before  each  titration. 
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A  summary  of  the  analytical  methods  is  presented  in  Table  2. 
Quality  control.   The  study  followed  the  following  protocol  on  quality 
control . 

1.  Precision:  To  determine  precision,  5%  of  all  water-quality  samples 
collected  were  treated  as  replicates.  The  results  demonstrated  a 
precision  for  all  measurements  of  ±10%. 

2.  Accuracy:  In  the  field,  the  accuracy  of  the  field  pH  measurement 
was  assured  by  calibration  of  the  electrode  with  standard  buffers 
before  each  measurement.  The  pH  meter  and  electrodes  were  periodi- 
cally calibrated  with  cold,  dilute  acid  standards. 

Table  2.  Analytical  techniques. 


Measurement 
parameter 


Instrumentation 


Technique  Summary 


Mg 

Ca 
Na 
K 

so4 

N03 

NH4 

CI 

SiO^ 

P0, 


Instrumentation  Laboratory  Atomic 
Absorption  Spectrophotometer,  Model  751 


Instrumentation  Laboratory  Atomic 
Absorption  Spectrophotometer,  Model  751 

Instrumentation  Laboratory  Atomic 
Absorption  Spectrophotometer,  Model  751 

Instrumentation  Laboratory  Atomic 
Absorption  Spectrophotometer,  Model  751 

Scientific  Instruments  Model  A200  with 
custom-designed  manifold 


Scientific  Instruments  Model  A200  with 
stock  manifold 

Scientific  Instruments  Model  A200  with 
stock  manifold 

Scientific  Instruments  Model  A200  with 
modified  manifold 

Scientific  Instruments  Model  A200  with 
Auto-Analyzer  with  stock  manifold 

Scientific  Instruments  Model  A200  with 
Auto-Analyzer  with  stock  manifold 


Conductivity   YSI  model  3403  Cell,  Beckman  RC-16C  meter 

pH  Corning  Model  3  Meter  with  Corning  Model  476182 

combination  electrode 

Alkalinity     Radiometer  auto  burette  with  Fisher  Accumet 
420  meter  and  Corning  Model  476182 
combination  electrode 


Lathanum  added,  aspirated  in  oxidizing 
flame  and  read  at  285.2  nm  using 
deuterium  background  correction 

Lathanum  added,  aspirated  in  oxidizing 
flame  and  read  at  422.7  nm 

Lithium  added,  aspirated  in  reducing 
flame  and  read  at  589.0  nm 

Lithium  added,  aspirated  in  reducing 
flame  and  read  at  766.5  nm 

Modification  of  the  thorin  technique 
developed  by  the  Norwegian  Inst,  for 
Air  Research  (NILU) 

Standard  Cadmium  Reduction  technique 


Standard  indophenol  blue  technique 


Standard  ferri cyanide  method  modified 
for  low  levels 

Standard  molybdenum  blue  technique 


Standard  molybdenum  blue  technique 


Standard  procedure 

Standard  two-point  calibration  with 
pH  7.00  and  pH  4.00  buffers 

Air-equilibrated  Grans  titration 
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The  water-chemistry  laboratory  at  the  University  of  Virginia  used  two 
techniques  to  determine  the  accurary  of  analytical  measurement.  Standard 
additions  were  performed  on  representative  samples.  In  addition,  interlabora- 
tory  calibration  with  EPA  (1982)  and  USGS  (1981)  was  conducted  by  the  analysis 
of  blind  replicates  for  all  ions.  The  results  of  these  procedures  showed  that 
the  accuracy  of  the  analytical  measurements  was  ±10%. 
Soils  and  Surficial  Materials 

William  Locke  (unpublished)  has  produced  a  general  map  of  major  soil 
(landfornrvegetation)  groups  in  Rocky  Mountain  National  Park  which  was 
available  for  our  use.  Within  each  basin  six  to  sixteen  sites,  representative 
of  the  major  soil  groups  and  providing  an  altitudinal  cross  section,  were 
selected  for  soil  sampling  (see  Appendix  A  map).  With  the  aid  of  Locke,  soil 
pits  were  prepared  and  a  soil-sampling  scheme  devised. 

Soil  samples  were  collected  by  excavating  soil  pits  by  hand  to  a  maximum 
depth  of  50  cm.  Samples  weighing  300  g  to  700  g  of  each  soil  horizon  exposed 
were  collected  and  described.  These  samples  were  assumed  to  be  representative 
of  the  soils  encountered  in  the  different  landform  types  in  the  study  area. 
Representative  specimens  of  rock  types  were  collected  and  identified. 
Rock- type  names  conformed  to  those  used  in  the  Park  by  Abbott  (1974),  O'Neill 
(1981)  and  Braddock  (unpublished). 

A  complete  sample  description  included  colors,  texture,  structure, 
rooting  depth,  and  estimates  of  coarse  material.  Slope,  size  of  area,  and 
surrounding  vegetation  was  also  recorded  along  with  photographs  of  the  sites. 
Approximately  80-90  samples  were  collected  during  the  1981  field  season. 

Air  dried  soil  samples  passed  through  a  2  mm  sieve  were  used  in  all 
laboratory  analyses.  (The  procedures  used  are  described  in  detail  in 
Appendix  B).  Briefly,  they  were  as  follows:  pH  was  determined  in  a  1:1  soil 
water  suspension.  Exchangeable  bases  were  extracted  with  1  N  NH.Ac  and 
determined  by  Atomic  Absorption  Spectroscopy.  Extractable  acidity  was 
determined  using  a  BaCl^-Triethanolamine  solution  at  pH  8.0,  followed  by  a 
titration.  Cation  exchange  capacity  was  determined  by  the  sum  of  exchangeable 
bases  and  extractable  acidity.  Organic  matter  was  estimated  by  loss  on 
ignition  at  500°C  and  by  wet  oxidation  in  K-CrpOy-hLSO,  solution  (Mebius 
method).  Percent  sand,  silt,  and  clay  was  determined  by  the  pipette  method 
after  dispersion  in  sodium  metaphosphate  solution. 
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Selected  samples  were  tested  for  sulfate  absorption  by  equilibrating  10 
grams  of  soil  with  25  ml  of  dilute  K^SCL  solution  at  a  pH  of  approximately 
4.1.  The  solutions  were  initially  at  1,  5,  and  10  mg  of  soil  per  liter. 
After  shaking  with  soil  and  separation  by  filtration,  solutions  were  analyzed 
for  sulfate  loss. 

Parent  and  surficial  materials.  Determination  of  the  parent  materials  of 
the  soils  was  accomplished  by  geologic  reconnaissance  of  the  area  adjacent  to 
each  soil  sample  and  by  binocular  microscope  examination  of  the  2-  to  4-mm 
fraction  of  each  soil  sample.  Glacial  till  is  the  parent  material  for  most  of 
the  soil  samples  analyzed  and,  in  turn,  is  derived  from  the  geologic  units 
upvalley.  Therefore,  the  composition  of  the  2-  to  4-mm  fraction  of  the  soil 
samples  is  representative  of  both  the  glacial  till  and  its  parent  material, 
the  bedrock  units  in  the  area.  The  areal  extent  of  mapped  geologic  units, 
surficial  and  bedrock,  in  each  basin  was  determined  using  a  planimetric 
digitizer  and  geologic  maps.  Weathering  characteristics  of  the  various 
bedrock  and  surficial  materials  were  studied,  since  they  directly  affect  the 
type  and  amount  of  soil  present. 

Landform  types  (Table  3)  are  helpful  in  determining  the  parent  material 
of  the  soils  developed,  the  thickness  of  overburden  and  the  existing  vegeta- 
tion. They  convey  surface  features  and  materials,  along  with  associated 
vegetation  that  characterize  mappable  areas.  This  information  can  be  obtained 

Table  3.  Representative  landform  types  (surface  features  and/or  materials) 
observed  in  research  area. 


Exposed  bedrock  Moraine  veneer 

Talus  slope  Wet  meadow 

Unglaciated  regolith  Dry  meadow 

Moraine-ground  Wet  tundra 

lateral  Dry  tundra 

end 
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by  analysis  of  aerial  photographs.   Stream  gradients  for  all  rivers  and 
streams  studied  were  determined  with  a  linear  digitizer. 

Seismic  refraction  surveys  were  conducted  at  seven  sites  to  determine  the 
thickness  and  types  of  surficial  materials  overlying  bedrock.  A  Nimbus  ES-125 
single-channel  signal-enhancement  seismograph  was  used.  The  surveys  were  made 
in  coniferous  forests  on  moraine  veneer,  wet  meadows  and  dry  meadows.  The 
seismic  velocities  obtained  permitted  grouping  of  geologic  materials  into  four 
categories:  soil,  alluvial  deposits,  glacial  till,  and  bedrock.  Numerous 
factors  influenced  these  values.  They  include  thickness  of  organic  material, 
water  content,  clay  content,  percent  of  boulders,  lithologic  makeup  of  glacial 
till  and  type  and  degree  of  weathering  and/or  jointing  of  bedrock. 

RESULTS 

Soils  and  Surficial  Materials 

Characteristics.  The  seismic  analysis  revealed  that  generally,  all  soils 
had  velocities  corresponding  to  moist,  loamy  or  silty  soils,  as  reported  by 
Redpath  (1973),  which  indicate  a  low  clay  content.  The  glacial  till 
velocities  correspond  to  values  obtained  by  Redpath  (1973)  for  dry  glacial 
moraine  deposits  in  the  Sierra  Nevada  of  California.  These  deposits  had  low 
clay  content  similar  to  those  encountered  in  the  RMNP  project  area.  Most 
bedrock  velocities  were  low,  indicating  a  fairly  large  degree  of  weathering 
and/or  jointing.  This  suggests  that  surface  water  penetrates  to  considerable 
depth. 

Seismic  profiling  revealed  the  thickness  of  the  soils  to  be  in  the  range 
of  0  to  1.9  m;  alluvial  materials,  3.0  to  3.7  m;  and  the  glacial  till,  0  to 
7  m  (Table  4).  The  high  values  represent  materials  in  the  center  of  the 
stream  valleys  with  decreasing  thickness  up  the  sideslopes.  Although  these 
depths  are  not  large  by  comparison  with  nonalpine  surficial  materials,  they 
could  be  quite  adequate  for  buffering  of  acid  precipitation  if  sufficient  clay 
and  organic  material  existed.  Stream  gradients  are  generally  indicative  of 
the  landform  type  and  are  easily  obtainable  from  topographic  maps.  The 
steepness  of  gradient  generally  is  indirectly  proportional  to  the  thickness  of 
both  glacial  and  alluvial  overburden  in  the  terrain  encountered. 

The  soils  tend  to  be  coarse,  stony,  and  shallow  with  very  little  evidence 
of  formation  of  secondary  minerals  or  eluviation.  (A  detailed  description  of 
the  soil   characteristics  of  each  sampling  site  may  be  obtained  from  the 
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author).  Slight  development  of  a  color  B  horizon  was  detectable  in  a  few  of 
the  lower  elevation  soils  developed  in  glacial  till.  The  major  variation  in 
soils  within  the  study  area  was  caused  by  thickness  of  deposit,  depositional 
mechanism  (local  alluvium,  till  or  colluvial  deposit)  and  amount  of  accumu- 
lated organic  matter.  Differences  in  mineralogy  of  the  parent  material  was 
not  evidenced  in  soil  chemical  or  physical  properties.  Elevation,  especially 
as  it  relates  to  vegetative  cover,  appears  to  be  important  in  determining  the 
organic  matter  content. 

Coarse  materials  (>2  mm  dia. )  comprised  much  of  the  soil  volume  in  most 
locations.  Some  of  the  meadows,  however,  were  essentially  stone-free.  The 
talus  slopes  and  most  other  locations  had  soils  with  10  to  85%  coarse 
fragments  in  the  samples  collected.  The  collection  process  itself  was  biased 
against  large  fragments  therefore  the  percentages  are  conservative. 

In  general  the  percent  base  saturation  and  pH  of  both  the  organic  layers 
and  mineral  soil  are  low.  All  pH  values  are  acid,  in  the  range  4.2  to  5.6, 
with  occasional  values  as  low  as  3.5  and  as  high  as  6.8.  The  basic  cation 
status  of  these  soils  is  generally  quite  low  due  to  the  scarcity  of  basic 
rocks  in  the  underlying  material.  The  %  base  saturation  of  the  soil  exchange 
complex  ranges  from  as  low  as  2%  upwards  to  90%,  but  very  few  samples  exceed 
70%.  The  average  for  each  of  the  basins  ranges  from  24  to  48%. 

The  pH  occurring  most  often  in  the  mineral  soil  horizons  was  5.2  with  a 
strong  central  tendency  around  5.2.  In  15  of  the  30  locations  the  mineral 
soil  pH  in  the  top  50  cm  was  within  the  range  5.0-5.3  (Figure  3).  The  surface 
organic  layers  had  a  pH  range  similar  to  the  mineral  soils,  but  the  distribu- 
tion was  shifted  slightly  toward  more  acid  values. 

The  clay  content  of  these  soils  was  low.  Only  two  sample  locations 
contained  horizons  with  more  than  20%  clay  (%  of  <2  mm  fraction).  It  is 
apparent  that  most  of  the  exchange  capacity  is  provided  by  the  organic  matter. 
Clay  mineralogy  was  not  determined  in  these  samples,  but  the  residual  CEC 
which  appears  to  be  due  to  the  clay  after  subtracting  the  CEC  commonly 
associated  with  the  organic  fraction,  indicates  a  material  high  in  CEC  such  as 
vermiculite.  Shroba  (1977)  reported  mica  alteration  and  formation  of  10-18A 
mixed  layer  clays  in  Alpine  soils  and  more  extensive  alteration  and  the 
presence  of  vermiculite  in  the  spruce-fir  region  in  the  Rocky  Mountains. 

Sampling  of  soils  was  limited  to  a  few  locations  within  each  watershed. 
The  variability  within  watersheds  was  great  and  masks  many  differences  that 
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Figure  3.  Distribution  of  pH  values  in  the  mineral  soil  and  the 
surface  organic  layer  of  all  four  watersheds. 


may  exist  between  drainage  basins  (Tables  5  through  8).  Fall  River  basin 
(Table  7)  contained  the  only  soil  samples  with  pH  in  excess  of  6.0.  These  two 
mineral  soil  samples  were  from  the  lowest  elevations  sampled  in  this  study. 
As  a  result,  the  average  soil  pH  in  that  basin  was  higher  than  that  in  the 
other  three.  Otherwise,  the  average  pH  in  Fall  River  basin  would  have  been 
5.3,  essentially  the  same  as  those  in  Upper  Colorado  and  East  Inlet  basins. 

The  samples  collected  in  Glacier  Gorge  Basin  (Table  8)  are  the  most  acid. 
Only  one  mineral  soil  sample  exceeded  a  pH  of  5.2.  The  mean  of  all  10  samples 
in  that  basin  was  5.0,  relatively  acid  for  young  soils. 

The  pH  of  surface  organic  layers  in  each  basin  was  usually  slightly  lower 
than  the  mineral  soil  underneath  it.  This  is  a  common  occurrence  in  soils 
with  a  surface  accumulation  of  organic  material  under  coniferous  forests. 

A  cross  section  of  soil  samples  were  selected  for  determination  of 
sulfate  absorption  capacity  which  can  contribute  greatly  to  a  system's 
resistance  to  the  impact  of  acidic  inputs.  Sixteen  soil  samples  were  tested 
at  1,  5  and  10  ug  S/ml  solution  for  their  ability  to  adsorb  sulfate.  None  of 
the  samples  tested  removed  detectable  quantities  of  sulfate  from  solution  at 
any  of  the  three  concentrations.  In  many  cases  the  soils  released  low  amounts 
of  sulfate  to  the  solution.  This  was  probably  released  from  the  organic 
matter  as  a  result  of  the  sample  drying  and  rewetting  and  microbial  activity 
during  storage  and  processing.  (Samples  for  sulfate  retention  determinations 
should  be  stored  moist  at  0°  C  until  the  test  is  performed.)  Based  on  these 
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tests,  it  is  obvious  that  these  soils  have  very  little  if  any  sulfate  absorp- 
tion capacity  at  S  concentrations  up  to  10  ug  S/ml  at  a  pH  of  4.1. 

Soil  buffering  capacity.  The  soils  of  the  Rocky  Mountain  National  Park 
at  the  higher  elevation  (above  2800  m)  are  only  slightly  developed.  They  have 
surface  organic  horizons  in  some  cases  and  usually  have  a  darkened  A  horizon 
underlain  by  slightly  weathered  material  derived  from  granites,  diorites  and 
other  rocks  low  in  bases.  The  soils  are  coarse,  low  in  clay,  low  in  basic 
cations  and  relatively  acid  (modal  pH  5.2).  The  organic  matter  provides  most 
of  the  cation  exchange  capacity  and  there  is  no  evidence  of  sulfate  adsorption 
capacity. 

The  physical  properties  of  these  soils — coarse,  high  in  sand,  low  in  clay 
and  steeply  sloping — promote  rapid  movement  of  water  to  the  streams  and  lakes. 
Residence  times  of  water  will  tend  to  be  short  especially  at  the  higher 
altitudes  where  there  is  less  vegetation,  thinner  soils  and  less  accumulated 
organic  matter.  Much  of  the  water  flow  will  be  rapid  and  through  coarse 
channels  such  that  only  a  small  percentage  will  approach  equilibrium  with  the 
soil  materials. 

The  chemical  properties,  low  base  saturation,  low  pH,  and  low  cation 
exchange  capacity  provide  little  opportunity  for  neutralizing  acidic  inputs  or 
providing  significant  alkalinity  to  the  water  system.  If  acidic  deposition  in 
the  region  reached  pH  values  in  the  low  4.0s  there  would  be  a  potential  for 
the  soil  to  buffer  the  pH  upward  slightly.  However,  since  most  water  moving 
through  these  coarse  soils  will  not  have  a  residence  time  sufficient  for 
equilibration  with  the  soil,  there  would  be  little  effect  on  acidity.  The 
current  situation  with  the  soils  more  acid  than  the  streams  indicates  that  the 
soils  are  not  controlling  the  aquatic  system  pH.  Since  the  dominant  soil  pH 
is  only  slightly  above  5  and  the  exchangeable  basic  cation  supply  is  low, 
these  soils  would  not  be  expected  to  provide  strong  buffering  of  the  ecosystem 
against  acidification. 

Even  though  aluminum  was  not  determined  in  this  study  the  present  pH  of 
these  soils  indicates  that  aluminum  concentrations  in  the  soil  solution  are 
not  high.  A  lowering  of  soil  pH  to  4.5-4.6  would  cause  a  significant  increase 
in  aluminum  in  soil  solution.  Some  additional  aluminum  would  reach  the 
aquatic  system,  but  much  more  information  would  be  needed  to  quantify  the 
expected  change  in  aluminum. 
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The  soil  properties  already  mentioned  plus  the  absence  of  sulfate  absorp- 
tion capacity  indicate  that  these  soils  should  be  considered  in  a  "sensitive" 
category.  They  are  among  the  soils  which  have  low  capacity  to  absorb  and 
neutralize  acidity  and  therefore  the  associated  aquatic  systems  would  not  be 
protected  from  pH  change  by  the  soil  ecosystem  (Galloway  et  al.  1983).  The 
soil  itself  is  probably  resistant  to  rapid  acidification.  The  abundance  of 
relatively  unweathered  minerals  provides  the  soil's  major  buffering  against 
lowering  of  the  pH  much  below  its  present  values. 
Chemistry  of  Surface  Waters 

The  lakes  and  streams  in  RMNP  are  characterized  by  low  ionic  strength. 
An  average  concentration  for  the  inorganic  constituents  are  presented  in 
Table  9.  The  raw  data  for  all  of  the  measured  constituents  of  the  more  than 
150  samples  are  presented  in  Appendix  C,  along  with  sampling  dates  for  each 
subbasin.  On  an  equivalent  basis,  the  relative  concentrations  of  the  base 
cations  are  Ca>>Na>Mg>>K  for  most  of  the  watersheds,  although  magnesium 
concentrations  are  greater  than  sodium  in  the  Upper  Colorado  and  Upper  Fall 
River  Basins. 

The  Ypsilon  Creek,  Tyndall  Gorge,  Loch  Vale,  and  Glacier  Gorge  subbasins 
had  similar  stream  chemistry.  These  watersheds  are  characterized  as  having 
lower  alkalinity  and  base  cation  concentrations  than  the  rest  of  the  study 
area.  Roaring  River  and  East  Inlet  had  slightly  higher  alkalinity,  base 
cation,  and  silica  concentrations.  The  Upper  Fall  River  Basin  had  signifi- 
cantly higher  concentrations  of  these  constituents,  while  the  Upper  Colorado 
River  Basin  had  the  highest  alkalinity,  base  cation,  and  sulfate  concentra- 
tions of  any  of  the  watersheds. 

DISCUSSION 

Relationship  of  Surface  Water  Chemistry  to  Atmospheric  Deposition, 
Elevation,  Soils,  Surficial  Materials  and  Geology 

Calculation  of  "excess"  cations.  The  streams  and  lakes  of  RMNP  have  very 
low  cation  and  anion  concentrations.  In  the  watersheds  underlain  by  granitic 
and  metamorphic  rock,  the  average  CR  (Na  +  K  +  Mg  +  Ca),  is  only  ^85  ueq/1 
(Table  9).  In  such  dilute  waters,  the  atmospherically  deposited  base  cations 
may  comprise  a  significant  fraction  of  the  stream  water  cation  concentration. 
Previous  researchers  have  corrected  surface  water  chemistry  data  for  the 
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atmospherically  deposited  base  cations  by  substracting  the  "seasalt"  contribu- 
tion (Henriksen  1980;  Aimer  et  al.  1978).  Although  RMNP  does  not  receive  sea 
salt  in  its  precipitation,  the  evaporite  basins  and  arid  regions  upwind  of  the 
park  result  in  the  wet  and  dry  deposition  of  base  cations  in  the  form  of  dust 
and  salts  (Gosz  1975).  In  an  effort  to  separate  "excess"  cations  (those 
released  as  a  result  of  acid  buffering  reactions)  from  atmospherically 
deposited  cations,  we  have  made  a  correction  for  atmospheric  deposition. 
Chloride,  a  conservative  element  with  no  terrestrial  source  (i.e.,  not  present 
in  the  bedrock)  is  used  to  make  this  correction.  The  calculation  is  made  by 
subtracting  from  each  cation  concentration  the  product  of  the  chloride 
concentration  at  that  sample  site  and  the  ratio  of  that  element  to  chloride  in 
precipitation  (NADP  1982).   The  result  is  the  precipitation  corrected  or 

"excess"  concentration — e.g.,  Na*  =  Na  -  (CI  x  Na    •  )/Cl    •    (* 

3  '  precip.    precip. 

indicates  "excess"  concentration).  Some  base  cations,  especially  potassium, 
are  taken  up  by  the  ecosystem.  As  a  result,  our  precipitation-corrected 
concentration  of  this  element  may  turn  out  to  be  negative.  If  the  calculation 
results  in  a  negative  value  for  an  element,  the  precipitation-corrected 
concentration  is  taken  as  zero.  The  results  of  these  calculations  are 
presented  in  Table  10.  This  method  of  calculating  the  "excess"  base-cation 
concentration  (Na*  +  K*  +  Mg*  +  Ca*  =  Cn*)  is  a  correction  for  the  maximum 
contribution  of  C,  (Na  +  K  +  Mg  +  Ca  =  CR)  by  wet  deposition. 

The  stochiometry  of  primary  dissolution  reactions  provides  a  check  on 
this  correction,  since  the  release  of  base  cations  by  mineral  dissolution  is 
accompanied  by  the  concurrent  release  of  hLSiO.,  e.g.: 

NaQ  7CaQ  3A11  3Si'2  70g  +  3.45  H20  +  1.3  H2C03  =  0.65  Al 2Si 205(0H)4 
(oligoclase)  (kaolinite) 

+1.4  H4Si04  +0.7  Na+1  +0.3  Ca+2  +1.3  HC0"1 


3  K2(Mg3Fe3)Al2Si6020(0H)4  +  24  H20  +  18  H2C03  = 
(biotite) 
2(Mg3Fe3)Al3Si5020(0H)4-8  H20  +  6  K+1  +  3  Mg+2  +  3  Fe+2  +  8  H4Si04  +  18  HC0"1. 
(vermicullite) 
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Table  10.  Comparison  of  mean  of  uncorrected  (U)  and 
"excess"  (E)  base  cations.3 


Watershed 

CB 

Na 

K 

Mg 

Ca 

Roaring  River  (U) 
(E) 

121.0 
72.3 

29.0 
21.9 

3.6 
0.0 

21.7 
11.7 

66.5 
39.3 

Ypsilon  Creek  (U) 
(E) 

86.2 

56.2 

22.9 
18.8 

4.1 
0.3 

15.2 
9.4 

43.6 
28.2 

Upper  Fall  River  (U) 
(E) 

222.0 
155.2 

40.2 
30.7 

7.5 

0.0 

67.3 
53.8 

106.9 
71.2 

Andrews  Creek  (U) 
(E) 

88.9 
57.3 

16.1 
11.9 

3.7 
0.0 

13.1 
7.3 

55.7 
38.1 

Loch  Vale  (U) 
(E) 

86.2 
47.9 

15.0 
9.6 

3.1 
0.3 

13.9 
6.1 

52.9 
32.4 

Glacier  Creek  (U) 
(E) 

73.8 
47.6 

14.2 
10.5 

3.0 
0.2 

10.0 
4.8 

46.3 
32.5 

Upper  Colorado  River 

(U) 
(E) 

357.1 
280.1 

45.6 
35.1 

7.7 
1.4 

77.8 
63.0 

218.1 
181.0 

East  Inlet  (U) 
(E) 

125.5 
91.5 

26.5 
22.5 

2.5 
0.0 

16.0 
8.8 

90.0 
60.6 

All  concentrations  in  ueq/1. 

We  can  assess  the  accuracy  of  our  estimate  of  CR*  by  comparing  the  corrected 

-4  *     .  -4 

Cn  to  silica  concentrations  (as  moles  of  SiO,  ).  Graphs  of  Cn*  vs.  SiO,  for 

Loch  Vale  and  Ypsilon  Creek  watersheds,  with  low  intercepts  show  that  our 
correction  is  appropriate  for  the  watersheds  with  lowest  alkalinity  (Figure 
4).  The  variation  in  the  slopes  of  the  regression  lines  among  the  watersheds 
reflects  a  difference  in  the  distribution  of  minerals.  Watersheds  with  a 
slightly  greater  percentage  of  biotite  will  have  larger  slopes  than  those  with 
more  oligoclase  in  their  drainage  according  to  the  different  stochiomatry  of 
the  weathering  reactions.  This  check  on  our  calculation  of  Cn*  is  not 
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Figure  4.  "Excess"  base  cations  vs.  silica  for  two  sensitive 
watersheds  in  RMNP. 


appropriate  in  watersheds  where  nonsilicious  minerals  and  soil  cation-exchange 
reactions  play  a  large  part  in  the  geochemistry  (i.e.,  the  Upper  Colorado 
River  and  parts  of  the  Fall  River  basins). 

Atmospheric  deposition  and  water  chemistry.  Researchers  have  shown  that 
the  predominent  direction  of  air  mass  movement  over  the  Front  Range  of  the 
Central  Rocky  Mountains  is  from  west  to  east  (Barry  1973),  with  periodic 
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upslope  movement  from  the  east  (Kelley  and  Stedman  1980).  Although  most  of 
the  precipitation  received  by  this  area  is  dropped  by  the  westerly  air  masses, 
the  upslope  transport  of  pollutants  from  sources  to  the  east  may  have  a 
profound  influence  on  the  deposition  chemistry  of  the  Front  Range  (Kelley  and 
Stedman  1980).  The  atmospheric  deposition  from  the  westerly  air  masses 
contains  little  sea  salt  (Junge  and  Werby  1958)  but  does  result  in  the  wet  and 
dry  deposition  of  airborne  dust  and  salts  from  the  dry,  windy  areas  to  the 
west  (Gosz  1975).  As  a  result,  RMNP  precipitation  has  a  relatively  high  base 
cation  concentration  [22.0  ueq/1  Ca  in  RMNP,  compared  with  13.0  peq/1  in  the 
Adirondacks  (NADP  1982,  Altwicker  and  Johannes  1980)]. 

The  nitrate  concentrations  in  RMNP  are  relatively  high,  comprising  35.6% 
of  the  acid  anions  in  precipitation,  compared  with  28.7%  at  Hubbard  Brook 
(NADP  1982).  These  high  nitrate  concentrations  are  probably  the  result  of 
upslope  transport  of  NO  from  the  Denver  metropolitan  area  (Kelley  and  Stedman 
1980).  The  bedrock  in  RMNP  contains  only  trace  amounts  of  chloride  (Lovering 
and  Goddard  1950),  while  sulfur-bearing  minerals  (e.g.,  chalcophyrite  and 
gypsum)  are  found  only  in  the  Upper  Fall  River  and  Colorado  River  basins.  As 
a  result,  atmospheric  deposition  is  expected  to  be  the  primary  source  of 
chloride  and  nitrate  in  all  of  the  streams  in  RMNP,  as  well  as  the  primary 
source  of  sulfate  in  the  waters  of  East  Inlet,  Glacier  Gorge,  Ypsilon  Creek, 
and  Roaring  River. 

Stream  water  chloride  and  nitrate  concentrations  are  fairly  uniform  over 
the  park,  although  East  Inlet,  on  the  west  side  of  the  continental  divide,  has 
lower  nitrate  concentrations  than  the  other  four  watersheds  (Table  9).  Two 
plausible  explanations  for  the  lower  nitrate  in  East  Inlet  are  (1)  less 
deposition  from  upslope  air  masses  bearing  NO  from  the  Denver  area  and/or  (2) 
more  biological  activity.  Chloride  and  sulfate  concentrations  show  little 
variability  with  elevation,  indicating  negligible  variation  in  the  atmospheric 
contribution  of  elements  with  elevation  (Figures  5  and  6).  Nitrate,  the  most 
biologically  active  anion,  shows  some  variation  over  elevation,  with  the 
highest  nitrate  concentrations  above  the  timberline,  where  biological  activity 
is  lowest  (Figure  7). 

Streamwater  chemistry  and  elevation.  Elevation  is  related  to  several 
parameters  that  may  effect  streamwater  chemistry.  The  thickness  of  surficial 
deposits  (glacial  till,  alluvium,  loess,  and  soils)  and  the  length  of  flow 
path  increase  with  decreasing  elevation  as  we  descend  from  headwater  lakes 
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Figure  5.  Graphs  of  CI  vs.  elevation  for  subbasins  in  RMNP. 
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Figure  6.  Graphs  of  SO.  vs.  elevations  for  subbasins  in 
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Figure  7.   Graphs  of  N03  vs.  elevation  for  subbasins  in  RMNP. 
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situated  on  scoured  bedrock  to  lower  elevation  lakes  surrounded  by  moraine 
veneer  with  thin  soil  cover.  In  addition,  watershed  slope  decreases  with 
decreasing  elevation,  further  increasing  the  contact  time  between  precipita- 
tion and  neutralizing  materials.  These  parameters  are  all  difficult  to 
measure.  Turk  and  Adams  (1983)  have  shown  that  elevation  used  as  a  surrogate 
for  these  and  other  unmeasured  watershed  parameters  (e.g.,  residence  time, 
length  of  flow  path,  effective  soil  and  bedrock  area),  successfully  predicts 
alkalinity.  Our  analysis  shows  that  alkalinity,  CR*  (excess  CR),  and  silicate 
are  all  inversely  related  to  elevation  in  subbasins  with  homogeneous 
mineralogy  and  low  alkalinities  (Glacier  Gorge,  Loch  Vale,  Ypsilon  Creek, 
Roaring  River,  and  East  Inlet)  (Figures  8  through  10). 

This  relationship  between  elevation  and  streamwater  chemistry  is  obscured 
when  the  data  from  more  than  one  subbasin  are  used  in  the  analysis.  This 
probably  results  from  the  differences  in  the  distribution  of  glacial  till  over 
elevation  between  subbasins  in  the  same  watershed.  The  relationship  between 
elevation  and  streamwater  chemistry  is  overwhelmed  by  the  variability  in 
mineralogy  and  subsequent  weathering  rates  in  the  Upper  Fall  River  and  Upper 
Colorado  River  basins  (Figure  8).  Although  elevation  contributes  to  the 
variability  of  stream  chemical  composition  in  areas  with  homogeneous  bedrock, 
variability  in  mineralogy  and  other  interwatershed  variables  can  overwhelm 
altitudinal  effects. 

Relationships  between  bedrock  geology,  surficial  materials,  soils 
and  stream  water  chemistry.  The  chemical  composition  of  natural  waters  is 
primarily  a  function  of  the  interactions  between  atmospheric  deposition, 
bedrock  geology,  and  surficial  deposits.  Previous  researchers  have  found  that 
an  area's  ability  to  neutralize  acidic  deposition  is  the  result  of  mineral 
weathering  and  soil  ion-exchange  reactions  (Norton  1980;  McFee  et  al.  1977; 
Johnson  and  Cole  1980).  Under  pristine  conditions,  the  chemical  dissolution 
of  minerals  by  carbonic  acid  results  in  a  loss  of  hydrogen  ion  (acidity)  and 
the  production  of  alkalinity  (acid-neutralizing  capacity),  base  cations,  and 
silica.  The  rates  of  these  reactions  vary  with  the  mineralogy  of  the  bedrock. 
Carbonaceous  minerals,  such  as  calcium  carbonate,  dissolve  very  rapidly, 
providing  "infinite"  buffering  capacity,  while  silicious  rocks,  such  as 
granite,  dissolve  very  slowly,  offering  little  buffering  capacity.  Surficial 
materials  also  play  an  important  role  in  buffering  acidic  deposition, 
as   soils  neutralize  acidic  deposition  through  cation  and  anion  exchange. 
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Figure  8.  Graphs  of  alkalinity  vs.  elevation  for  subbasins  in  RMNP, 


37 


200  — i 


ROARING  RIVER 


o  o 


o  o 


o 


0 


I  I  I  I  I  I  I  I  I  I  I  I  I  I 


2500  3000  3500 

ELEVATION   (m) 


4000 


200  — l 


-5    100  - 

*     CO 


0 


YPSILON  CREEK 


o 


o   o 


o 


o 


o    do 


1 — I — I — I — I — I — 

2500  3060 

ELEVATION   (m) 


8 


3580 


500  -i 


UPPER   FALL   RIVER 

o 


o 

K 

o  o 

I   I  I  I  I  I  I  I  I  I  I  I  I 


00,0 


0 


o 


2500  3000  3500  4000 

ELEVATION   (m) 


200—1 


5    iW-1 


0 


LOCH  VALE 


•o 


0    °o     o0    % 


I  I  I  I  I   I  I   I  I  I  I  I  I  I 

2750     3000     3250     3500 

ELEVATION  fm) 


100  — 1 


GLACIER  CREEK 


o 


o 


o 


Po 


o 
o 


o 


o 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

2500     3000     3500     4000 

ELEVATION  (m) 


500  — i 


0  — 

0            UPPER 

COLORADO  RIVER 

oX 

*    °0 

00        u 

0 
0 

0 

1    1    ■ 

T"-1-! 

200      400 

ELEVATION  (m) 


600 


150  — 

EAST 

INLET 

: 

0 

0 

^    100  ~ 

cr 
<u 

CO 

°0 
0 

°cP° 

50 

1   1   ' 

1     '     1 

1 

"1 

2600   2800   3000   3200   3400 

ELEVATION  (m) 


Figure  9.  Graphs  of  excess  base  cations  vs.  elevation  for 
subbasins  in  RMNP. 
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Figure  10.   Graphs  of  silicate  vs.  elevation  for  subbasins  in  RMNP. 
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Surface  water  chemistry,  particularly  alkalinity  is  an  integration  of  these 
mechanisms.  As  a  result,  alkalinity  is  frequently  used  as  an  index  of  sensi- 
tivity. 

The  geology  in  the  Park  ranges  from  Precambrian  granite  to  tertiary 
intrusive  and  extrusive  rocks.  The  difference  in  the  weathering  rates  of  the 
minerals  in  these  formations  is  reflected  by  the  streamwater  chemistry.  The 
highest  concentrations  of  CR*,  alkalinity,  and  silicate  occur  in  the  Upper 
Colorado  River  Basin,  an  area  underlain  by  highly  weatherable  ash  flow  tuff 
and  andesite.  The  CR*  and  alkalinity  concentrations  in  Glacier  Creek,  a 
watershed  underlain  by  Silver  Plume  granite,  are  an  order  of  magnitude  lower 
than  those  in  the  Upper  Colorado  River  Basin  (Tables  9  through  10). 

R-mode  factor  analysis  was  used  to  investigate  the  relationships  among 
the  water  chemistry  variables  and  the  bedrock  and  soil  parameters.  This 
method  of  analysis  has  been  used  by  several  authors  (Dawdy  and  Feth  1967; 
Miller  and  Drever  1977a;  Reeder  et  al.  1972)  to  interpret  natural  water  chem- 
istry in  terms  of  its  geochemical  origin.  In  these  analyses,  factor  analysis 
is  used  to  reveal  the  geochemical  processes  responsible  for  the  stream  chemis- 
try composition.  A  brief  review  of  factor  analysis  follows. 

Factor  analysis  is  a  statistical  data  reduction  technique  that  rearranges 
the  information  contained  in  the  correlation  matrix  for  a  set  of  variables 
into  a  smaller  set  of  independent  factors.  These  factors  are  linear  combina- 
tions of  the  original  variables.  The  first  step  in  the  analysis  calculates 
the  principal  components,  or  factors  that  explain  the  greatest  amount  of  the 
variances  and  co-variances  in  the  correlation  matrix.  In  the  varimax  solution 
used  in  this  analysis,  these  principal  components  are  rotated  orthogonally  to 
achieve  a  simple  structure.  This  rotation  produces  a  set  of  uncorrelated 
factors,  so  that  the  factor  loadings  (the  extent  to  which  each  factor  is 
associated  with  a  particular  variable)  tend  towards  unity  or  zero.  The  re- 
sulting factors  are  interpreted  as  source  variables  accounting  for  the  rela- 
tionships between  the  original  variables. 

The  data  set  from  RMNP  was  split  into  two  geologic  groups  for  this 
analysis.  One  group  represents  all  samples  from  areas  underlain  by  granite 
and  gneiss;  the  other  represents  the  samples  collected  from  watersheds  con- 
taining tertiary  volcanics.  Each  group  was  analyzed  with  and  without  a  bed- 
rock geology  and  soil  parameter.   The  analysis  of  the  stream  chemistry  alone 
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allowed  a  larger  sample  size.  Variables  that  had  many  zero  values  (K*  and 
NhL)  were  not  included  in  the  analysis.  All  factors  with  eigen  values 
greater  than  0.51  are  reported. 

The  results  of  the  factor  analysis  of  the  granite  and  gneiss  group  show 
that  four  factors  account  for  nearly  93%  of  the  variance  in  stream  water 
chemistry  data  (Table  11).  The  first  factor,  accounting  for  almost  50%  of  the 

Table  11.  R-mode  varimax  factor  matrix  of  chemical  data 
for  88  lake  and  stream  samples  underlain  by  granite  and 
biotite  gneiss  and  schist  in  Rocky  Mountain  National  Park. 


Factor 

Variable 

1 

2 

3 

4 

Na+' 

0.934 

0.243 

-0.185 

-0.042 

Mg 

Ca+2 

-1 

0.827 

0.346 

0.092 

-0.127 

0.325 

0.919 

0.144 

0.068 

Cl" 

•1 

-0.031 

0.082 

0.112 

0.975 

N03 

0  2 

-0.150 

-0.025 

0.945 

0.078 

SO4 

0.259 

0.437 

0.609 

0.479 

Si'04 

0.937 

0.156 

-0.029 

0.180 

Alk 

0.538 

0.656 

-0.409 

0.277 

%  of  variance 

1 

explained 

by 

factor 

49.5 

23.6 

12.0 

7.7 

Cumulative 

:  % 

of 

variance 

49.5 

73.1 

85.1 

92.8 

variance,  has  high  loadings  for  Na*,  Mg*,  SiO,,  and  alkalinity.   This  factor 
is  interpreted  as  an  oligoclase  and  biotite  weathering  factor.   These  two 

minerals  weather  to  kaolinite  and  vermicullite,  respectively,  releasing  Na  , 

+2         -1  +2 

Mg  ,  Si(L,  HCOt,  and  a  small  amount  of  Ca  . 
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NaQ  7CaQ  3A11  3Si'2  ?0g  +  3.45  H20  +  1.3  H2C03  =  0.65  Al 2Si 205(0H )4 
(oligoclase)  (kaolinite) 

+1.4  H4Si04  +0.7  Na+1  +0.3  Ca+2  +1.3  HCO"1 


3  K2(Mg3Fe3)A12Si6020(0H)4  +  24  H20  +  18  H2C03  = 
(biotite) 
2(Mg3Fe3)Al3Si5020(0H)4-8  H20  +  6  K+1  +  3  Mg+2  +  3  Fe+2  +  8  H4Si04  +  18  HCO"1. 
(vermicul lite) 

The  high  percentage  of  the  variance  explained  by  this  oligoclase  and 
biotite  weathering  factor  suggests  that  primary  mineral  weathering  is  one  of 
the  processes  dominating  stream  water  chemistry. 

The  second  factor,  accounting  for  23.6%  of  the  variance,  is  less  easily 
interpreted.  This  factor  loads  heavily  on  Ca  ,  HC03,  and  S04.  Ca  '  is  not  a 

major  bedrock  element,  although  a  small  percentage  is  contained  in  the 

+2 
oligoclase.   Ca   in  precipitation,  however,  is  quite  high  and  is  also  con- 
tributed to  by  dry  deposition  of  dust  and  salts.   This  factor  can  be  inter- 
preted as  a  dry  deposition  factor,  caused  by  the  deposition  of  CaS04,  and 
possibly  CaC03> 

The  third  factor,  accounting  for  12%  of  the  variance,  is  an  acid  deposi- 

-1-2  • 

tion  factor.   The  high  loadings  of  N03  ,  S04  ,  with  an  inverse  relationship 

between  the  acid  anions  and  alkalinity,  represents  the  titration  of  lake  and 
stream  HC03  by  anthropogenic  acids.  The  fact  that  N03  loads  higher  on  this 
factor  suggests  that  HN03  may  be  a  more  important  component  of  acid  precipita- 
tion than  hLS04.  The  fourth  factor,  explaining  7.7%  of  the  variance  repre- 
sents atmospheric  deposition  but  does  not  appear  to  involve  any  buffering 
mechanisms. 

The  addition  of  percent  unreactive  rock  (quartz  +  microcline)  and  percent 
soil  organic  matter  (%  O.M.)  to  the  factor  analysis  results  in  a  decrease  in 

sample  size  and  the  loss  of  the  acid  deposition  factor  (Table  12).  The  first 

+2 
factor  still  represents  oligoclase  and  biotite  weathering,  although  Ca   loads 
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Table  12.  R-mode  varimax  factor  matrix  of  stream  chemistry,  %  soil 
organic  matter,  and  %  unreactive  rock  for  40  lake  and  stream  samples 
underlain  by  granite  and  biotite  gneiss  and  schist  in  Rocky  Mountain 
National  Park. 


Variable 

Factor 

1 

2 

3 

4 

5 

+  1 
Na 

0.958 

-0.064 

0.111 

-0.066 

0.164 

Mg 
Ca 

i 

0.917 

0.208 

0.084 

-0.090 

-0.154 

0.558 

0.727 

0.354 

-0.042 

-0.048 

Cl~ 

1 

0.009 

0.273 

0.212 

0.923 

-0.082 

NOl 

12 
S04 

A 

-0.226 

0.851 

-0.212 

0.286 

-0.093 

0.275 

0.701 

-0.172 

0.543 

-0.116 

Si04 

0.917 

0.053 

0.024 

0.238 

0.168 

Alk 

0.741 

-0.066 

0.591 

0.140 

0.049 

%  O.M. 

0.101 

-0.107 

0.064 

-0.087 

0.978 

%  rock 

-0.118 

0.085 

-0.953 

-0.112 

-0.062 

%  of  variance 
explained  by 

1 

factor 

41.7 

25.1 

12.4 

9.0 

6.1 

Cumulative 
variance 

!  % 

of 

41.7 

66.8 

79.2 

88.2 

94.3 

somewhat  higher  than  in  the  previous  analysis.  The  second  factor  still  repre- 
sents atmospheric  deposition  of  CaSO*  in  dust  and  salts,  although  NG\  also 
loads  heavily  on  this  axis.  The  third  factor  represents  the  inverse  relation- 
ship between  HCOo  and  unreactive  rock,  further  evidence  that  primary  mineral 
weathering  plays  a  dominant  role  in  stream  chemistry.  The  fourth  factor 
represents  atmospheric  deposition.  The  fifth  factor  represents  %  O.M. ,  which 
is  unrelated  to  any  other  variable,  suggesting  that  soil  ion-exchange 
processes  play  a  very  minor  role  in  stream  chemistry.  The  factor  analysis  of 
stream  samples  from  areas  underlain  by  tertiary  volcanics  shows  that  four 
factors  account  for  86.4%  of  the  variance  (Table  13). 

The  first  factor,  accounting  for  over  40%  of  the  variance,  indicates  that 
the  weathering  processes  in  these  watersheds  are  very  different  from  those  in 
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Table  13.   R-mode  varimax  factor  matrix  of  stream  chemistry 
for  33  lake  and  stream  samples  underlain  by  tertiary 
intrusive  bedrock  in  Rocky  Mountain  National  Park. 


Factor 

Variable 

1 

2 

3 

4 

Na+' 

+  2 

Mg 
Ca+2 

0.122 

0.911 

-0.225 

-0.076 

0.835 

0.152 

-0.264 

0.129 

0.772 

-0.410 

0.138 

0.115 

Cl" 

1 

0.277 

0.033 

-0.025 

0.955 

NO, 

12 

so4 

-0.036 

-0.219 

0.946 

-0.025 

0.853 

0.121 

0.185 

0.095 

Si04 

0.006 

0.884 

-0.066 

0.121 

Alk 

0.881 

0.127 

-0.184 

0.247 

%  of  variance 

\ 

explained 

by 

factor 

40.3 

26.5 

10.0 

9.6 

Cumulative 

!  % 

of 

variance 

40.3 

66.8 

76.8 

86.4 

-1    -2    +2       +2 
the  granite  and  gneiss  areas.  HC(L  ,  SO.  ,  Ca  ,  and  Mg   all  load  heavily  on 

this  factor,  representing  the  weathering  of  mafic  minerals  low  in  silica,  as 

well  as  the  dissolution  of  CaSO,  in  the  Pierre  shale  of  the  Upper  Colorado 

Basin.   The  second  factor  explains  26.5%  of  the  variance  and  loads  heavily  on 

+1       -4 
Na  '  and  Si0»  .   This  factor  can  be  described  as  oligoclase  weathering.  The 

third  and  fourth  factors,  accounting  for  19.6%  of  the  variance,  each  explain  a 

single  variable  and  cannot  be  interpreted  in  terms  of  buffering  mechanisms. 

The  addition  of  soil  and  bedrock  parameters  in  this  analysis  reduced  the 

sample  size  to  16  observations,  too  few  to  be  used  with  factor  analysis. 

The  results  of  the  analyses  of  the  two  geologic  groups  suggest  that 

primary  mineral  weathering  is  the  major  mechanism  underlying  the  stream 

chemistry.   The  weathering  of  oligoclase  and  biotite,  the  dry  deposition  of 

CaSO.  dust,  and  the  titration  of  bicarbonate  by  anthropogenic  acids  appear  to 

be  the  major  geochemical  processes  in  the  watersheds  underlain  by  granite  and 

biotite  gneiss  and  schist.  An  assemblage  of  low-silica,  high-sulfate  minerals 
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appears  to  be  responsible  for  much  of  the  stream  chemistry  in  areas  underlain 
by  tertiary  intrusives.  Oligoclase  weathering  plays  an  important  role  in  the 
chemistry  of  these  watersheds,  while  no  evidence  of  acid  titration  of  bicarbo- 
nate was  found. 

Discriminant  function  analysis  (DFA)  was  used  to  examine  the  differences 
in  stream  chemistry  between  the  areas  underlain  by  tertiary  volcanic  bedrock 
and  those  surrounded  by  granite  and  biotite  gneiss  and  schist.  DFA  dis- 
tinguishes between  two  populations  on  the  basis  of  observations  of  multiple 
variables.  DFA  may  also  be  used  to  classify  data  on  the  basis  of  observed 
variables.  In  this  two-group  discriminant  analysis,  one  discriminant  function 
is  formulated  by  the  analysis.  This  function  is  a  linear  combination  of  the 
measured  variables,  i.e., 

L  =  B1X1  ♦  B2X2  +  .  .  .  +  BnXn  , 

where  L  is  the  discriminant  function,  X  are  the  measured  variables,  and  B 

'  n  n 

are  the  discriminant  function  coefficients.  The  discriminant  function  L  is 
formulated  to  achieve  the  maximum  discrimination  between  the  two  groups;  i.e., 
the  variance  in  L  within  each  group  is  much  less  than  the  variance  in  L 
between  the  two  groups.  The  discriminant  function  analysis  also  indicates 
which  variables  differ  most  between  the  two  groups.  The  correlation  between 
the  measured  variables  and  the  discriminant  function  (ranges  between  -1.0  and 
+1.0)  gives  an  indication  of  the  importance  of  that  variable  in  differen- 
tiating between  the  two  groups. 

In  this  application,  DFA  is  used  to  differentiate  between  the  two  major 
geologic  groups  (tertiary  volcanics  vs.  granite  and  biotite  gneiss  and  schist) 
on  the  basis  of  the  water  chemistry  parameters.  The  results  of  the  DFA  show 
that  the  two  groups  are  significantly  different  (Table  14).  The  correlations 

between  the  water  chemistry  parameters  and  the  discriminant  function  show  that 

+2    +2        -2 
alkalinity,  Mg  ,  Ca  ,  and  SO,   are  the  parameters  responsible  for  the 

difference  between  the  two  geologic  groups.  These  are  the  same  variables  that 

make  up  the  first  factor  in  the  factor  analysis  of  the  tertiary  volcanic 

group.   This  result  is  consistent  with  the  belief  that  the  differences  in 

stream  chemistry  between  the  two  geologic  groups  are  caused  by  mineral 

weathering  processes.   The  reclassification  of  the  data  set  on  the  basis  of 
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Table  14.  Discriminant  analysis  of  surface 
water  chemistry  from  117  lake  and  stream 
samples  from  Rocky  Mountain  National  Park. 


Variable 


Alkalinity 
Mg+2 
Ca+2 
SO"2 

cr1 

Na+1 
NO"1 


Correlation 

0. 

854 

0. 

823 

0. 

615 

0. 

558 

0. 

348 

0. 

337 

0. 

297 

-0. 

173 

the  discriminant  function  resulted  in  a  95.7%  correct  classification 
(Table  15). 

The  results  of  this  analysis  show  that  water  chemistry  characteristics  of 

areas  underlain  by  similar  bedrock  geology  can  be  estimated  on  the  basis  of 


Table  15.   Reclassification  results. 


Predicted  group  membership 
Actual  group  No.  of  cases      


1 


Group  1  (granite  and  gneiss)        88  88  0 

100.0  0 

Group  2  (tertiary  intrusives)        29  5  24 

17.2  82.8 

Percent  of  grouped  cases  correctly  classified  =  95.73 
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geologic  type.  Areas  with  tertiary  intrusive  rocks  present  in  their  catch- 
ments can  be  expected  to  have  significantly  higher  alkalinity,  magnesium, 
calcium,  and  sulfate  than  areas  underlain  by  granite  and  gneiss. 

Summary.  The  chemical  composition  of  the  waters  in  RMNP  is  a  product  of 
the  interactions  between  atmospheric  deposition  and  bedrock  geology  and 
surficial  materials.  Soils  play  a  minor  role.  The  previous  discussions  show 
that  bedrock  mineralogy,  atmospheric  deposition,  and  elevation  all  signifi- 
cantly affect  the  stream  chemistry  in  RMNP. 

Primary  mineral  weathering  appears  to  be  the  dominant  mechanism  deter- 
mining the  concentrations  of  base  cations,  silica,  and  alkalinity  throughout 
the  park.  Factor  analysis  shows  that  the  primary  mineral  weathering  of 
oligoclase  and  biotite  account  for  almost  50%  of  the  variance  in  stream 
chemistry  in  areas  underlain  by  granite  and  biotite  gneiss  and  schist.  In 
areas  that  contain  tertiary  volcanic  bedrock,  the  mineral  weathering  of  mafic 
minerals,  sulfur-bearing  minerals,  and  oligoclase  account  for  more  than  50%  of 
the  variance  in  stream  chemistry.  Soils  appear  to  have  little  effect  on 
stream  chemistry.  The  soils  are  highly  permeable,  low  in  clays,  and  very 
thin.  Soil  organic  matter  accounts  for  most  of  the  soil  CEC,  but  %  O.M.  does 
not  form  a  factor  with  any  of  the  stream  chemistry  variables. 

Alkalinity,  C$,  and  sulfate  are  inversely  related  to  elevation  in  the 
Glacier  Creek,  Loch  Vale,  Ypsilon  Creek,  Roaring  River,  and  East  Inlet  sub- 
basins.  This  is  a  result  of  the  deeper  glacial  till  (larger  flow  path)  and 
gentle  gradients  (larger  residence  time)  at  lower  elevations. 

Atmospheric  deposition  is  the  primary  source  of  chloride  and  nitrate  in 
the  streams  of  RMNP.  Atmospheric  deposition  is  also  the  primary  source  of 
sulfate  in  RMNP  waters,  with  the  exception  of  the  Upper  Colorado  and  Upper 
Fall  River  basins,  where  the  weathering  of  sulfur-bearing  minerals  is  a  source 
of  sulfate.  Chloride  and  sulfate  are  relatively  constant  with  elevation, 
while  nitrate  concentrations  are  highest  above  the  timberline,  where  biolog- 
ical activity  is  lowest.  Atmospheric  deposition  is  also  a  significant  source 
of  stream  water  CR,  a  result  of  the  deposition  of  airborne  dust  and  salts  from 
the  dry,  windy  regions  upwind  of  RMNP. 

Discriminant  function  analysis  of  the  water  chemistry  data  shows  that 
bedrock  mineralogy  can  be  used  to  estimate  the  water  chemistry  of  RMNP. 
Bedrock  mineralogy  will  be  used  in  later  sections  to  assess  the  sensitivity 
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(as  measured  by  alkalinity)  of  similar  geologic  areas  in  the  central  Rocky 

Mountain  region. 

Alkalinity  (Sensitivity)  of  RMNP  Lakes  and  Streams 

An  area's  ability  to  neutralize  acidic  deposition  through  interactions 
with  bedrock  and  surficial  geology  determines  its  vulnerability  to  acidic 
deposition.  The  lithological  characteristics  of  a  watershed,  combined  with 
its  water  chemistry  data,  provide  the  basis  for  assessing  the  area's  sensi- 
tivity to  acid  rain.  Alkalinity,  an  integrator  of  watershed  buffering  mecha- 
nisms, is  often  used  as  an  index  of  sensitivity.  Hendrey  et  al.  (1980)  define 
sensitive  waters  as  those  with  alkalinities  lower  than  200  ueq/1 ,  a  level  low 
enough  to  be  neutralized  by  acidic  deposition  and  runoff.  Using  this  defini- 
tion, we  find  that  much  of  the  RMNP  is  sensitive  to  acidic  deposition. 

As  we  have  seen  in  the  previous  sections,  the  sensitivity  of  the  water- 
sheds in  RMNP  is  primarily  determined  by  the  bedrock  geology  of  the  water- 
sheds. Within  each  watershed,  elevation  may  be  used  to  further  divide  the 
basins  into  areas  of  different  sensitivities.  The  classification  scheme  used 
for  assessing  the  sensitivity  of  RMNP  to  acidic  deposition  defines  sensitive 
waters  as  those  with  alkalinities  less  than  200  ueq/1,  while  waters  with 
alkalinities  less  than  100  ueq/1  will  be  considered  very  sensitive.  Lakes  and 
streams  with  alkalinities  below  50  ueq/1  are  classified  as  extremely  sensi- 
tive. Following  this  classification  scheme,  we  have  evaluated  the  sensitivity 
above  and  below  3300  m  in  each  subbasin  (Figure  11).  The  results  show  the 
following  classifications: 

Extremely  sensitive  Glacier  Gorge 

(alkalinity  S  50  ueq/1)  Ypsilon  Creek 

Very  Sensitive  Roaring  River 

(50  %   alkalinity  ^  100  ueq/1)         Upper  East  Inlet 

Sensitive  Lower  East  Inlet 

(100  ueq/1  alkalinity  £  200  ueq/1)     High  elevations  of 

Upper  Fal 1  River 

Nonsensitive  Lower  elevations  of 

(alkalinity  £  200  ueq/1)  Upper  Fall  River 


Upper  Colorado  River 


Current  Acidification  Status  of  RMNP 


Researchers  have  observed  acidic  deposition  in  the  Rocky  Mountains  since 
1967  (Lewis  and  Grant  1979).   These  authors  report  a  decrease  in  stream 
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(50  <  alkalinity  <  100)      iSi]  (100  <  alkalinity  <  200 


□ 


Non-sensitive 
(alkalinity  >  200  yeq/i) 


Figure  11.      Spatial   distribution  of  sensitivity  in  RMNP 
(see  map  Figure  1). 
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bicarbonate  in  the  Como  Creek  watershed  of  north-central  Colorado  over  a 
150-week  period.  This  record,  however,  is  too  short  to  be  considered  evidence 
of  increasing  acidification.  Lewis  (1982)  reports  an  average  decrease  in 
alkalinity  from  22  mg/1  to  18  mg/1  as  CCL  [equivalent  to  a  decrease  from  1000 
ueq/1  to  818.2  ueq/1  (Lowenthal  and  Marias  1978)]  for  104  lower-elevation 
lakes  surveyed  by  Pennak  from  1938-1942  and  resurveyed  by  Lewis  in  1979. 
Although  Lewis  attributes  this  loss  to  acidic  deposition,  the  concurrent 
decrease  in  total  residue  suggests  that  hydrologic  variability  is  responsible 
for  this  decrease  in  alkalinity.  The  high  discharge  in  1979  (~30%  above 
average)  and  the  similar  decrease  in  total  residue  and  alkalinity  (22%  for 
alkalinity,  21%  for  total  residue)  suggests  that  this  decrease  is  mainly 
caused  by  dilution.  In  addition,  a  180  ueq/1  loss  of  alkalinity  is  much 
greater  than  any  alkalinity  loss  observed  in  acidified  regions  of  the  north- 
eastern United  States  (Hendrey  et  al.  1980).  Although  several  areas  in  the 
Rocky  Mountains  have  been  shown  to  be  potentially  sensitive  to  acidic  deposi- 
tion (Harte  et  al.,  submitted;  Dodson  1981;  Turk  and  Adams  1983),  no  con- 
vincing evidence  of  acidification  has  been  presented. 

At  present,  RMNP  is  receiving  much  less  acidic  deposition  than  acidified 
regions  in  the  northeastern  United  States.  A  comparison  of  the  (volume- 
weighted  average)  concentrations  of  the  major  ions  in  precipitation  show  that 
the  precipitation  at  RMNP  contains  much  less  acid  than  that  at  Hubbard  Brook 
(Table  16)  (NADP  1982).  During  1981  the  average  precipitation  pH  at  RMNP  was 
5.07,  compared  with  an  average  of  4.33  at  Hubbard  Brook.  Sulfate  concentra- 
tions at  Hubbard  Brook  are  also  higher  than  those  observed  at  RMNP  (48.9  ueq/1 
at  Hubbard  Brook,  compared  with  35.0  ueq/1  at  RMNP),  while  the  Cn  concentra- 
tion at  RMNP  is  more  than  double  the  concentration  at  Hubbard  Brook.  However, 
the  nitrate  concentrations  observed  at  these  two  locations  are  very  similar, 
21.9  ueq/1  at  RMNP,  compared  with  21.2  ueq/1  at  Hubbard  Brook.  These  results 
suggest  that  much  of  the  acidity  observed  in  RMNP  precipitation  may  be  caused 
by  nitrate,  while  the  sulfate  is  probably  a  combination  of  anthropogenic 
emissions  of  sulfur  dioxide,  along  with  sulfate  associated  with  Cn  from  air- 
borne dust  and  salts. 

Since  no  historical  water-chemistry  data  are  available  for  RMNP,  we  base 
our  assessment  of  its  present  acidification  status  on  the  current  composition 
of  its  waters.  According  to  electroneutral ity  conditions: 
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Table  16.   Comparison  of  Rocky  Mountain  National  Park 
and  Hubbard  Brook  precipitation  chemistry  (National 
Atmospheric  Deposition  Program  1982). 


RMNP 


Hubbard  Brook 


Vol.  wt. 

Yearly 

average 

deposition 

Ion 

(ueq/1) 

(meq/m2) 

H+ 

8.6 

3.0 

S042 

35.0 

12.2 

NO"1 

21.9 

7.5 

cr1 

4.9 

1.7 

Ca+2 

21.9 

7.7 

Mg+2 

8.1 

2.9 

K+1 

5.8 

1.1 

Na+1 

3.0 

2.1 

Vol.  wt. 

average 

(ueq/1) 

Yearly 

deposition 

(meq/m2) 

46.5 

71.1 

48.9 

75.8 

21.2 

32.5 

5.3 

8.3 

6.5 

10.4 

3.6 

5.6 

0.4 

0.7 

5.1 

7.9 

CB  =  HC031  +  S042  +  N031  +  CI  l 


H 


+1 


metals 


+n 


Under  acid  rain  conditions,  anthropogenic  acids  (H^SO,  and  HNO^)  reduce  the 
relative  concentration  of  alkalinity  in  this  equation  by  titration.  We  can 
use  the  "excess"  concentration  of  base  cations  to  provide  an  estimate  of  the 
alkalinity  replaced  by  acid  anions  (i.e.,  amount  of  acidification).  After 
correcting  for  the  deposition  of  salts  and  excluding  H  and  metals,  which  are 
negligible  at  the  pH  of  water  in  RMNP,  the  ion  balance  equation  is  reduced  to: 


C  * 
LB 


HCO"1  =  CA* 


■  so,2* 


'A 


'31 


The  results  of  these  calculations  for  the  most  sensitive  watersheds 
(Glacier  Creek,  Loch  Vale,  East  Inlet,  and  Ypsilon  Creek)  show  that  these 
waters  have,  at  the  most,  suffered  a  small  loss  of  alkalinity  (Table  17).   It 
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Table  17.  Calculation  of  acidification  status  for  Rocky  Mountain 
National  Park  watersheds  according  to  ion  balance  considerations. 


Watershed 


C* 


B 


(ueq/1) 


Alkal inity 
(ueq/1)  " 


Calculated 

Acidification 

(ueq/1) 


Roaring  River 
Ypsilon  Creek 
Upper  Fall  River 


72.3 

56.2 

155.2 


74.4 

48.0 

180.5 


None 

8.2 

None 


Andrews  Creek 
Loch  Vale 
Glacier  Creek 


57.3 
47.9 
47.6 


38.8 
41.2 
40.3 


18.5 
6.7 
7.3 


Upper  Colorado  River 


East  Inlet 


280.1 


91.5 


319.0 
85.5 


None 
6.0 


is  so  small  that  the  combined  effect  of  analytical  error  and  the  error 
associated  with  the  correction  for  atmospheric  salts  is  probably  as  large  as 
the  calculated  acidification.  The  Upper  Colorado  River,  Upper  Fall  River,  and 
Roaring  River  show  no  evidence  of  acidification. 

A  similar  method  for  evaluating  the  current  acidification  status  of  low 
ionic  strength  lakes  is  the  "predictor  nomograph"  developed  by  Henri ksen 
(1980).  This  model  was  empirically  derived  to  evaluate  the  impact  of  in- 
creased acid  precipitation  on  lakes  in  Norway  but  may  also  be  used  to  assess 
current  acidification  status.  The  model  is  based  on  electroneutral  ity  condi- 
tions, assumes  no  increase  in  weathering  of  base  cations,  and  considers 

-2 
atmospheric  deposition  to  be  the  only  source  of  SO,  .   The  nomograph  is  a  plot 

■  p  ■  p  —  y 

of  the  "excess"  Ca   and  Mg   vs.  "excess"  S0A   in  lake  waters.   "Excess" 

-2 
means  nonmarine  in  origin.   The  lake  S0A  concentrations  have  been  correlated 

-2 
with  pH  and  SO,  concentrations  in  rain  by  a  regression  analysis  of  data  from 
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719  Norwegian  lakes.  The  graph  has  been  divided  into  three  sections  to 
represent  three  stages  of  acidification.  The  first  stage  represents  waters 
that  still  contain  enough  alkalinity  to  buffer  incoming  acid  precipitation, 
"bicarbonate  lakes."  The  second  stage  of  the  nomograph  is  representative  of 
"transition  lakes."  These  waters  experience  rapid  fluctuations  in  pH  because 
of  their  low  alkalinities.  The  third  stage  of  acidification  are  the  "acid 
lakes",  characterized  by  low  pH  and  increased  aluminum  concentrations.   The 

nomograph  is  used  to  predict  acidification  status  by  plotting  precipitation  pH 

-2  +2     +2 

or  excess  lake  SO,  *  with  "excess"  Ca   +  Mg    in  ueq/1.   This  simple 

empirical  model  has  been  successfully  applied  to  lakes  in  the  Adirondack, 
Canada,  Scandinavia,  and  Scotland  (Wright  et  al.  1980).  Although  there  are 
some  differences  between  these  areas  and  RMNP  (notably  the  soils),  the  simi- 
larity of  geochemical  processes  should  allow  our  use  of  the  nomograph  with 
some  changes. 

We  have  modified  the  Henri ksen  nomograph  to  use  as  a  tool  to  evaluate  the 

+2      +2 
current  and  future  effects  of  acid  rain  in  the  Rocky  Mountains.   Ca   and  Mg 

+1     +2 
are  the  major  base  cations  in  Henriksen's  study  area,  while  Na   +  Mg   are 

the  dominant  cations  in  RMNP.   As  a  result,  we  have  replaced  the  "excess" 

+2     +2 

Ca   +  Mg   used  by  Henriksen  as  the  y_  axis  in  the  nomograph  with  Cg*.   As 

previously  stated,  RMNP  receives  minimal  sea-salt  deposition  but  does  receive 
atmospheric  deposition  of  salts  and  dust.  CR*  is  our  best  estimate  of 
"excess"  base  cations. 

Since  much  of  the  acid  deposition  in  this  region  is  in  the  form  of  HN03 

(NADP  1982;  Lewis  and  Grant  1979;  Kelley  and  Stedman  1980),  we  have  further 

-2     -1 
modi  fed  the  Henriksen  nomograph  by  using  C«*  (SO  **  +  N03  *)  as  the  x  axis  in 

the  nomograph.  The  waters  in  RMNP  contain  a  large  indeterminate  "natural,"  or 

-2 
background,  concentration  of  S0A  ,  a  result  of  the  atmospheric  deposition  of 

-2 
dust  and  salts.   We  expect  that  only  part  of  the  SO.   in  RMNP  lakes  is  the 

result  of  acid  deposition.   Aimer  et  al.  (1978)  estimate  the  background 

concentration  of  SO,,   in  Scandinavian  lakes  to  be  20-60  ueq/1.   To  estimate 

-2 
acidic  deposition,  we  will  use  Aimer's  lowest  estimate  of  background  SO-,  to 

_2  * 

calculate  a  maximum  "excess"  SO.  *: 

S0~2*  =  S0~2  -  20.0  ueq/1   . 

In  this  scenario,  CA*  =  NO"1*  +  S0~2*. 
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Plotting  this  information  on  the  nomograph  shows  that  most  of  the  lakes 
in  RMNP  can  be  classified  as  "bicarbonate",  while  a  few  approach  the 
"transition"  stage  (Figure  12). 
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Figure  12.  Status  of  23  lakes  in  RMNP. 

The  high  cation  concentrations  in  RMNP  precipitation  indicate  that  the 
park  is  undoubtedly  receiving  "natural"  deposition  of  sulfate  in  the  form  of 
salts  and  dust  from  the  surrounding  arid  regions.  However,  the  low  pH  of 
Rocky  Mountain  precipitation  indicates  that  some  of  the  sulfate  and  nitrate  is 
being  deposited  as  acid  deposition.  It  is  likely  that  the  waters  in  RMNP  have 
suffered  some  loss  of  alkalinity  because  of  acidic  deposition,  but  this  loss 
(a  maximum  of  ~10  peq/1)  is  minor  compared  with  that  experienced  in  lakes  in 
the  northeastern  United  States. 
Future  Impacts  of  Acidic  Deposition 

As  we  have  seen  in  the  previous  sections,  the  Central  Rocky  Mountains  are 
extremely  sensitive  to  acidic  deposition  but  have  as  yet  suffered  little  or  no 
acidification  on  a  regional  basis.  An  increase  in  acidic  deposition  could 
have  some  serious  consequences.  The  development  of  major  coal  and  oil  shale 
resources  upwind  of  this  region  and  the  subsequent  increase  in  acidic  deposi- 
tion could  have  a  serious  impact  on  the  surface  water  acidity  in  the  region. 
Using  the  Henriksen  predictor  nomograph  described  in  the  previous  section,  we 
can  estimate  the  best-  and  worst-case  consequences  of  an  increase  in  the 
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current  rate  of  deposition,  to  that  experienced  in  the  northeastern  United 
States.  The  scenario  for  increased  acidic  deposition  in  RMNP  assumes  a 
decrease  in  precipitation  pH  to  that  experienced  in  the  northeastern  United 
States.  This  involves  an  increase  in  H  ,  from  pH  5.1  to  4.3,  which  will  be 
accompanied  by  an  increase  in  lake  C»*  of  approximately  80  peq/1. 

In  the  worst-case  scenario,  the  increase  in  acidic  precipitation  will  not 
increase  CR*  but  will  only  result  in  decreased  alkalinity.  Using  our  estimate 
of  current  C«*,  (SO,  *  +  NO^  *  -  20  peq/1)  a  worst-case  prediction  using  the 
nomograph  indicates  that  most  of  the  lakes  in  RMNP  will  reach  "acid"  status, 
while  the  remainder  will  be  classified  as  "transition"  lakes  (Figure  13). 

In  the  best-case  scenario,  the  increase  in  acidic  deposition  will  be 
accompanied  by  an  increase  in  CR*  of  0.4  peq/1  per  1.0  peq/1  C.*.  This 
increase  in  CR*  with  increased  acidic  deposition  was  empirically  calculated  by 
Henriksen  using  data  from  low-alkalinity  lakes  in  several  areas.  The  best- 
case  scenario  predicts  that  a  few  lakes  will  remain  "bicarbonate,"  the 
majority  will  become  "transition"  lakes,  and  several  lakes  will  still  reach 
"acid"  status  (Figure  13). 

Henriksen' s  evaluation  of  data  sets  from  Norway,  Sweden,  Canada,  and  the 
United  States  shows  that  there  is  an  increase  in  base  cation  weathering  with 
acid  rain  for  some,  but  not  all  watersheds.  Increases  in  CR  will  probably 
occur  in  regions  that  have  soils  with  high  levels  of  exchangeable  bases.  When 
these  areas  receive  acidic  deposition,  base  cations  on  the  soil  exchange  sites 
will  be  replaced  with  H  by  mass  action,  increasing  surface  water  CR*.  In 
regions  with  poorly  buffered  soils,  the  increase  in  CR*  with  acid  deposition 
will  be  smaller,  since  fewer  exchange  sites  contain  base  cations  for  ion 
exchange.  Other  researchers  have  shown  that  the  rate  of  primary  mineral 
weathering  in  areas  with  granitic  and  metamorphic  rock  does  not  increase  under 
acid  rain  conditons  (Johnson  et  al.  1981).  The  low  ion-exchange  capacity  of 
the  soils  in  RMNP  and  the  bedrock  composition  indicate  that  the  increase  in 
Cg*  with  increased  acidic  deposition  in  RMNP  will  be  low  and  that  with  in- 
creasing acidic  deposition  many  of  the  lakes  will  shift  from  bicarbonate 
towards  acid  status. 
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Figure  13.  Predictor  nomographs  projecting  the  best  and  worst-case 
estimates  of  the  effects  of  increased  acidic  deposition  on  lakes  in  RMNP. 


56 


YELLOWSTONE  NATIONAL  PARK 

INTRODUCTION 

Located  in  northwestern  Wyoming  at  the  northern  edge  of  the  potentially 
sensitive  Rocky  Mountain  region,  Yellowstone  National  Park  (YNP)  is  charac- 
terized by  extreme  variability  in  geology,  geologic  history,  and  water  chem- 
istry. The  3,742  square  mile  park  contains  four  large  lakes  (Yellowstone, 
Shoshone,  Lewis,  Heart),  ahd  numerous  streams,  rivers,  and  small  back-country 
lakes,  a  number  of  which  support  trout  populations.  As  in  the  case  of  Rocky 
Mountain  National  Park,  the  general  objectives  of  the  study  in  Yellowstone 
National  Park  were  to  determine  the  sensitivity  of  waters  to  acidification  by 
current  levels  of  acid  precipitation  and  to  determine  whether  this  will  affect 
fish  populations.  The  accomplishment  of  these  objectives  is  based  primarily 
on  the  evaluation  of  a  water  chemistry  data  base  developed  over  the  last  30 
years. 

DESCRIPTION  OF  GEOLOGY 

Topographically,  Yellowstone  is  dominated  by  a  high  plateau  from  about 
1800  to  2800  m  in  altitude.  The  plateau  is  bounded  on  three  sides  by 
mountain  ranges:  the  Gallatin  and  Beartooth  mountains  to  the  north,  the 
Absaroka  range  to  the  east,  and  uplands  from  the  Teton  and  Washakie  ranges  in 
the  south  (Cox  1973).  The  continental  divide  transects  the  park  along  the 
southwestern  edge — approximately  20  percent  of  the  park  area  drains  west  into 
the  Snake  River  basin,  and  80  percent  drains  east  into  the  Missouri  River 
drainage. 

Cox  (1973)  divided  Yellowstone  into  seven  hydrologic  units  following 
geologic  and  geographic  boundaries.  These  are:  Rhyolite  plateau  (RP), 
Gallatin  (GT),  Beartooth  (BT),  Absaroka  (AB),  Falls  River  (FR),  Snake  River 
(SR),  and  West  Yellowstone  (WY)  (Figure  14).  The  plateau  area  (RP)  is  under- 
lain by  Tertiary  and  Quaternary  rhyolite  flows,  while  north  and  south  of  the 
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Figure  14.   Regional-geological  map  of  Yellowstone  National 
from  Cox  1973).   (AB)  Absaroka  region  is  predominantly  andes 
flows  and  breccia,  with  basalt,  and  some  occurrence  of  rhyol 
and  limestone;  (BT)  Beartooth  region  is  a  mix  of  Precambrian 
Paleozoic  and  Mesozoic  sandstones  and  shales,  and  Tertiary/Q 
volcanics;  (FR)  Falls  River  region  is  Quaternary  rhyol ite  an 
frequently  overlain  by  alluvial  and  glacial  deposits;  (GT)  G 
region  is  Precambrian  granites,  Paleozoic  and  Mesozoic  limes 
sandstones,  and  shales,  Tertiary/Quaternary  volcanics;  (RP) 
plateau  region  is  predominantly  Tertiary  and  Quaternary  rhyo 
(SR)  Snake  River  region  is  Paleozoic  and  Mesozoic  limestones 
and  shales,  with  some  outcroppings  of  Tertiary  rhyolite  and 
West  Yellowstone  region  is  rhyolite  overlain  by  alluvial,  gl 
lacustrine  deposits. 
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plateau  are  a  heterogeneous,  heavily  faulted  mix  of  PreCambrian  granites, 
Paleozoic  and  Mesozoic  limestones,  shales,  and  sandstones,  and  more  recent 
(Tertiary,  Quaternary)  volcanic  deposits  (GT,  BT,  SR).  Between  the  plateau 
region  and  the  Gallatin  range  along  the  extreme  western  park  border  lies  a 
series  of  deep  (c.200  feet)  alluvial  and  glacial  deposits  underlain  by 
volcanic  rocks  (WY).   East  of  the  plateau  is  a  mountainous  area  composed 
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mostly  of  Tertiary  andesitic  lava  flows  and  breccia  (AB),  while  basalt  flows 
occur  in  the  extreme  southwest  (FR). 

Rhyolite,  located  in  the  southwest  and  central  portion  of  the  park  as 
well  as  in  outcroppings  throughout  the  rest  of  the  region,  is  clearly  the 
dominant  bedrock  type.  Extremely  rich  in  silica,  rhyolite  is  chemically 
nearly  equivalent  to  granite  (Bryan  1979).  Although  fairly  uncommon  through 
the  world,  nearly  all  geysers  are  associated  with  rhyolite  formations,  since 
the  rock  provides  the  pressure-tight  "piping"  necessary  for  geyser  creation. 
The  rhyolite  in  Yellowstone  consists  primarily  of  lava  and  welded  tuff  with 
assorted  deposits  of  breccia,  ash,  and  glass.  All  of  the  types  are  chemically 
similar,  although  tuff  may  be  slightly  more  weatherable  than  lava  due  to 
differing  extrusion  and  deposition  processes  (Herzog  1982;  Cox  1973). 
Rhyolite  may  be  altered  by  hot  water  and  gas  near  hot  springs. 

Weatherabi 1 ity  of  the  bedrock  types  may  be  roughly  ranked  by  chemical 
characteristics  as  follows  (Herzog  1982  and  personal  communication;  Loughnan 
1969): 


low 


high 


granite,  rhyolite 

ash  flow  tuff 

metamorphics — biotite  gneisses  and  schists 

andesite 

basalt 

limestones,  shales,  sandstones 


In  addition  to  a  heterogeneous  geology,  the  park  bears  the  ubiquitous  remains 
of  three  major  glaciations  originating  in  the  Absaroka  Range  and  Beartooth 
Mountains.  Most  of  the  current  park  deposits  date  to  the  last  major  glacia- 
tion,  the  Pinedale  (10-30,000  years  b.p.),  which  covered  90%  of  the  park. 
Abundant  till  and  kame  deposits  throughout  Yellowstone  attest  to  the  influence 
of  these  geologic  events. 

DETERMINATION  OF  SURFACE  WATER  CHEMISTRY 

Surface  water  quality  studies  in  Yellowstone  date  to  as  early  as  the  19th 
century,  when  Gooch  and  Whitfield  (1888)  published  the  first  chemical  data  for 
Yellowstone  Lake.  Until  comparatively  recently,  however,  surveys  concentrated 
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on  the  four  largest  lakes,  all  of  which  are  chemically  insensitive  to  acidi- 
fication by  acid  precipitation  (alkalinity  >200  ueq/1).  Chemical  surveys  of 
the  more  vulnerable  back  country  lakes  were  initiated  in  1964  by  the 
Yellowstone  National  Park  Fishery  Management  Investigation,  and  continue  to 
the  present. 

Field  pH  and  alkalinity  were  measured  with  a  Hach  Field  engineers  kit 
until  1974,  when  the  Hellige  kit  was  introduced  for  pH  measurements.  During 
1969  and  1970,  laboratory  analysis  of  water  samples  were  performed  by  the 
Bureau  of  Indian  Affairs  Soil  Laboratory  in  Gallup,  New  Mexico.  (No  labora- 
tory analyses  were  performed  from  1971  to  1973.)  Beginning  in  1974  and  con- 
tinuing to  the  present,  lab  analysis  of  samples  was  performed  by  Orlando 
Laboratories,  Orlando,  Florida.  Almost  all  samples  were  refrigerated,  mailed 
within  two  days,  and  processed  within  a  week.  Analytical  methods  used  by  both 
labs  were  in  accordance  with  then-current  Standard  Methods  for  the  Analysis  of 
Water  and  Wastewater. 

All  pH  values  reported  in  this  study  were  obtained  colorimetrically  in 
the  field.  Burns  et  al.  (1981)  found  that  measurements  from  a  Hellige  kit 
agreed  to  within  0.15  pH  unit  with  potentiometric  values,  although  Pfeiffer 
and  Festa  (1980)  report  a  systematic  bias  expressed  by  the  following 
relationship: 

Old  (Hellige)  =  0.6639  (pH  meter)  +  2.534 

r2  =  0.91 

All  al kal i ni ties  reported  were  determined  in  the  laboratory  by  colori- 
metric  titration  to  pH  4.6  (Standard  Methods  #403),  except  during  the  years 
1965-1966  and  1971-1973,  when  only  field  measurements  were  taken.  Field 
al kal i  ni  ties  were  determined  colorimetrically  during  this  period  with  the  Hach 
field  engineers  kit  (Model  DR-EL  or  AL-36-P).  Colorimetric  titrations  do 
overestimate  alkalinity  by  the  amount  of  free  hydrogen  ion  in  solution  at  the 
endpoint;  at  titration  to  pH  4.6,  this  is  25  ueq/1  alkalinity  for  all  samples. 
Analytical  methods  for  major  cations  and  anions  are  given  in  Table  18  for  the 
years  1969-1970  and  1974-present.  The  current  (15th  edition)  Standard  Methods 
is  referenced  in  Table  18,  although  it  is  assumed  that  the  most  recently 
available  version  was  followed  for  historical  samples.  Quality  assurance  was 
determined  by  calculating  ion  balances  for  all  those  lakes  with  complete 
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Table  18.  Analytical  methods,  Yellowstone  National  Park. 

Parameter  Method,  1969-70  Method  1974-present 

Alkalinity     Standard  method  #403.  Colorimetric  titration  using 

phenolphthalein  (end  point  pH  8.3)  and  methyl  orange  (end 
point  pH  4.6)  indicators 

Calcium        Standard  method  #306C  Hach  Method,  p.  122 

EDTA  titri metric  Titri metric 

Magnesium      Standard  method  #313C;  calculation  from  EDTA  hardness 

Sodium        Standard  method  #325B;  flame  photometric 

Potassium      Standard  method  #322B;  flame  photometric 

Sulfate       Standard  method  #427A  EPA  #375.4 

Gravimetric  with  ignition         Turbidimetric 
residue 

Chloride       Standard  method  #407A;  argentometric 

(Standard  methods  1980;  EPA  1979;  Hach  1978) 

chemical  records  (see  Appendix  D).  Lakes  with  ion  balances  with  absolute 
values  less  than  or  equal  to  20%  are  considered  usable  for  the  purposes  of 
this  study. 

The  106  lakes  discussed  in  this  report  represent  a  significant  sample  of 
all  the  lakes  in  YNP.  All  the  major  lakes  have  been  surveyed,  including  those 
of  special  recreational  or  scientific  significance.  Every  region,  major 
geological  formation  and  geochemical  type  of  lake  has  been  sampled  in  the 
survey.  Because  of  its  large  size,  Yellowstone  Lake  was  surveyed  at  four 
different  locations  in  the  lake.  The  chemistry  of  the  lakes  not  surveyed  may 
be  inferred  by  the  chemistry  of  neighboring  surveyed  lakes. 

RESULTS  AND  DISCUSSION 

As  stated  in  the  introduction,  the  RMAS  project  relied  on  available  data 
from  Yellowstone  National  Park.  There  was  thus  no  opportunity  to  design  a 
sampling  program  to  test  hypotheses  relating  lake  and  stream  sensitivity  to 
elevation,  soils,  and  geology,  as  was  the  case  for  Rocky  Mountain  National 
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Park.  YNP  data  were  available  for  lakes  only,  and  there  had  been  no  attempt 
to  obtain  data  on  complete  watersheds  or  to  relate  sampling  locations  to 
elevation,  soils,  or  geology.  In  addition,  water  chemistry  data  were  obtained 
over  the  years  using  different  analysis  methods  and  samples  were  analyzed  by 
different  laboratories.  All  of  these  factors  made  it  impossible  to  carry  out 
the  extensive  evaluations  reported  in  the  chapter  for  Rocky  Mountain  National 
Park.  The  evaluation  of  information  from  YNP  focuses  on  sensitivity  or 
alkalinity,  which  are  the  critical  factors  in  assessing  the  health  of  fish 
populations.  In  addition,  the  general  YNP  findings  can  be  compared  with  those 
from  RMNP  to  verify  the  ability  to  extrapolate  results  from  one  region  to 
another. 
Lake  Alkalinities  (Sensitivities) 

A  spatial  sensitivity,  map,  using  alkalinity  as  an  index  of  vulnerability 
to  acidification,  was  created  for  the  park  with  recent  alkalinity  data 
(Figure  15).  The  selection  of  alkalinity  as  the  best  index  of  sensitivity  to 
acid  precipitation  is  based  on  its  physical  significance  as  the  emergent  sum 
of  many  acid-neutralizing  processes  occurring  in  the  watershed,  and  its  well- 
documented  relation  to  pH  (Henriksen  1979).  Using  Hendrey  et  al.'s  (1980) 
convention,  those  lakes  with  alkalinity  values  <200  peq/1  are  considered 
potentially  sensitive  to  long-term  inputs  of  acid  precipitation.  For 
Yellowstone,  lakes  with  alkalinity  values  reported  as  230  or  less  are  con- 
sidered in  this  category,  due  both  to  the  underestimation  implicit  in  colori- 
metric  titrations,  and  to  account  for  imprecision  in  the  analysis.  One  fourth 
(30  of  106)  the  surveyed  lakes  may  be  considered  "sensitive"  by  this  crite- 
rion. Thirteen  of  these  lakes  show  alkalinities  ^100  peq/1,  although  at  least 
six  are  influenced  by  thermal  springs  or  humic  acids  (YNP  Fishery  Management 
Investigations  reports). 

Excluding  those  lakes  that  are  located  in  the  midst  of  major  geothermal 
areas,  the  lowest  alkalinity  measured  was  40  peq/1.  Yellowstone  lakes  on  the 
average  apparently  show  considerably  greater  levels  of  alkalinity  than  those 
determined  for  Rocky  Mountain  National  Park  lakes,  probably  due  primarily  to 
differing  geology  between  the  two  areas  (most  of  the  Rocky  Mountain  National 
Park  samples  are  collected  in  areas  underlain  by  Precambrian  granites  and 
metamorphic  rocks).  Of  the  four  basins  considered  in  the  RMNP  study,  the 
Upper  Colorado  River  watershed  is  geologically  the  most  similar  to  Yellowstone, 
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Figure  15.  Alkalinity  map,  Yellowstone  National  Park. 

Yellowstone,  as  it  shows  a  significant  (  ca.  8%)  accumulation  of  volcanic 
rocks.  This  watershed  shows  the  highest  alkalinity  and  pH  values  of  the  KMNP 
study  sites  [average  alkalinity  equals  248  ueq/1  (this  study)]. 

Although  the  low-alkalinity  (<230~ueq/l)  lakes  in  Yellowstone  show  no 
clear  geographical  (spatial)  patterns,  a  geological  pattern  is  evident.  Most 
occur  within  the  large  rhyolite  flow  which  rises  from  the  southwest  and 
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spreads  along  the  central-west  and  central  portions  of  the  park.  Lakes  with 
slightly  higher  alkalinity  values  (  ca.  300  peq/1)  are  found  in  the  extreme 
southwestern  Fall  River  basalt  formation  and  the  northwestern  part  of  the 
rhyolite  plateau.  Most  of  the  northernmost  lakes,  dominated  by  andesites  and 
basalts,  show  very  high  alkalinity  (>1000  ueq/1)  and  are  not  sensitive  to  acid 
deposition.  The  exception  to  this  is  a  group  of  five  high-altitude  lakes  of 
the  Specimen  Creek  drainage  in  the  northwestern  corner  of  the  park.  These 
small  headwater  lakes  have  an  andesi te-basal t  bedrock  but  show  very  low 
alkalinities  (<200  peq/1).  The  water  supply  of  these  lakes  is  primarily  from 
snowmelt  (YNP  1965-1981))  and  lakewater  chemistry  appears  to  be  influenced 
more  by  this  dilute  source  than  by  reaction  with  the  surrounding  bedrock. 

Field  pH  measurements  were  used  to  create  a  spatial  pH  map  to  complement 
alkalinity  in  determining  geographical  distributions  of  sensitivity 
(Figure  16).  The  pH  map  for  Yellowstone  shows  most  lakes  with  circumneutral 
pH's  of  6.5-7.49,  and  most  of  the  rest  of  the  lakes  slightly  to  strongly 
alkaline.  Six  of  the  106  lakes  have  pH  values  less  than  6.5.  Of  these  six, 
two  are  unquestionably  influenced  by  thermals  and  one  probably  so,  and  three 
are  dystrophic  lakes.  As  with  the  high  alkalinity  lakes,  the  greatest  pro- 
portion of  high  pH  lakes  is  found  in  the  northern  part  of  the  park. 

A  similar  map  denoting  sulfate  concentrations  for  the  Yellowstone  lakes 
(Figure  17),  shows  a  cluster  of  high  sulfate  lakes  in  the  north  and  central 
east  areas;  these  lakes  are  generally  characterized  by  high  alkalinity  (>1600 
peq/1)  and  in  some  cases  are  found  in  low-lying  marshy  areas.  Sulfate  in 
these  lakes  probably  has  a  biological  and/or  bedrock  component  (YNP 
1965-1981).  Sulfate  concentration  in  lakes  throughout  the  rest  of  the  park 
shows  no  particular  trends;  most  lakes  show  sulfate  values  <200  peq/1,  while  a 
majority  of  lakes  in  the  rhyolite  bedrock  have  sulfate  concentrations  <100 
peq/1. 

The  alkalinity,  pH,  and  sulfate  maps  for  Yellowstone  must  be  viewed  with 
a  number  of  factors  in  mind.  First,  all  lakes  for  which  "recent"  (1965  or 
later)  data  were  available  were  used  to  create  the  maps,  regardless  of  ion 
balance.  Some  of  the  values,  therefore,  may  be  in  question  from  a  purely 
analytical  view.  Second,  there  exists  a  16-year  spread  in  chemical  analyses 
for  the  maps:  Lakes  sampled  once  in  1965  are  not  distinguished  from  those 
sampled  in  1980.   In  all  cases,  the  most  recent  chemical  measurements  are  used 
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Figure  16.  pH  map,  Yellowstone  National  Park. 

for  each  lake.  Third,  the  water  sources  for  the  lakes  vary,  and  reflect 
different  residence  times  through  different  watersheds.  Only  two  of  the  lakes 
with  alkalinities  <230  ueq/1  have  significant  inlets,  and  most  lack  outflow. 
The  primary  source  of  water  for  these  dilute  lakes  is  groundwater  and  snowmelt 
(R.  E.  Gresswell ,  personal  communication).  A  number  of  the  less  sensitive 
lakes  have  significant  flow-through.   Finally,  humic  lakes,  geothermally 
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Figure  17.  Sulfate  map,  Yellowstone  National  Park. 

influenced  lakes,  and  those  fed  by  subterranean  seeps  through  calcareous 
deposits  are  not  distinguished  from  those  lakes  which  more  readily  reflect  the 
influence  of  the  surficial  geology.  Indeed,  many  lakes  in  this  geologically 
active  area  may  be  influenced  by  unknown  factors.  Figure  18  shows  the  loca- 
tion of  the  major  thermal  springs  and  geysers  in  the  park  (Waring  1965): 
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Figure  18.  Location  of  hot  springs  and  geysers,  Yellowstone  National 
Park  (from  Waring  1965). 

springs  and  geysers  in  the  park  (Waring  1965):  These  are  concentrated  in 
geyser  basins  along  the  Firehole  and  Gibbon  rivers,  as  well  as  north  of 
Yellowstone  Lake.  Major  hydrothermal  activity  in  the  park  generally  follows 
the  southwest-northeast  sweep  of  the  rhyolite  plateau. 

Examination  of  alkalinity  and  pH  data  indicates  that  most  of  the  lakes 
are  well  buffered  against  potential  increases  in  rainfall  acidity.  One  fourth 
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of  the  surveyed  Yellowstone  lakes  do  report  alkalinity  values  <230  ueq/1 
(defined  as  "sensitive"  in  this  report),  although  a  number  of  these  (at  least 
7  of  30)  receive  significant  internal  acidity  from  humic  acids  or  hydrothermal 
springs.  All  but  four  of  the  lakes  in  the  <  230  peq/1  class  are  naturally 
barren  of  fish. 
Factors  Influencing  Alkalinity  (Sensitivity) 

Numerous  investigators  (e.g.,  Galloway  and  Cowling  1978)  have  considered 
surficial  geology  as  the  critical  factor  in  determining  water  chemistry  of 
many  lakes  and  streams.  Commonly,  surficial  geology  is  correlated  with  water 
chemistry  when  soils  are  derived  from  the  underlying  bedrock,  surface  water 
originates  as  flow  of  incident  precipitation  through  a  watershed,  and  chemical 
constituents  are  derived  primarily  through  ion  exchange/mineral  weathering 
reactions  occurring  in  the"  soil. 

As  a  whole,  lakes  in  Yellowstone  are  influenced  by  far  more  complex 
factors.  Three  major  glaciers  caused  the  deposition  of  extensive  allocthonous 
rubble  in  some  areas.  The  chemical  characteristics  of  this  transported  rock 
material  may,  in  some  cases,  dominate  lake  chemistry.  The  ash  and  lava  flows 
covering  large  areas  of  the  park  are  dotted  with  hydrothermal s  which  provide 
internal  sources  of  sulfate  and  other  chemicals.  Underground  springs, 
especially  in  the  north,  may  contribute  large  concentrations  of  dissolved  ions 
to  lakes  in  "unreactive"  bedrocks.  Finally,  watersheds  frequently  show  the 
effects  of  a  number  of  different  geologic  events,  making  simple  separation  of 
most  lakes  into  "dominant"  bedrock  type  subject  to  error. 

Considering  only  lakes  with  good  ion  balances  and  no  known  humic  or 
thermal  influence,  a  separate  variance  t-test  was  used  to  determine  if  signif- 
icant differences  existed  in  surface  water  chemistry  between  the  volcanic  and 
the  mixed-geology  regions.  The  results  of  the  pairwise  comparison  indicate 
that  the  alkalinity  and  the  base  cation  sum  are  significantly  different 
between  the  volcanic  (FR,  RP)  and  the  older  sedimentary-volcanic  (GT,  AB,  BT) 
regions  of  the  park,  with  p_  <  0.01  (Table  19).  Both  alkalinity  and  base 
cations  are  significantly  lower  in  the  volcanic  regions.  There  appears  to  be 
no  statistically  significant  difference  in  lake  sulfate  concentration  between 
the  two  regions,  indicating  that  local  sulfur  sources  are  not  confined  to  a 
single  region,  but  may  influence  lakes  in  many  areas  of  the  park. 

The  statistical  results  indicate  that,  with  the  best  available  separation 
of  "poor"  data  (bad  ionic  balance,  thermal  or  humic  influence),  lakes  in  the 
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Table  19.  Regional  alkalinity  of  Yellowstone  National  Park  lakes. 

Mean  alkalinity 
Region  ,        ,,,  No.  of  samples        Range  (ueq/1) 

RPa  182  8  80-360 

FR  320  2  320 

GT  1468  8  640-2280 

AB  1493  3  240-3520 

BT  1907  3  160-4800 

FR  =  Ralls  River  hydrologic  unit,  RP  =  Rhyolite  Plateau  hydrologic  unit,  GT  = 
Gallatin  hydrologic  unit,  AB  =  Absaroka  hydrologic  unit,  BT  =  Beartooth 
hydrologic  unit. 

Rhyolite  Plateau  and  Fall  River  basalt  region  (i.e.,  those  most  influenced  by 
volcanic  bedrock)  are  potentially  more  sensitive  to  acidification  than  those 
in  the  non-volcanic  northern  and  eastern  ranges.  The  term  "sensitive"  is,  as 
always,  relative,  for  the  more  sensitive  lakes  in  Yellowstone  show  similar 
geology  to  the  least  sensitive  lakes  in  Rocky  Mountain  National  Park  (this 
study). 

Alkalinity  is  commonly  observed  to  decrease  with  increasing  elevation,  as 
soil  depth  and  development,  and  watershed  residence  time  are  decreased.  A 
fairly  strong  alkalinity-elevation  correlation  exists  for  a  number  of  water- 
sheds in  Rocky  Mountain  National  Park,  with  only  a  weak  trend  for  the 
volcanic-andesite  dominated  Upper  Colorado  watershed  (Figure  8).  Results  for 
the  Yellowstone  lakes  dominated  by  rhyolitic  bedrock  are  intermediate 
(Figure  19),  showing  a  fairly  strong  correlation,  except  for  two  lakes  located 
in  deep  depressions  at  the  edge  of  the  formation.  The  Yellowstone  rhyolite  is 
a  somewhat  more  homogeneous  bedrock  than  the  Upper  Colorado  geology;  this  may 
explain  the  clearer  gradient. 

An  evaluation  of  the  relationship  of  sensitivity  to  basin  characteristics 
was  attempted  by  statistical  analysis  of  water  quality  data  of  lakes  from 
differing  hydrologic  regions.   Lakes  in  regions  dominated  by  volcanic  bedrocks 
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Figure  19.  Alkalinity  versus  elevation,  lakes  draining  rhyolite  bedrock, 
Yellowstone  National  Park. 


show  the  lowest  alkalinities  and  base  cation  concentrations.  On  the  whole, 
these  lakes  may  be  considered  the  most  sensitive  in  the  park  to  acidification. 
A  wide  spread  in  chemistry  values  for  the  other  lakes  indicates  that  deep 
springs  or  local  geologic  deposits  may  greatly  influence  water  chemistry. 
Historical  Changes  in  Surface  Water  Chemistry 

Acidification  trends  may  be  determined  in  two  major  ways:  Analysis  of 
historical  data,  and  use  of  (empirical  or  mechanistic)  predictive  models.  The 
major  difficulties  with  historical  data  are  that  (1)  water  chemistry  may 
change  daily  and  seasonally  in  relation  to  biological  activity,  (2)  method- 
ology may  change  over  the  years,  and  (3)  hydrological  and  meteorological 
conditions  may  not  be  comparable.  Minimization  of  daily  and  seasonal  fluctua- 
tions is  best  accomplished  by  holding  these  parameters  as  constant  as 
possible,  and  by  comparing  mean  values  of  a  number  of  lakes.  All  water 
chemistry  data  used  in  this  study  are  from  the  summer  months  of  June  through 
September,  but  daily  fluctuations  in  the  chemistry  were  not  determined,  since 
lakes  were  sampled  at  different  times  of  the  day  from  one  year  to  the  next. 
Changes  in  methodology  may  be  accounted  for  by  application  of  "correction 
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factors"  if  earlier  methods  are  systematically  biased  (e.g.,  Burns  et  al. 
1981;  Pfeiffer  and  Festa  1980).  When  complete  chemical  analyses  are  avail- 
able, ion  balances  may  serve  as  a  quality  check.  For  the  Yellowstone  data, 
all  alkalinities  were  determined  using  colorimetric  titrations,  and  all  pH's 
with  either  the  Hach  or  Hellige  kit.  Since  a  single  lab  performed  analyses  of 
Yellowstone  Park  water  since  1974  and  methods  have  not  changed  considerably, 
data  collected  since  that  year  are  highly  comparable. 

The  final  difficulty  with  historical  data  is  related  to  differing  hydro- 
logical  and  meteorological  conditions.  Ionic  concentrations  may  be  substan- 
tially reduced  by  dilution  during  a  wet  period  and  increased  by  concentration 
during  a  dry  period.  While  absolute  concentrations  have  changed,  relative 
concentrations  (with  respect  to  the  sum  of  all  ionic  constituents)  have  not. 
Figures  20  through  23  attempt  to  account  for  simple  hydrologic  variations  in 
comparing  historical  with  recent  data.  The  dotted  line  in  each  figure 
connects  the  origin  with  the  ion  sum  for  each  year.  The  position  of  each 
individual  ion  with  respect  to  that  line  indicates  whether  its  relative  con- 
centration has  increased  (above  the  line)  decreased  (below  the  line)  or 
remained  constant  over  the  period  of  time  indicated  (from  Henriksen  1982).  In 
Figures  20  through  23,  a  solid  45  degree  line  is  included  for  comparison.  By 
indicating  relative  changes  in  surface  water  concentration,  the  plots  reduce 
interpretive  error  due  to  differences  in  hydrologic  conditions.  Figures  20 
through  23  represent  available  historical  data  for  lakes  with  alkalinity 
values  of  <500  peq/1 ,  including  those  with  poor  ion  balances  (i.e.,  >20% 
variations). 

Of  the  seven  lakes  shown  in  the  historical  comparison  plots,  three  show  a 
relative  decrease  in  alkalinity  (High,  Crescent,  Grebe)  and  increase  in 
sulfate  over  the  10-year  period,  two  show  the  reverse  trend  (Wolf,  Ice),  and 
two  show  no  change  in  alkalinity  (Crag,  Cascade).  Thus,  the  overal 1  trend  for 
the  seven  lakes  appears  random.  Two  of  the  lakes  showing  alkalinity  declines 
are  located  in  the  extreme  northwestern  corner  of  the  park  in  a  group  of  five 
small  snowmel t-seepage  lakes  in  the  Specimen  Creek  watershed.  These  lakes 
show  uniformly  low  alkalinities;  however,  ion  balances  are  poor.  More 
sampling  would  be  necessary  to  make  definitive  statements  about  the  sensi- 
tivity of  these  lakes.  Of  the  three  lakes  which  show  an  increase  in  lake 
sulfate,  the  relative  sulfate  concentration  changes  (averaging  80  ueq/1)  are 
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CHEMICAL  COMPOSITION  --  1963/1979 
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CHEMICAL  COMPOSITION  --  1969/1980 
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CHEMICAL  COMPOSITION  --  1970/1979 
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Figure  22.  Historical  chemistry  comparisons,  Ice  and  Crag  Lakes.  All 
concentrations  in  peq/1. 
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CHEMICAL  COMPOSITION  —  1969/1978 
CASCADE  LAKE 
YELLOWSTONE  NATIONAL  PARK 
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1969  --  NOT  AVAILABLE 
1978 14.7 


Figure  23.  Historical  chemistry  comparisons,  Cascade  Lake.  All 
concentrations  in  ueq/1. 


considerably  greater  than  the  sulfate  load  in  precipitation  (approximately  33 
peq/1).  These  sulfate  changes  are  probably  too  great  to  be  caused  by  acid 
precipitation  and  likely  reflect  the  influence  of  local  natural  sources  of 
sulfate. 

The  relative  concentrations  of  sulfate  and  alkalinity  were  also  computed 
for  the  seven  lakes  and  are  shown  in  Table  20.  Comparing  data  in  terms  of 
relative  concentrations,  or  percent  of  each  ion  in  relation  to  the  total  ion 
sum,  is  another  method  to  minimize  the  influence  of  sampling  under  different 
hydrologic  conditions,  t-tests  on  historical  versus  recent  relative  concen- 
trations of  alkalinity  and  sulfate  for  the  seven  lakes  show  no  significant 
changes  in  these  parameters  over  time. 

Predictive  models,  such  as  that  developed  by  Henri ksen  (1979),  generally 
assume  that  bicarbonate  lost  in  water  acidified  by  atmospheric  deposition  is 
stoichiometrically  replaced  by  sulfate,  as  strong  acids  from  anthropogenic 
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Table  20.  Relative  concentration  of  alkalinity  and 
sulfate,  lakes  with  historical  data. 


Lake 

Alkal 

inity 

Sulfate 

Historical 
(1963-70) 

Recent 
(1978-80) 

Historical 
(1963-70) 

Recent 
(1978-80) 

Crescent 

47.2 

27.8 

11.3 

25.2 

High 

36.1 

20.1 

8.4 

20.9 

Cascade 

33.3 

30.0 

12.7 

16.9 

Ice 

18.6 

36.2 

28.8 

9.1 

Crag 

40.2 

41.2 

24.4 

10.7 

Wolf 

27.5 

35.6 

17.3 

12.3 

Grebe 

45.3 

35.8 

5.1 

13.2 

X 

35.5 

32.4 

15.4 

15.5 

Relative  concentration  =  [ion  of  interest]/ion  sum  x  100. 

sources  titrate  existing  alkalinity  or  replace  bicarbonate  as  a  major  weather- 
ing ion.  This  assumption  implies  a  negligible  increase  in  base  cation  release 
from  soils  with  acidification.  A  second  assumption  is  that  "internal"  natural 
sources  of  sulfate  are  negligible.  However,  Figures  20  through  23  indicate 
that  historical  changes  in  alkalinity  are  generally  balanced  by  changes  in 
both  calcium  (cation  compensation  for  acidity  increases)  and  sulfate,  in 
violation  of  the  first  assumption.  The  second  assumption  is  also  invalid  for 
Yellowstone,  which  is  heavily  influenced  by  numerous  geothermal  sources  of 
sulfur.  This  may  be  illustrated  by  comparing  the  concentrations  of  chloride 
and  sulfate  among  the  Yellowstone  lakes:  chloride  (Figure  24),  assumed  to  be 
a  conservative  ion,  remains  relatively  constant  in  nearly  all  the  Park  lakes. 
In  this  plot,  four  of  the  five  outliers  represent  large,  old  lakes  (three 
samples  from  Yellowstone  Lake  and  one  from  Lewis  Lake);  these  are  probably 
influenced  by  long-term  concentration  by  evaporation.   In  contrast,  sulfate 
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CHLORIDE  VS.   ALKALINITY 
YELLOWSTONE   NATIONAL   PARK  LAKES 
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Figure  24.  Chloride  versus  alkalinity,  Yellowstone  National  Park  lakes. 

concentration  throughout  the  Park  lakes  is  highly  variable  (Figure  25).  In 
Figure  25,  a  cluster  of  17  low  sulfate,  low  alkalinity  lakes  (mean  sulfate 
concentration  =  35.8  ueq/1)  is  evident  in  the  Park,  along  with  a  wide  scat- 
tering of  higher  sulfate  lakes.  Thus,  the  unique  characteristics  of 
Yellowstone  make  use  of  a  Henri ksen- like  predictive  acidification  model 
inappropriate. 

Detecting  historical  changes  in  surface  water  chemistry  for  Yellowstone 
is  confounded  by  natural  variability  in  water  chemistry,  internal  sources  of 
acidity,  and  a  scarcity  of  reliable  historical  data.  Yellowstone  precipita- 
tion is  currently  not  highly  acid,  nor  do  the  lakes  appear  acidified  by 
anthropogenic  activity.  No  significant  overall  change  in  water  chemistry  was 
noted  from  lakes  sampled  once,  then  resampled  several  years  later.  The  dilute 
lakes  in  the  rhyolite  plateau  and  also  in  the  far  northwestern  corner  of  the 
Park  may  be  the  most  subject  to  future  acidification. 
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SULFATE  VS.   ALKALINITY 
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Figure  25.  Sulfate  versus  alkalinity,  Yellowstone  National  Park  lakes. 

Current  Status  and  Future  Trends  in  Surface  Water  Chemistry 

Current,  historical,  and  future  acidification  scenarios  depend  upon 
changing  composition  of  precipitation.  Yellowstone  receives  an  average  of 
42  cm  of  rainfall  annually  (NOAA  1970),  but,  because  of  the  Park's  great  size 
and  varied  terrain,  precipitation  distribution  is  not  homogeneous.  The 
Continental  Divide  intercepts  eastward  moving  storms  from  the  Pacific,  causing 
the  greatest  amount  of  precipitation  near  the  Divide  in  the  southwest,  and  the 
lowest  in  the  low  altitudes.  Up  to  1.2  m  of  snow  may  accumulate  in  the  Park 
on  average  winters.  Figure  26  indicates  how  weekly  averaged  pH  and  sulfate 
concentrations  over  the  NADP  deposition-monitoring  site  vary  with  season  from 
June  1980  to  January  1982.  Although  no  overall  trends  are  evident,  chemistry 
is  quite  variable,  with  a  few  sharp  sulfate  peaks  over  the  year.  These  high 
peaks  represent  the  lowest  volume  periods  of  the  record,  and  are  not  corre- 
lated with  especially  low  pH. 

In  1981,  weighted  average  precipitation  pH  at  the  NADP  collecting  station 
was  5.2  and  weighted  average  annual  sulfate  concentration  was  33.5  peq/1. 
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PRECIPITATION  CHEMISTRY 
YELLOWSTONE  NATIONAL  PARK 
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Figure  26.   Precipitation  pH  and  sulfate  concentration,  Yellowstone  National 
Park,  WY,  1980-1981  (from  NADP  1981). 

During  the  same  year,  over  a  "typically"  acidified  area  of  the  eastern  U.S., 
Hubbard  Brook,  New  Hampshire,  precipitation  pH  averaged  4.4  and  precipitation 
sulfate  averaged  48.9  ueq/1 .  A  comparison  of  these  two  sites  indicates  that 
the  northwestern  Wyoming  area  is  currently  not  threatened  by  highly  acid 
precipitation. 

The  major  anionic  contributor  to  strong  acidity  is  the  sulfate  ion. 
Although  sulfate  mobility  can  be  affected  by  soil  adsorption  and  biological 
uptake  (Johnson  and  Cole  1980),  inputs  and  outputs  of  this  ion  have  been  found 
approximately  equal  in  a  number  of  lakes,  such  as  those  in  the  Canadian  Shield 
region  (NRCC  1981).  Long-term  increases  in  sulfate  deposition  to  sensitive 
lakes  in  Yellowstone  could  result  in  increased  sulfur  concentration  and 
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potential  acidification.  As  discussed  previously,  the  mean  sulfate  concentra- 
tion for  the  low  sul fate-low  alkalinity  group  of  lakes  was  approximately  36 
ueq/1.  This  value  corresponds  well  to  the  estimated  average  annual  precipita- 
tion input  of  about  33  ueq/1  as  derived  from  recent  NADP  data  and  may,  as  an 
extreme,  represent  a  100%  contribution  from  rainfall  (no  sulfate  derived  from 
soils,  bedrock,  or  nearby  thermal  sources).  This  simplification  represents  a 
maximum,  and  ignores  any  potential  soil  sulfate  adsorption  and  residence  time 
considerations. 

Yellowstone  does  not  currently  receive  acid  precipitation  comparable  to 
the  eastern  United  States.  The  low  sulfate  concentration  of  36  peq/1  falls  in 
the  range  of  "background"  sulfate  for  the  Canadian  Shield  lakes  (30-60  peq/1), 
and  probably  is  also  almost  entirely  of  natural  origin.  This  is  further 
corroborated  by  the  relatively  high  rainfall  pH.  In  contrast,  lakes  and 
streams  in  acidified  regions  throughout  the  U.S.  and  Canada  show  average 
sulfate  ranges  of  3  to  4  times  that  (Table  21).  The  higher  values  of  sulfate 
found  in  some  Yellowstone  lakes  is  assumed  due  to  groundwater,  soil,  bedrock, 
or  geothermal  sources.  Yellowstone  probably  experiences  little  current  acidi- 
fication from  acid  precipitation.  Considering  the  highly  variable  lake 
sulfate  concentrations,  however,  local  sulfate  sources  undoubtedly  pervade  the 
park.  Given  the  existence  of  geothermal,  local  atmospheric,  bedrock,  and 
underground  spring  sources  of  sulfate,  it  is  impossible  to  state  how  anthro- 
pogenic acidification  will  affect  these  lakes. 

Table  21.  Sulfate  levels  in  selected  acidified  waters. 

Area  Sulfate  range  (peq/1)         Reference 

Southern  Ontario  (lakes)  160-220  Dillon  et  al.  (1978) 

Nova  Scotia  (lakes)  100-140  NRCC  (1981) 

New  Hampshire  (streams)  129-142  Likens  et  al.  (1977) 


It  is  possible,  however,  to  approximate  how  Yellowstone  would  be  affected 
if  rainfall  over  the  park  were  to  change  composition  to  approximate  the 
acidity  of  typically  eastern  precipitation.  The  monitoring  station  at  Hubbard 
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Brook  experimental  forest  provides  the  longest  continuous  record  of  precipita- 
tion chemistry  in  the  United  States.  Over  the  ten  year  period  from  1963  to 
1974,  the  weighted  annual  mean  concentration  of  sulfate  in  precipitation  was 
60.3  peq/ 1 ,  and  for  nitrate,  23.7  ueq/1.  A  far  shorter  record  (1.5  years)  is 
available  for  the  evaluation  of  precipitation  chemistry  of  Yellowstone 
National  Park,  however,  a  value  of  approximately  30  ueq/1  sulfate  and  10  ueq/1 
nitrate  is  reasonable  (NADP  1981).  Thus,  a  doubling  of  the  sulfate  and 
nitrate  concentration  of  Yellowstone  precipitation  from  approximately  40  to  80 
ueq/1  would  roughly  yield  the  mean  precipitation  chemistry  for  these  constit- 
uents in  the  east.  Two  extremes  exist  in  the  response  of  watersheds  to  this 
chemical  change.  As  a  minimum,  all  of  this  new  sulfate  and  nitrate  would  be 
incorporated  into  the  watershed,  i.e.,  adsorbed  onto  soil  particles  or  tied  in 
biotic  cycling.  As  a  maximum,  all  of  the  new  strong  acid-derived  anions  would 
be  introduced  into  the  lake  ecosystem.  Thus,  the  range  of  acidification  of 
Yellowstone  lakes  would  be  from  0-40  ueq/1.  The  maximum  loss  of  alkalinity 
expected,  under  these  conditions,  is  40  peq/1 ;  the  true  loss  is  probably 
somewhat  less.  A  loss  of  40  ueq/1  alkalinity  from  each  lake  in  Yellowstone 
would  shift  the  number  of  "sensitive"  (alkalinity  <230  ueq/1 ;  no  significant 
hot  springs)  lakes  from  24  to  33.  Only  one  of  these  new  lakes  has  a  signifi- 
cant trout  population,  the  rest  are  historically  barren. 


81 


FISH  POPULATIONS 

INTRODUCTION 

The  loss  of  fish  populations  in  acidified  lakes  and  streams  has  been  the 
most  publicized  effect  of  acidic  deposition.  Although  it  is  believed  that 
other  aquatic  biota  are  also  unfavorably  affected  by  increasing  acidity, 
effects  on  fish  have  been  the  most  widely  studied  and  documented.  The  basis 
for  the  RMAS  was  to  determine  if,  in  fact,  there  was  any  evidence  that  the 
acidification  of  waters  in  the  most  sensitive  areas  of  the  Rocky  Mountain 
region  had  progressed  to  a  point  that  might  endanger  any  of  the  fish  species. 
The  evaluation  relied  on  the  acidity  and  alkalinity  data  from  RMNP  and  YNP  and 
published  and  unpublished  fisheries  information  for  the  Rocky  Mountain  region. 
Assessment  Approach 

Assessment  of  the  potential  responses  of  fish  populations  in  the  Rocky 
Mountain  region  to  acidification,  resulting  from  increased  atmospheric  strong 
acid  deposition,  was  approached  in  the  following  manner: 

1.  Areas  potentially  sensitive  to  acidification  were  defined  on  the 
basis  of  relationships  between  water  quality,  geology,  and  soils, 
developed  in  representative  watersheds  in  RMNP  and  YNP. 

2.  Information  on  predominant  fish  species,  fisheries,  and  management 
policies  prevailing  in  potentially  sensitive  areas  were  obtained 
from  the  published  literature  and  unpublished  fishery  survey  data. 

3.  The  possible  effects  of  acidification  on  the  species  at  risk  were 
evaluated  from  published  information  on  life  history  patterns,  in 
relation  to  expected  changes  in  water  quality.  Potentially  critical 
life  history  stages  were  identified  on  the  basis  of  documented 
effects  on  species  with  comparable  life  histories  and  habitat  re- 
quirements, in  regions  where  acidification  has  already  occurred. 

4.  The  implications  of  possible  acidification  effects  on  fish  popula- 
tions for  State  and  Federal  management  policies  were  considered, 
critical  information  gaps  were  identified,  and  future  research 
direction  recommended. 


82 


Evaluations  in  this  study  of  water  chemistry  data  obtained  from  RMNP  and 
YNP  and  data  from  other  studies  in  the  Rocky  Mountains  have  not  revealed  any 
instances  of  chronic  acidification  at  levels  that  would  be  detrimental  to  fish 
survival.  Chronic  acidification  of  surface  water  occurs  where  rates  of  base 
supply,  derived  from  weathering  processes  in  the  drainage  basins,  are  exceeded 
by  rates  of  strong  acid  input  from  atmospheric  depositions.  Such  conditions 
have  been  observed  in  sensitive  areas  of  eastern  North  America  and  Scandinavia 
(Wright  et  al.  1980),  where  chronically  acidified  lakes  and  streams  exhibit  an 
absence  of  bicarbonate  buffering,  low  pH  (<5.0),  and  increased  levels  of 
potentially  toxic  metals.  Fish  populations  are  generally  absent  or  extremely 
stressed  at  these  chronic  levels  of  acidification  (Schofield  1976;  Muniz  and 
Leivestad  1980).  However,  it  is  important  to  recognize  that  perturbations 
leading  to  the  observed  decimation  of  fish  populations  in  these  chronically 
acidified  waters  probably  were  initiated  at  much  earlier,  transitional  stages 
of  acidification,  prior  to  the  complete  loss  of  bicarbonate  buffering  in  the 
systems  (Henri ksen  1980;  Schofield  1982). 

As  indicated  in  the  previous  sections  on  water  chemistry  and  geology,  a 
large  number  of  headwater  drainages  in  the  geologically  sensitive  areas  of  the 
Rocky  Mountains  exhibit  extremely  low  rates  of  base  supply,  as  indicated  by 
the  low  alkalinity  and  cation  levels.  These  drainage  systems  border  on  a 
transitional  stage  of  acidification,  at  current  levels  of  atmospheric  strong 
acid  deposition.  Even  minor  excursions  of  increased  strong  acid  loading,  as 
might  occur  during  the  early  stages  of  snow-melt,  could  episodically  acidify 
these  systems.  The  presently  available  water  chemistry  data  represent  late 
snow-melt  or  baseflow  conditions,  hence  the  contemporary  occurrence  of  such 
events  is  unknown.  The  probable  consequences  of  such  perturbations  for 
indigenous  fish  populations,  should  they  occur  in  these  sensitive  drainage 
systems,  are  considered  in  subsequent  discussions  of  fish  responses  to  water 
quality  change. 

RESULTS  AND  DISCUSSION 

Fish  Resources  in  Potentially  Sensitive  Areas  of  the  Rocky  Mountains 

The  general  areas  in  the  Rocky  Mountains  potentially  sensitive  to  acidi- 
fication are  situated  in  the  alpine  (>3,200  m)  and  upper  montane  (2,500-3,200 
m)  altitudinal  zones.   Sensitivity  of  specific  drainage  basins  is  tempered  by 
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local  variations  in  bedrock  and  surficial  geology,  as  discussed  previously. 
Fish  species  richness  and  standing  crops  are  relatively  low  in  these  headwater 
drainages  and  many  were  originally  devoid  of  native  fish  populations,  because 
of  natural  barriers  to  colonization  in  the  form  of  falls  and  stream  gradients 
impassible  to  fish.  However,  most  of  the  originally  barren  headwater  lakes 
were  subsequently  stocked  with  either  native  or  exotic  salmonids,  beginning  in 
the  early  1880s  (Pennak  1963).  The  cutthroat  trout  (Salmo  clarki )  is  the  only 
native  trout  in  the  area  and  it  is  represented  by  several  subspecies  re- 
stricted to  specific  drainage  systems.  The  only  other  native  species  locally 
present  in  these  sensitive  watersheds  are  the  grayling  (Thymal lus  signifer) , 
several  species  of  the  mountain  sucker  (Pantosteus  spp.),  sculpin  (Cottus 
spp.),  and  the  mountain  whitefish  (Prosopium  wi 1 1 iamsoni ).  Introduced,  non- 
native  species  that  have  become  regionally  established  include  the  brook  trout 
(Salvel inus  fontinal is) ,  rainbow  trout  (Salmo  gairdneri ) ,  and  brown  trout 
(Salmo  trutta).  The  intentional  and  inadvertent  introduction  of  these 
species,  particularly  the  brook  trout  and  rainbow  trout,  has  been  most  detri- 
mental to  the  native  cutthroat  trout  stocks,  both  because  of  competitive 
exclusion  and  hybridization.  Additionally,  indiscriminate  hatchery  plantings 
of  cutthroat  have  so  blurred  the  genetic  integrity  of  this  species,  that  only 
in  a  few  high  lakes  and  streams  can  the  original  subspecies  still  be  distin- 
guished morphologically  (Pennak  1963).  The  problems  associated  with  the 
preservation  and  maintenance  of  these  subspecies  will  be  discussed  further  in 
the  section  on  wild  trout  management,  with  specific  reference  to  the  water- 
sheds in  RMNP. 

The  extensive  fisheries  data  base  for  YNP  provides  both  a  detailed  chron- 
ology of  the  changes  in  fish  distribution  and  management  policies  for  this 
region  since  the  1800s.  The  comprehensive  summary  of  fish  stocking  activities 
in  YNP  by  Varley  (1981)  provides  the  most  recent,  detailed  information  on  the 
status  of  fish  populations  in  the  Park.  An  earlier  publication  by  Fromm 
(1941)  provides  a  more  anecdotal  account  of  the  early  hstory  of  fish  surveys. 
Although  eighteen  species  of  fish  (12  native  species)  are  currently  recorded 
as  being  present  in  YNP  waters,  their  distribution  is  quite  limited.  Many  of 
the  lakes  and  streams  in  YNP  were  and  still  are  fishless,  primarily  because  of 
the  physical  limitations  to  colonization  described  earlier.  Lakes  and  streams 
in  approximately  40%  of  YNP  (Figure  27)  were  estimated  to  be  devoid  of  fish  in 
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Figure  27.  Cross-hatching  shows  the  area  in  Yellowstone  National  Park 
found  barren  of  fishes  by  Jordan  in  1889,  with  the  exception  of  sculpins 
in  the  Gibbon  River  above  Gibbon  Falls  (from  Fromm  1941). 

the  1800 ' s  (Jordan  1889).  Early  fishery  management  efforts  sought  to  estab- 
lish fishes  in  these  barren  waters  through  extensive  plantings  of  both  native 
and  exotic  species.  Varley  (1981)  provides  a  detailed  chronology  of  these 
stocking  efforts.  Many  of  the  previously  fishless  headwater  streams  and  lakes 
received  introductory  stockings  during  the  period  1920-1935.  Although  fish 
populations  were  probably  established  temporarily  in  most  of  these  waters  by 
stocking  the  majority  of  the  initially  fishless  lakes  are  presently  barren 
once  again.  However,  this  is  not  the  case  for  streams.  Most  of  these 
historically  fishless  streams  are  currently  supporting  viable  fish  populations 
(Table  22).  The  difference  probably  reflects  the  lack  of  suitable  spawning 
habitat  in  many  of  the  small  headwater  lakes.  In  the  previous  section  23 
potentially  sensitive  lakes  were  identified  in  YNP  (alkalinity  <200  ueq/1) 
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Table  22.  Current  and  historical  fish  population  status  of  lakes 
and  streams  in  Yellowstone  National  Park  (from  Varley  1981). 


Category 


Number 


Historical 

Current 

Fishless 

Fishless 

Fishless 

Fish 

Fish 

Fish 

Unknown 

Total 

Lakes 


Streams 


29 

18 
11 
13 
70 


2 
38 
17 
16 
73 


from  the  water  chemistry  data  available  for  107  lakes.  Information  on  histor- 
ical and  current  fish  status  obtained  for  11  of  the  sensitive  lakes  is  sum- 
marized in  Table  23.  Only  one  of  the  lakes  supported  cutthroat  trout  histor- 
ically and  the  only  other  species  currently  present  is  the  introduced  brook 
trout  in  three  of  these  lakes. 
Fisheries  Management  Policies  in  Potentially  Sensitive  Areas 

State  and  federal  policies  which  must  be  considered  in  assessing  poten- 
tial acidification  impacts  on  fish  populations  include  state  level  wild  trout 
management  policies,  National  Park  Service  aquatic  resource  management  policy, 
and  U.S.  Fish  and  Wildlife  Service  implementation  of  the  Endangered  Species 
Act,  including  cooperative  agency  efforts. 

The  NPS  policies  governing  the  management  of  aquatic  ecosystems  in  desig- 
nated Natural  Zones  are  particularly  relevant  to  the  sensitive  areas  identi- 
fied in  RMNP  and  YNP.  The  primary  goal  of  resource  management  programs  in 
natural  zones  is  the  preservation  and  restoration  of  native  aquatic  eco- 
systems, including  those  waters  originally  barren  of  fish.  Park  waters 
falling  into  this  category  are  allowed  (or  in  some  cases  "rehabilitated")  to 
revert  to  their  original  fishless  condition.  This  would  include  the  majority 
of  sensitive  waters  in  YNP  and  a  large  number  in  RMNP  as  well.  Strict  inter- 
pretation and  adherence  to  this  policy  would  seem  to  make  the  question  of 
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Table  23.  Total  alkalinity  and  fish  population  status  for  sensitive 
lakes  in  Yellowstone  National  Park  (from  this  study,  and  Varley  1981). 


Alkalinity 
(Meq/1)  " 

Fish  Status 

Lake  name 

Historical 

Current 

Wrangler 

40 

Fishless 

Fishless 

Summit 

60 

Fishless 

Fishless 

Shelf 

80 

Fishless 

Fishless 

Mt.  Everts 

160 

Fishless 

Fishless 

Ice 

200 

Fishless 

Fishless 

Ranger 

160 

Fishless 

a 

Obsidian 

80 

Fishless 

brook  trout 

High 

170 

Fishless 

cutthroat 

Forest 

192 

Fishless 

cutthroat3 

Trilobite 

200 

Fishless 

brook  trout 

Robinson 

100 

cutthroat 

brook  trout 

Status  questionable 

potential  acidification  one  of  largely  academic  concern,  with  little  relevance 
for  fisheries  programs.  However,  there  are  two  important  exceptions  to  this 
policy  which  should  prioritize  concern  for  the  potential  acidification  of 
these  waters.  In  some  previously  fishless  waters,  populations  of  either 
native  or  non-native  fishes  have  become  well  established.  If  these  estab- 
lished populations  represent  distinctive  sub-species  (particularly  of  endan- 
gered or  threatened  indigenous  species)  or  valuable  genotypes  of  non-native 
species,  the  populations  may  be  designated  as  "naturalized"  and  managed  as 
integral  components  of  the  ecosystem.  Similarly,  should  research  indicate  the 
need  to  stock  or  re-locate  threatened  or  endangered  species,  fishless  or 
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formerly  fishless  waters  may  be  utilized  as  refugia  for  these  species  (Watson 
1980;  Jones  1980). 

The  relatively  recent  evolution  and  implementation  of  "Wild  trout  manage- 
ment" programs  (most  notably,  the  Colorado  Division  of  Wildlife  Commission's 
Wild  Trout  Policy,  adopted  in  1982)  in  the  Rocky  Mountain  states  is  indicative 
of  the  widespread  concern  for  the  depletion  of  wild  trout  stocks.  The  objec- 
tives of  these  programs  are  diverse  and  although  they  are  user  oriented,  they 
share  a  common  objective  with  NPS  and  FWS  programs  for  the  protection  of 
native  species.  Cooperative  efforts  are  particularly  evident  and  vital  in 
programs  developed  for  the  protection  and  restoration  of  endangered  or 
threatened  species  (Johnson  and  Rinne  1982;  Behnke  and  Zorn  1976).  Genetic 
swamping,  habitat  destruction,  and  competition  with  introduced  species  are 
primary  reasons  for  the  widespread  demise  of  native  trout  in  the  Rocky 
Mountain  Area  (Behnke  1979).  The  Endangered  Species  Act  provides  a  vehicle 
for  the  rastoration  of  these  western  salmonids  and  recovery  action  programs 
are  currently  in  place  throughout  the  region.  Implementation  of  these 
programs  first  requires  listing  of  the  species  as  endangered  or  threatened, 
which  involves  extensive  study  and  documentation  (Johnson  and  Rinne  1982).  In 
addition  to  affording  the  listed  species  protection,  under  the  conditions  of 
the  Act,  habitats  necessary  for  their  survival  are  also  afforded  protection 
from  adverse  alteration  or  destruction.  How  this  condition  might  apply  to 
federal  programs  that  influence  the  potential  for  acidification  of  these 
habitats  is  unclear  at  this  time.  The  general  components  of  recovery  plans 
entail  definition  of  the  genetic  purity  of  the  population(s) ,  survey  and 
monitoring  of  existing  populations,  evaluation  of  habitat  in  existing  and 
candidate  reintroduction  waters,  and  reestabl ishment  of  the  native 
populations. 

Endangered  or  threatened  Salmo  species  present  in  potentially  acid  sensi- 
tive regions  of  the  Rocky  Mountains  include  several  subspecies  of  the  cut- 
throat trout.  These  fish  and  their  listings  are  given  in  Table  24.  The  life 
history  patterns  of  these  species  and  the  potential  problems  that  acidifica- 
tion might  impose  on  restoration  programs  in  sensitive  watersheds  are  consid- 
ered in  the  following  section. 
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Table  24.  Endangered  and  threatened  Salmo  genotypes  in  the 
Rocky  Mountain  Region. 


Common  name 


Scientific  name 


States 


Listing  status 


Greenback  cutthroat 
Colorado  cutthroat 

Westslope  cutthroat 
Eastslope  cutthroat 
Snake  River  cutthroat 
Rio  Grande  cutthroat 

Gila  trout 
Arizona  trout3 


Salmo  clarki  stomias 


Colorado 


Salmo  clarki  pleuriticus  Colorado, 

Wyoming 


Salmo  clarki  subsp. 
Salmo  clarki  subsp. 
Salmo  clarki  subsp. 


Salmo  gilae 
Salmo  apache 


Montana 
Montana 
Wyoming 


Salmo  clarki  virginal  is   Colorado, 

New  Mexico 


New  Mexico 
Arizona 


Endangered 
Threatened 

Endangered 
Endangered 
Rare 
Endangered 

Endangered 
Threatened 


Indicates  nationally  recognized  as  endangered. 

Potential  Impacts  of  Acidification  on  Fish  Populations 

The  few  remaining  native  trout  populations  in  interior  regions  of  the 
Rocky  Mountains  persist  only  in  small,  isolated  headwater  drainages.  These 
native  cutthroat  populations  prefer  and  function  best  at  lower  temperatures 
than  other  species.  With  few  exceptions,  the  cutthroat  coexist  and  dominate 
introduced  species  only  in  these  cold,  headwater  situations  (Behnke  1979). 
Unfortunately,  the  displacement  of  these  rare  and  endangered  genotypes  to 
headwater  drainages  also  makes  them  most  susceptible  to  potential  acidifica- 
tion in  these  sensitive  habitats.  Because  of  the  current  endangered  status  of 
these  fish,  the  Salmo  clarki  complex  must  be  considered  the  primary  species  at 
risk  and  priority  is  given  here  to  a  consideration  of  potential  acidification 
impacts  on  these  populations. 

Given  the  restricted  distributions  of  the  subspecific  cutthroat  popula- 
tions and  the  geologic  heterogeneity  in  watershed  sensitivities  to  acidifica- 
tion in  the  Rocky  Mountains,  current  juxtaposition  of  specific  populations 
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significantly  determines  their  relative  sensitivity  to  any  regional  increase 
in  acidic  deposition.  The  risk  of  possible  extinction,  resulting  from  acidi- 
fication, could  certainly  be  lessened  if  management  recovery  teams  had  prior 
knowledge  of  the  distribution  of  the  populations  relative  to  specific  water- 
shed sensitivities  to  acidification.  In  the  process  of  selecting  new  sites 
for  reestabl ishment  of  cutthroat  populations,  watershed  acid  neutralizing 
capacity  should  also  be  considered  among  the  selection  criteria.  The  current 
distributions  of  the  Colorado  River  cutthroat  (S.  c.  pleuriticus)  and  the 
greenback  cutthroat  (S.  c.  stomias) ,  relative  to  the  marked  differences  in 
watershed  sensitivity  observed  in  RMNP  exemplify  this  point.  The  Colorado 
River  trout  is  native  to  the  upper  Colorado  River  basin  (Figure  28)  and 
efforts  are  being  made  to  reestablish  this  fish  in  headwater  sections  located 
in  RMNP.  The  greenback  occupies  the  headwaters  of  the  Arkansas  and  South 
Platte  drainages  (Figure  29),  on  the  eastern  side  of  the  Continental  Divide. 
Comparable  efforts  have  been  made  to  reintroduce  this  form  into  these  areas  of 
the  Park.  The  higher  alkalinities  in  the  upper  Colorado  drainage  would  cer- 
tainly favor  the  maintenance  of  S.  c.  pleuriticus  in  the  advent  of  increased 
acid  deposition,  whereas  the  low  alkalinities  on  the  other  side  of  the  Divide 
would  probably  not  inhibit  acidification  of  the  greenback's  prime  habitat. 
Behnke's  and  Zorn's  (1976)  prophetic  suggestion  that  the  greenback  trout  may 
be  the  most  vulnerable  of  all  western  trouts  to  extinction,  would  likely  be 
realized  with  acidification. 

In  addition  to  distribution,  there  are  species  specific  life  history 
characteristics  that  must  be  considered  in  assessing  potential  sensitivity  to 
acidification.  All  of  the  interior  western  trouts  of  the  genus  Salmo  have 
basically  similar  life  histories.  They  spawn  in  the  spring  when  water  temper- 
atures reach  5.5-9.0°C,  which  can  be  anywhere  from  early  April  to  June  or 
July,  depending  on  latitude  and  elevation.  All  are  obligatory  stream  spawners 
and  fry  emergence  occurs  in  early  to  mid-summer.  Growth,  maturation,  and 
fecundity  are  variable,  depending  on  prevailing  temperature  regimes  and  pro- 
ductivity of  the  local  habitats  (Behnke  and  Zorn  1976).  The  basic  life  his- 
tory pattern  outlined  above  is  markedly  different  from  that  exhibited  by 
salmonid  populations  inhabiting  waters  of  eastern  North  America  and  Scandi- 
navia, where  acidification  impacts  have  been  described  (Schofield  1976;  Muniz 
and  Leivestad  1980).   The  predominant  salmonid  species  in  these  areas 
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Figure  28   Indigenous  distribution  of  Salmo  clarki  pleuriticu* 
(from  Behnke  and  Zorn  1976).         LaLLL  H      cu' 


(SalveVmus  fontinalis,  Salvelinus  namaycush.  Salmo  trutta.  and  Sajmo  saUr) 
are  all  fall  spawners  and  either  stream  or  lake  spawning  may  be  locally  preva- 
lent. Fry  emergence  for  these  species  is  in  early  spring,  often  coinciding 
with  snowmelt  periods  when  water  quality  is  very  poor  in  acidified  areas  The 
fry  (part,cularly  during  and  shortly  after  hatching)  of  most  of  these  species 
have  also  been  found  to  be  physiologically  more  sensitive  to  acidification 
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Figure  29.  Indigenous  distribution  of  Salmo  clarld  stpjnias 
(from  Behnke  and  Zorn  1976). 

than  either  e.bryos  or  older  fish  (Spry  et  al.  1981;  Baker  and  Schofield 
1980)  For  these  reasons,  the  early  life  history  stages  of  develop.ent  have 
been  identified  as  critical  periods  for  survival  in  acidified  habitats.  Given 
the  different  life  history  pattern  of  the  western  Salmo  sp.  and  uncerta, 
water  quality  conditions  that  .ight  prevail  during  the  early  life  tnstory  of 
these  fish,  under  an  acidification  regime,  it  is  difficult  to  extrapolate 
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these  findings  from  eastern  North  America  and  Scandinavia  to  the  Rocky 
Mountain  region.  Water  quality  conditions  would  be  more  favorable  for  cut- 
throat fry  if  episodic  acidification  were  to  occur  before  hatching. 

A  more  likely  critical  period  for  the  headwater  populations  of  this 
species  might  be  during  late  winter  and  early  spring,  prior  to  or  during 
spawning.  Studies  by  Johnson  and  Webster  (1977)  demonstrated  a  marked  avoid- 
ance of  acidic  water  by  spawning  brook  trout  and  Flick  et  al.  (1982)  noted  a 
tendency  for  emigration  of  brook  trout  populations  from  lakes  during  episodes 
of  acidic  snowmelt.  These  latter  observations  are  particularly  relevant  in 
terms  of  the  known  propensity  for  emigration  by  cutthroat  populations  inhabit- 
ing headwater  streams  subject  to  unfavorable  winter  temperature  extremes 
(Bjornn  1971).  According  to  Behnke  and  Zorn  (1976),  this  was  also  the  cause 
of  a  failed  transplant  of  greenback  trout  in  the  North  Big  Thompson  River  in 
RMNP,  where  all  the  fish  migrated  downstream  over  a  barrier  during  the  winter 
months.  A  dense  brook  trout  population  below  the  barrier  made  it  doubtful 
that  the  cutthroat  population  could  sustain  itself  there.  These  behavioral 
responses  to  adverse  environmental  conditions  suggest  a  subtle,  but  potential- 
ly devastating  mechanism  whereby  even  relatively  minor  acidification  excur- 
sions in  headwater  trout  refugia  could  lead  to  population  extermination. 
Behavioral  studies  would  be  needed  to  define  thresholds  of  acidification  that 
induce  avoidance  in  the  form  of  downstream  migration.  It  is  quite  likely  that 
these  thresholds  (e.g.,  in  terms  of  pH  change)  would  be  much  lower  than  those 
determined  by  classical  bioassay  for  definition  of  dose-response  functions, 
where  response  is  death  or  acute  physiological  stress. 

As  noted  above,  the  relevance  of  any  discussion  of  species  specific 
dose-response  functions,  as  usually  defined,  is  somewhat  questionable  at  this 
point,  given  the  uncertainties  in  determining  critical  life  history  stages  and 
population/community  level  responses.  However,  there  are  some  potentially 
important,  physiological  level  questions  that  need  to  be  considered  in  defin- 
ing the  sensitivity  of  this  species  to  acidification.  Most  comparisons  of 
relative  tolerance  to  acidity  among  salmonid  species  indicate  that  the  rainbow 
trout  (which  is  closely  related  to  the  cutthroat)  is  the  most  sensitive  to  low 
pH  (Haines  1981).  However,  no  studies  have  yet  been  conducted  with  the  cut- 
throat to  define  its  relative  tolerance.  Intraspecific  variation  in  acid 
tolerance  to  acidity  has  been  observed  in  brook  trout  (Flick  et  al.  1982)  and 
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brown  trout  (Gjedrem  1976),  and  the  marked  polytypic  character  of  S.  clarki 
(Trojnar  and  Behnke  1974)  suggests  that  it  might  also  exhibit  significant 
variation  in  acid  tolerance  among  the  defined  subspecies.  Determination  of 
the  extent  of  this  variation  in  acid  tolerance  among  the  extant  populations  in 
the  Rocky  Mountain  region  would  be  potentially  useful  for  rehabitation  pro- 
grams, in  the  event  of  acidification. 

The  potential  for  aluminum  mobilization  by  acidification  in  Rocky 
Mountain  watersheds  is  uncertain,  but  probably  not  as  great  as  in  eastern 
forested  systems  where  soils  exhibit  marked  accumulations  of  amorphous  forms 
of  aluminum  that  are  readily  mobilized  by  acidic  deposition  (Cronan  and 
Schofield  1979).  The  enhanced  toxicity  of  acidified  waters  containing  alum- 
inum is  well  documented  (Schofield  and  Trojnar  1980)  and  if  aluminum  mobili- 
zation is  not  an  integral  facet  of  the  acidification  process  in  Rocky  Mountain 
soils,  then  comparisons  of  biological  responses  between  the  two  areas  would  be 
further  complicated  by  significant  differences  in  solution  chemistry. 

Another  important  difference  in  solution  chemistry  between  Rocky  Mountain 
waters  and  those  of  eastern  North  America  are  the  levels  of  dissolved  calcium. 
The  low  calcium  levels  in  the  headwaters  of  the  Rockies  are  more  similar  to 
the  alpine  waters  of  Scandinavia,  than  those  of  eastern  North  America.  The 
significance  of  calcium  as  a  mediator  of  gill  membrane  permeability  and  acid 
stress  was  noted  by  Brown  (1981).  Below  1  mg/1  of  calcium  (typical  for  many 
of  the  Rocky  Mountain  headwaters)  trout  would  be  susceptible  to  acid  induced 
osmoregulatory  stress  at  much  higher  pH  levels  than  populations  inhabiting 
higher  (>2  mg/1,  typical  for  eastern  waters)  calcium  waters.  However,  poten- 
tial adaptation  (acclimation)  of  resident  trout  populations  to  low  calcium 
environments  might  ameliorate  expected  acid  stress  responses  (Guthrie  1981). 
Again,  it  is  difficult  to  generalize  dose-response  function  when  dealing  with 
fish  populations  that  have  evolved  under  unique  environmental  conditions.  The 
intraspecif ic  genetic  diversity  of  Salmo  clarki  is  quite  remarkable,  but  the 
potential  adaptability  of  the  genotypes  of  acidified  environments  remains  to 
be  determined. 
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CONCLUSIONS 


INTRODUCTION 

This  section  summarizes  the  conclusions  reached  concerning  the  four  major 
objectives  outlined  in  the  introduction.  The  first  two  sections  are  specific 
to  the  two  areas  studied,  RMNP  and  YNP.  The  third  section  evaluates  the 
results  of  the  water  chemistry  studies  in  terms  of  effect  on  several  species 
of  fish  important  to  this  region.  Finally,  the  assessment  of  sensitivity  of 
lakes  and  streams  to  acidification  in  the  Rocky  Mountain  Region  in  general  is 
based  on  studies  in  each  of  the  two  Parks. 

ROCKY  MOUNTAIN  NATIONAL  PARK 
Soils  and  Geology 

1.  The  soils  of  the  Rocky  Mountain  National  Park  at  the  higher  eleva- 
tions (above  9,000  feet)  are  only  slightly  developed.  They  have 
surface  organic  horizons  in  some  cases  and  usually  have  a  darkened 
A-horizon  underlain  by  slightly  weathered  material  derived  from 
granite,  diorites  and  other  rocks  low  in  bases.  The  soils  are 
coarse,  low  in  clay,  low  in  base  cations  and  relatively  acid  (modal 
pH  5.2).  The  organic  matter  provides  most  of  the  cation  exchange 
capacity  and  there  is  no  evidence  of  sulfate  adsorption  capacity. 

2.  The  physical  properties  of  these  soils — coarse,  high  in  sand,  low  in 
clay,  and  steeply  sloping — promote  rapid  movement  of  water  to  the 
streams  and  lakes. 

3.  The  low  base  saturation,  low  pH  and  low  cation-exchange  capacity 
provide  little  opportunity  for  neutralizing  the  acidic  inputs  or 
providing  significant  alkalinity  to  the  water  system. 

4.  The  soil  itself  is  probably  resistant  to  rapid  acidification  due  to 
the  abundance  of  relatively  unweathered  minerals. 

Surface  Water  Chemistry 

1.  Primary  mineral  weathering  appears  to  be  the  dominant  mechanism 
determining  the  concentrations  of  base  cations,  silica  and  alka- 
linity throughout  the  Park.   Factor  analysis  shows  that  the  primary 
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mineral  weathering  of  oligoclase  and  biotite  accounts  for  almost  50% 
of  the  variance  in  stream  chemistry  in  areas  underlain  by  granite 
and  biotite  gneiss  and  schist.  In  areas  that  contain  tertiary 
intrusive  bedrock,  the  mineral  weathering  of  mafic  materials,  sulfur 
bearing  minerals,  and  oligioclase  account  for  more  than  50%  of  the 
variance  in  stream  chemistry. 

2.  Atmospheric  deposition  is  the  primary  source  of  chloride  and  nitrate 
in  the  streams  of  the  Park.  Atmospheric  deposition  is  also  the 
primary  source  of  sulfate  in  the  Park  waters,  with  the  exception  of 
the  Upper  Colorado  and  Upper  Fall  River  basins  where  the  weathering 
of  sulfur  bearing  minerals  is  a  source  of  sulfate.  Chloride  and 
sulfate  are  relatively  constant  with  elevation,  while  nitrate  con- 
centrations are  highest  above  the  timberline,  where  biological 
activity  is  lowest.  Atmospheric  deposition  is  also  a  significant 
source  of  streamwater  CR,  a  result  of  deposition  of  airborne  dust 
and  salts  from  the  dry,  windy,  regions  upwind  of  the  Park. 

3.  Most  of  the  waters  in  the  study  watershed  have  alkalinities  of  <100 
peq/1.  The  alkalinities  are  lowest  at  higher  elevations.  The 
waters  of  Glacier  Gorge  and  Ypsilon  Creek  are  extremely  sensitive. 
Those  of  Roaring  River  and  Upper  East  Inlet  watersheds  are  very 
sensitive  with  alkalinities  of  <100  ueq/1.  The  waters  of  Lower  East 
Inlet  and  the  higher  elevations  of  Upper  Fall  River  watersheds  are 
slightly  higher  with  alkalinities  <200  ueq/1.  The  waters  in  the 
lower  elevations  of  Upper  Fall  River  watershed  are  non-sensitive 
with  alkalinities  >200  ueq/1. 

4.  Based  on  the  current  concentrations  of  sulfate,  nitrate,  base 
cations  and  alkalinity  we  estimate  that  the  waters  of  Glacier  Gorge, 
Loch  Vale,  East  Inlet  and  Ypsilon  Creek  watersheds  may  have  suffered 
a  small  loss  of  alkalinity  (<10  peq/1).  We  estimate  that  the  waters 
in  the  Upper  Colorado,  Upper  Fall  River  and  Roaring  River  watersheds 
have  suffered  no  loss  of  alkalinity. 

5.  If  precipitation  in  Rocky  Mountain  National  Park  becomes  as  acidic 
as  in  the  eastern  United  States,  we  estimate  that  in  the  worst  case 
(i.e.,  if  the  increase  in  acidic  deposition  does  not  cause  an 
increase  in  base  cation  loss  from  the  watershed),  most  of  the  lakes 
in  Rocky  Mountain  National  Park  will  become  acidified  below  pH  4.7. 
At  the  other  extreme,  if  for  every  peq/1  increase  in  acid  sulfate  in 
the  waters  there  is  a  0.4  peq/1  increase  of  base  cation  in  the 
waters,  we  estimate  using  the  Henri ksen  nomograph  that  the  majority 
of  lakes  will  become  transitional  and  only  a  few  will  reach  acid 
status.  The  low  ion-exchange  capacity  of  the  soils  in  the  Park  and 
the  resistance  of  the  bedrock  to  chemical  weathering  indicates  that 
the  increase  in  base  cation  concentrations  with  increased  acidic 
deposition  in  the  Park  will  be  low  and  that  with  increasing  acidic 
deposition  many  of  the  lakes  will  shift  from  a  bicarbonate  towards 
an  acid  status. 


96 


YELLOWSTONE  NATIONAL  PARK 

1.  The  region-chemistry  t- tests  and  alkalinity-elevation  plot  showed 
that,  after  other  influences  are  accounted  for,  surficial  geology 
and  elevation  are  correlated  to  water  chemistry.  Rhyolite,  the 
dominant  bedrock  formation  in  the  park,  appears  to  be  the  most 
sensitive  (conferring  the  least  buffering  capacity),  although  the 
subalpine  (mean  elevation  =  2651  m)  lakes  located  in  the  andesitic- 
volcaniclastic  rocks  of  the  northwestern  corner  also  may  be  poten- 
tially sensitive. 

2.  Lakes  found  in  other  regions  of  the  park--the  andesitic  Absaroka 
Mountains,  the  mixed  metamorphic  Gallatin  and  Beartooth  ranges- 
appear  not  immediately  sensitive  to  acidification  by  acid  precipita- 
tion (alkalinities  generally  above  500  peq/1).  These  findings  are 
corroborated  by  a  streamwater  chemistry  study  in  the  Absaroka 
Mountains  (Miller  and  Drever  1977b),  where  mean  alkalinity  for  14 
stream  samples  in  the  Shoshone  River  Basin  was  found  to  be  650 
ueq/1 . 

3.  Finally,  no  clear  trends  in  surface  water  chemistry  appear  from 
examination  of  historical  data.  It  is  suggested  that  if  any  chem- 
istry changes  have  occurred  at  all,  these  are  limited  to  the  north- 
west corner  lakes.  Since  estimated  "baseline"  (precipitation- 
derived)  sulfate  levels  in  Yellowstone  lakes  are  still  fairly  low  in 
comparison  to  those  levels  in  known  acidified  lakes  throughout  the 
continent,  any  acidification  that  has  occurred  is  minimal.  Most  of 
the  park  lakes  appear  well  protected  from  acidification  in  the 
future.  For  the  present,  the  dilute  lakes  of  the  rhyolite  bedrock 
and  Fall  River  basalt  region,  as  well  as  the  Specimen  Creek  drainage 
basin  lakes,  are  most  vulnerable  to  changes  in  chemistry  by  acid 
rain. 


EFFECTS  OF  CURRENT  AND  FUTURE  ACIDIFICATION  LEVELS  IN  FISH  POPULATIONS 
IN  YELLOWSTONE  AND  ROCKY  MOUNTAIN  NATIONAL  PARKS 

1.  Currently,  there  is  no  evidence  of  chronic  acidification  in  Rocky 
Mountain  waters  that  would  be  detrimental  to  fish  survival. 

2.  The  very  low  base  cation  concentrations  observed  in  the  headwater 
drainages  of  RMNP  suggests  extreme  sensitivity  to  acidification. 
Fish  populations  present  in  these  low  calcium  waters  may  be  particu- 
larly susceptible  to  osmoregulatory  stress  from  episodic  acidifica- 
tion. 

3.  The  few  remaining  native  trout  (Salmo  clarki )  populations  located  in 
interior  regions  of  the  Rocky  Mountains  persist  only  in  small, 
isolated  headwater  drainages.  The  cutthroat  coexists  with  and 
dominates  introduced  species  only  in  cold,  headwater  situations. 
Displacement  of  these  rare  and  endangered  genotypes  to  headwater 
drainages  also  makes  them  most  susceptible  to  potential  acidifica- 
tion in  these  sensitive  habitats. 
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4.  Although  a  large  proportion  of  the  headwater  lakes  in  YNP  and  other 
areas  are  historically  fishless,  and  many  are  still  in  that  condi- 
tion, these  waters  represent  potentially  invaluable  refugia  for  the 
reestabl ishment  of  endangered  species. 

5.  Comparisons  of  alkalinity  and  bedrock  geology  in  the  upper  Colorado 
and  eastern  Divide  drainages  in  RMNP,  suggests  that  the  greenback 
trout  (Sal mo  clarki  stomias) ,  which  occupies  the  latter  basins,  is 
at  potentially  greater  risk  than  the  Colorado  River  cutthroat  (Salmo 
clarki  pleuriticus). 

6.  A  comparison  of  the  life  history  patterns  of  the  western  Salmo  with 
eastern  salmonids,  in  relation  to  seasonal  changes  in  acidity, 
indicates  that  different  life  history  stages  may  be  impacted  by 
acidification.  Hatching  and  fry  development  have  been  identified  as 
critical  periods  for  the  fall  spawning  eastern  salmoids.  However, 
early-late  summer  emergence  of  fry  in  the  western  Salmo  populations 
(spring  spawners)  indicates  this  may  be  a  less  critical  life  history 
stage. 

7.  Avoidance  response,  in  the  form  of  downstream  emigration,  to  epi- 
sodic acidification  in  headwaters  occupied  by  S.  clarki  is  suggested 
as  a  more  subtle,  but  potentially  devasting  impact  of  acidification 
in  the  Rocky  Mountain  watersheds. 

8.  The  relative  sensitivities  of  the  cutthroat  genotypes  to  acidifica- 
tion stress,  capacities  for  adaption  in  low  calcium  water,  and 
dose-response  functions  are  currently  unavailable. 

SENSITIVITY  EVALUATION  OF  THE  CENTRAL  ROCKY  MOUNTAIN  REGION 

The  examination  of  the  geochemistry  of  Rocky  Mountain  National  Park  has 
shown  that  many  areas  in  RMNP  are  sensitive  to  acidic  deposition  and  that  this 
sensitivity  is  primarily  determined  by  bedrock  geology.  In  addition,  sensi- 
tivity varies  inversely  with  elevation  in  watersheds  with  consistent  geology. 
The  results  of  the  evaluation  of  the  existing  data  in  YNP,  though  not  so 
definitive,  corroborate  the  general  findings  in  the  RMNP  studies.  An  evalua- 
tion of  the  sensitivity  of  the  Central  Rocky  Mountains  (Colorado  and  Wyoming), 
using  geologic  maps  and  elevation,  can  be  provided  based  on  this  information. 

The  analyses  in  RMNP  show  that  watersheds  underlain  by  granite  and  bio- 
tite  gneiss  and  schist  are  equally  sensitive  to  acidic  deposition.  The  lakes 
and  streams  in  these  watersheds  had  alkalinities  <200  ueq/1  ,  while  the  waters 
at  higher  elevations  (>3300  m)  were  very  sensitive  (alkalinity  i  100  ueq/1). 
The  Upper  Colorado  River  Basin  and  the  Upper  Fall  River  Basin  contain  tertiary 
intrusive  rocks  in  their  drainage,  resulting  in  low  (alkalinity  >200  ueq/1) 
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sensitivity.   In  YNP,  low-alkalinity  or  sensitive  lakes  were  found  in  regions 
underlain  by  rhyolite  flows  or  basalt  and  andesi tic-basal  tic  flows  (USGS  1977; 

Tweto  1979). 

The  analysis  of  sensitive  aquatic  systems  in  the  Central  Rocky  Mountain 
Region  has  been  accomplished  by  extrapolating  the  results  from  RMNP  and  thus 
delineates  areas  underlain  by  granite  biotite  gneiss  and  schist  and  similar 
gneisses  and  schists  (Figures  30  and  31).  Areas  underlain  by  these  formations 
are  classified  as  sensitive  (alkalinity  <200  ueq/1),  lakes  and  streams  located 
at  higher  elevations  (>3300  m)  can  be  classified  as  very  sensitive  (alkalinity 
<100  ueq/1).  Although  areas  of  YNP  have  a  limited  number  of  moderately  sensi- 
tive lakes  (>200  ueq/1 ),  areas  underlain  by  tertiary  intrusive  rocks  are 
generally  classified  as  nonsensitive  (alkalinity  >200  ueq/1). 


•  DENVER 
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Figure  30.   Spatial  distribution  of  sensitivity  in  the  central  Rocky 
Mountain  Region:  Colorado. 
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Figure  31.  Spatial  distribution  of  sensitivity  in  the  central  Rocky 
Mountain  Region:  Wyoming. 

Mountain  ranges  underlain  by  rock  formations  similar  to  those  in  RMNP 
include  the  Sawatch  Range,  the  Front  Range  (Lovering  and  Goddard  1950),  and 
the  Sangre  de  Cristo  Mountains  (Litsey  1958)  in  Colorado  (Figure  30),  and  the 
Medicine  Bow  Mountains  (Hill  et  al.  1968)  and  Wind  River  Range  (Bassett  and 
Giletti  1963)  in  Wyoming  (Figure  31).  All  of  these  mountains  contain  granite, 
gneisses,  schists,  and  tertiary  intrusive  rocks. 

Although  this  method  of  evaluation  is  adequate  for  a  regional-scale 
sensitivity  assessment,  sensitivity  evaluations  for  individual  lakes  and 
streams  must  be  carried  out  on  a  drainage  by  drainage  basis.  Small  outcrops 
of  highly  weatherable  rocks,  too  small  to  be  seen  on  regional-scale  maps,  may 
have  an  overwhelming  influence  on  the  sensitivity  of  a  watershed.  This  is  the 
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case  in  the  Upper  Fall  River  Basin  in  RMNP.  The  bedrock  in  this  basin  is 
primarily  granite  and  biotite  gneiss  and  schist.  However,  a  small  deposit  of 
tertiary  intrusive  rock  at  the  head  of  the  watershed  gives  the  Fall  River  a 
relatively  high  alkalinity  (>200  ueq/1).  Differences  in  hydrologic  flow  path 
and  soil  development  may  also  dominate  sensitivity  on  a  local  scale. 


101 


RECOMMENDATIONS  FOR  RESEARCH  AND  MANAGEMENT 

INTRODUCTION 

While  these  studies  and  most  others  would  indicate  that  areas  of  the 
Rocky  Mountain  West  were  not  experienceing  significant  impacts  from  acid 
deposition,  they  also  demonstrate  that  many  headwater  lakes  and  streams  are 
very  sensitive  (al kalinities  <  200  peg/e).  It  is  therefore  believed  prudent 
to  undertake  some  long  term  programs  to  more  fully  assess  the  current  status 
of  both  deposition  and  surface  water  chemistries  and  to  develop  some  long  term 
measurement  programs.  In  addition  strategies  should  be  developed  to  protect 
indigenous  fish  populations  in  the  advent  that  this  area  does  experience 
increased  sulfate  and  nitrate  deposition.  The  following  recommendations  are 
formulated  to  address  these  issues. 

ROCKY  MOUNTAIN  NATIONAL  PARK  AND  ROCKY  MOUNTAIN  REGION 

1.  Long  term  watershed  experiments  should  be  established  at  several 
points  in  the  Rocky  Mountain  Region.  At  minimum,  routine  sampling 
(on  at  least  a  bi-weekly  basis  and  in  the  spring,  every  2  to  3  days) 
should  be  performed.  The  samples  should  be  analyzed  for  all  major 
inorganic  constituents,  DOC,  aluminum,  alkalinity  and  total  acidity. 
Based  on  the  research  of  this  project,  we  recommend  that  the  Glacier 
Gorge  watershed  in  RMNP  be  considered  for  long-term  monitoring.  (As 
a  result  of  this  recommendation  the  National  Park  Service  is  now 
conducting  a  long-term  study  in  the  Loch  Vale  subbasin  of  the 
Glacier  Gorge  watershed.) 

2.  Additional  surveys  of  water  chemistry  should  be  performed  in 
mountainous  areas  not  only  in  the  Rocky  Mountain  region  but  also  in 
other  mountainous  areas  of  the  western  United  States.  It  is  most 
probable  that  the  acidity  of  precipitation  will  increase  in  future 
years  and  it  will  be  to  our  benefit  to  obtain  background  data  at 
this  time. 

3.  The  results  from  the  watershed  studies  should  be  made  available  to 
research  groups  having  models  on  watershed  response  to  acidic  depo- 
sition. This  will  ensure  that  if  there  are  pecularities  about 
western  watersheds,  that  the  models  will  be  developed  with  those 
pecularities  taken  into  account. 
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4.  Since  one  of  the  differences  between  the  Rocky  Mountain  region  and 
other  regions  experiencing  acid  rain  is  the  greater  importance  of 
nitric  acid  relative  to  sulfuric  acid,  watershed  studies  should  be 
designed  to  take  this  into  account. 

5.  A  greater  effort  must  be  made  to  determine  the  rate  of  dry  deposi- 
tion of  neutral  salts  (e.g.,  CaSCL). 

YELLOWSTONE  NATIONAL  PARK 

1.  Future  studies  in  Yellowstone  National  Park  should  concentrate  on 
those  regions  which  this  report  has  determined  as  potentially  sensi- 
tive to  acidification.  In  particular,  reliable  water  chemistry  from 
the  Specimen  Creek  drainage  basin  is  necessary  to  determine  if  these 
lakes  are  indeed  as  dilute  as  indicated. 

2.  A  more  intensive  sampling  of  the  low  alkalinity  lakes  (i.e.,  several 
alkalinity  measurements  over  the  course  of  a  year)  would  indicate 
seasonal  fluctuations  in  alkalinity  and  provide  a  more  complete 
estimate  of  sensitivity.  It  is  suggested  that  the  analytical  tech- 
nique used  to  determine  alkalinity  of  these  lakes  be  one  designed 
specifically  for  low  alkalinity  water,  i.e.,  Gran's  plot  or  double 
endpoint  potentiometric  titration. 

3.  Finally,  a  survey  of  the  headwater  streams  of  the  Park,  particularly 
those  in  sensitive  regions  or  with  important  fisheries,  is  necessary 
for  a  complete  understanding  of  the  response  of  Yellowstone  to 
potential  acidification. 

FISHERY  RESEARCH  AND  MANAGEMENT 

1.  The  potential  sensitivity  to  acidification,  of  watersheds  currently 
occupied  by  endangered  or  threatened  populations  of  S.  clarki  or  of 
watersheds  considered  as  potential  candidate  sites  for  reintroduc- 
tion,  should  be  determined  and  given  consideration  in  recovery  plans 
for  the  species. 

2.  Any  further  experiments  designed  to  develop  dose-response  functions 
for  interior  western  Sal  mo  clarki  populations,  should  consider  the 
possibility  of  significant  variations  in  tolerance  of  this  polytypic 
species.  Additionally,  these  responses  to  acidification  need  to  be 
determined  in  very  low  calcium  media,  typical  of  the  headwater 
habitats  of  this  species. 

3.  The  behavioral  responses  and  emigration  tendencies  of  S.  clarki 
populations  exposed  to  episodes  of  acidification  should  be  deter- 
mined experimentally. 

4.  The  potential  for  aluminum  mobilization  in  Rocky  Mountain  watersheds 
exposed  to  increased  acid  deposition  should  be  determined. 
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5.  Future  studies  of  winter  and  spring  lake  and  stream  chemistry  in 
headwater  catchments  should  also  evaluate  movement  of  trout  popula- 
tions, in  response  to  chemical  change. 
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APPENDIX  A 

SOIL  AND  WATER  CHEMISTRY  SAMPLING  SITE  MAPS, 

ROCKY  MOUNTAIN  NATIONAL  PARK 
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Figure  A-l. 
RMNP. 


Soil  and  chemistry  sampling  sites  in  Upper  Colorado  basin, 
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Figure  A-3.      Soil   and  chemistry  sampling  sites   in  Glacier  Gorge  basin,    RMNP. 
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Figure  A-4.  Soil  and  chemistry  sampling  sites  in  Fall  River  basin,  RMNP. 
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APPENDIX  B 
SOIL  LABORATORY  PROCEDURES,  PROFILE  DESCRIPTIONS  AND  CHEMICAL  ANALYSES 

SOIL  LABORATORY  ANALYSIS  PROCEDURES 

Cation  Exchange  Capacity 

Ex  tract  able  acidity.  To  10  g  of  air  dry  soil  add  100  ml  of  extraction 
solution  (BaCl2,  triethanolamine).  Stopper,  shake  well,  and  let  stand  over- 
night. Run  blank.  Filter,  wash  with  extracting  solution.  Bring  up  to  250 
ml.  Add  indicator  (bromcresol  green  and  methyl  red)  and  titrate  with  0.2  N 
HC1  to  pink  orange  end-point.  Calculate.  (Black  1965) 

Base  extraction.  Use  macerated  paper  in  60  ml  fritted  glass  funnels. 
Extract  with  a  total  of  100  ml  1  N  ammonium  acetate  pH  7.0  (allow  to  sit 
overnight  in  30  ml  extracting  solution).  Run  Ca,  Mg,  K,  and  Na  on  extract. 
Calculate  C.E.C.  and  report  in  meq.  Calculate  percent  base  saturation. 
(Black  1965) 

Potassium:  1.5  g  soil  in  15  ml  of  1  N  ammonium  acetate,  pH  7.0.  Shake  5 
mm  at  200  opm.  Filter.  Read  on  AA.  Compare  with  standards.  (Carson 
1975) 

Calcium  and  magnesium:  To  1.5  g  soil  add  15  ml  IN  ammonium  acetate  pH 
7.0.  Shake  5  min  at  200  opm  and  filter.  Dilute  with  Lanthanum  Chloride 
(1500  ppm).  Read  on  AA  Spectrophotometer.  Compare  with  standards. 
(Walsh  1971). 

Sodium:  To  5  g  soil  add  15  ml  IN  ammonium  acetate  pH  7.0.  Shake  5  min 
at  200  opm.  Do  not  filter.  Read  supernatant  liquid  on  AA  Spectropho- 
tometer. Compare  with  standards.  (Walsh  1971). 

Organic  Matter 

Oxidize  with  potassium  dichromate  and  cone,  sulfuric  acid  by  standing  for 
30  min;  add  water,  phosphoric  acid,  sodium  flouride,  and  diphenylamine  indi- 
cator. Titrate  with  0.5  N  ferrous  ammonium  sulfate.   (Black  1965). 

£H 

1:1  soil  to  water  volume.  Mix  5  seconds.  Let  stand  10  min.  Stir,  read 
pH  with  pH  meter.  (McLean  1975). 
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Particle  Size 

Samples  are  sieved  to  remove  particles  larger  than  2  mm.  The  remaining 
sample  is  dispersed  with  sodium  hexametaphosphate,  placed  in  settling  con- 
tainers and  subsampled  at  appropriate  intervals  with  a  pipette.  The  subsample 
is  dried  at  105  C  and  weighed.   (Franzmeier  et  al.  1977). 

Loss  on  Ignition  (LOI) 

Ovendried  (105  L)  samples  are  ignited  in  a  muffle  furnace  at  550  C  and 
the  loss  in  weight  is  determined.  In  soils  with  low  clay  contents,  the  loss 
is  a  valid  estimate  of  organic  matter.  (Black  1965). 

Sulfate  Absorption 

Soil  samples  (10  g)  were  equilibrated  with  25  ml  of  dilute  K2S04  solution 
at  a  pH  of  4.1.  The  solutions  were  initially  at  1,  5  and  10  mg  sulfur  per 
liter.  After  filtration  the  solutions  were  analyzed  for  remaining  sulfate 
utilizing  an  indirect  Atomic  Absorption  method.   (0ien  1979). 
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Table 

B-l. 

Particle  size 

distribution 

of  soil 

samples. 

Depth 
(cm) 

>2  mm 
(%) 

Percent  of 

less  than  2 

mm 

Textural 
class 

Location 

Sand 

Silt 

Clay 

TG-1 

0-11 

55 

68 

27 

5 

SL 

TG-1 

11-25 

60 

82 

13 

5 

S 

TG-2 

0-17+ 

0 

35 

55 

10 

SIL 

TG-4 

0-5 

1 

21 

66 

13 

SIL 

TG-4 

5-11 

2 

21 

60 

19 

SIL 

TG-4 

11-24 

8 

35 

46 

19 

L 

AC-1 

0-8 

— 

52 

39 

9 

L 

GC-1 

Talus 

66 

91 

6 

3 

LS 

GC-2 

24-50+ 

56 

83 

14 

3 

S 

GC-3 

0-6 

3 

23 

70 

7 

SIL 

GC-3 

9.5-30 

19 

75 

15 

10 

SL 

EI-1 

Talus 

56 

65 

27 

8 

SL 

EI-2 

0-9 

0 

20 

75 

5 

SIL 

EI-3 

3-13 

58 

39 

45 

16 

L 

EI-3 

13-28+ 

76 

18 

70 

12 

SIL 

EI-5 

0-16 

8 

43 

35 

22 

L 

EI-5 

16-31+ 

10 

43 

34 

23 

L 

EI-6 

0-4 

23 

51 

41 

8 

L 

BX-1 

0-51+ 

53 

50 

35 

15 

L 

HG-2 

0-5 

7 

34 

54 

12 

SIL 

HG-2 

5-35+ 

70 

64 

23 

13 

SL 

MN-1 

Talus 

43 

70 

23 

7 

SL 

IL-2 

6-0 

— 

43 

54 

3 

SIL 

L80-1 

5-0 

— 

23 

73 

4 

SIL 

L80-1 

0-19 

— 

77 

19 

4 

S 

L80-1 

19-54 

— 

93 

4 

3 

LS 

L80-1 

54-80 

— 

84 

12 

4 

S 

L80-1 

80-85+ 

— 

80 

16 

4 

S 

L80-6 

0-18 

— 

63 

28 

9 

s 

L80-6 

18-50 

— 

49 

46 

5 

SL 

L80-6 

50-80+ 

— 

85 

11 

4 

s 

L80-11 

0-35+ 

— 

85 

11 

4 

s 

L80-13 

Muck 

— 

37 

56 

7 

SIL 

L80-13 

Gravel 

— 

77 

16 

7 

S 

L81-13 

0-1 

— 

55 

37 

8 

SL 

L81-13 

1-18 

— 

61 

32 

7 

SL 

L81-13 

18-46 

— 

62 

24 

14 

SL 

L81-13 

46-60+ 

— 

73 

17 

10 

SL 

L81-16 

31-42+ 

— 

32 

62 

6 

SIL 

L81-18 

0-24 

— 

68 

25 

7 

SL 

L81-18 

24-36 

— 

78 

17 

5 

S 

L81-18 

36-59+ 

— 

68 

27 

5 

SL 

YCS-1 

4-11 

— 

62 

31 

7 

SL 

YCS-1 

11-21 

— 

83 

13 

4 

S 

YCS-1 

21-36 

— 

YCS-1 

45-60+ 

-- 

82 

15 

3 

s 

See  maps  Appendix  A. 
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APPENDIX  C 
WATER  CHEMISTRY  DATA,  ROCKY  MOUNTAIN  NATIONAL  PARK 

The  raw  water  chemical  data  collected  in  Rocky  Mountain  National  Park  are 
presented  with  the  elevation  of  each  sample  site.  Missing  data  are  denoted  by 
an  asterisk.  Sample  location  maps  are  included  in  Appendix  A.  The  sample 
identification  code  follows: 

GA  Roaring  River 

GB  Ypsilon  Creek  Fall  River  Basin 

GC  Upper  Fall  River 

MA  Andrews  Creek 

MC  Loch  Vale  Glacier  Gorge 

MD  Glacier  Gorge 

V  Upper  Colorado  River  Basin  MW  East  Inlet 

The  sample  code  indicates  whether  the  sample  was  taken  from  a  stream  (S), 
lake  inlet  (I),  lake  outlet  (0),  or  lake  surface  (L).  Replicate  samples  are 
denoted  by  an  'R'  after  the  sample  number. 

The  lake  inlet  and  stream  samples  were  used  to  calculate  the  mean  concen- 
trations of  the  major  anions  and  cations  for  each  subbasin. 
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Table  Ol.   Lake  and  stream  pH  and  alkalinity  data 
for  Rocky  Mountain  National  Park. 


Sample 
ID 

Eleva 

F  pHb 

Temp 

Condd 

L  pHe 

Alkf 

Roaring  River 

GA01I 

3511.30 

6.23 

2.0 

12.7 

6.45 

26.0 

GA02S 

3499.10 

6.95 

8.0 

19.2 

6.98 

96.0 

GA04L 

3511.30 

6.63 

13.0 

9.6 

6.81 

47.0 

GA050 

3511.30 

6.50 

11.0 

10.2 

6.66 

50.0 

GA080 

3511.30 

6.50 

11.5 

10.6 

6.66 

44.0 

GA10I 

3364.99 

7.05 

11.0 

19.2 

7.09 

98.0 

GA12L 

3352.80 

6.03 

14.0 

15.8 

6.70 

81.0 

GA160 

3352.80 

6.19 

14.0 

14.3 

6.89 

83.0 

GA18S 

3291.84 

6.08 

14.0 

10.0 

6.98 

83.0 

GA20S 

3194.30 

6.65 

13.0 

17.1 

6.96 

78.0 

GA24S 

3169.92 

6.47 

13.0 

16.0 

6.78 

75.0 

GA26S 

3017.52 

6.30 

13.5 

15.1 

6.96 

95.0 

GA28S 

2901.70 

6.40 

13.0 

13.5 

6.87 

77.0 

GA30S 

2804. 16 

6.09 

13.0 

13.2 

7.01 

72.0 

GA32S 

2621.28 

6.14 

8.0 

15.4 

6.97 

92.0 

GA32SR 

2621.28 

6.14 

8.0 

15.4 

6.87 

97.0 

GA46S 

2926.08 

6.00 

14.0 
Ypsilon  Creek 

12.3 

6.83 

57.0 

GB02I 

3413.76 

5.90 

5.0 

— 

6.43 

31.0 

GB060 

3413.76 

5.63 

8.0 

21.3 

6.56 

31.0 

GB08I 

3352.80 

5.96 

9.0 

20.9 

6.56 

35.0 

GB09I 

3352.80 

6.20 

6.0 

— 

6.78 

66.0 

GB10L 

3352.80 

6.05 

13.0 

19.9 

6.65 

42.0 

GB120 

3352.80 

5.94 

11.0 

19.9 

6.66 

42.0 

GB14I 

3279.65 

6.27 

13.0 

21.9 

6.68 

60.0 

GB160 

3279.65 

6.50 

11.0 

21.4 

6.75 

54.0 

GB18S 

3108.96 

— 

— 

11.8 

6.69 

63.0 

GB19I 

3462.53 

6.41 

11.0 

16.5 

6.09 

16.0 

GB220 

3462.53 

6.41 

8.0 

19.1 

6.05 

19.0 

GB260 

3462.53 

6.58 

12.0 

20.4 

6.53 

26.0 

GB30L 

3462.53 

6.83 

14.0 

22.2 

6.58 

34.0 

GB320 

3462.53 

6.78 

13.0 

19.0 

6.61 

38.0 

GB34S 

3352.80 

6.12 

8.0 

22.7 

6.56 

33.0 

GB35S 

3352.80 

6.17 

8.0 

20.4 

6.33 

28.0 

GB38I 

3218.69 

6.45 

8.0 

15.9 

6.65 

33.0 

GB40L 

3218.69 

7.00 

13.0 

15.1 

6.59 

42.0 

GB40LR 

3218.69 

7.00 

13.0 

15.1 

6.65 

36.0 

GB420 

3218.69 

6.83 

13.0 

15.5 

6.56 

37.0 
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Tabl 

e  C-l.  (continued) 

Sample 

p> 

b 

r 

d 

p 

f 

ID 

Eleva 

F  pHD 

Temp 

Conda 

L  pHe 

AlkT 

Upper  Fall  River 

GC02S 

3413.76 

0.00 

7.0 

42.2 

7.32 

377.0 

GC02SR 

3413.76 

0.00 

7.0 

42.2 

7.61 

358.0 

GC03S 

3413.76 

0.00 

6.0 

26.5 

7.12 

172.0 

GC04S 

3535.68 

0.00 

6.0 

10.0 

6.88 

69.0 

GC06S 

3523.49 

6.99 

5.0 

15.0 

7.13 

138.0 

GC08S 

3523.49 

6.89 

5.0 

39.6 

7.47 

357.0 

GC10S 

3511.30 

6.90 

5.0 

24.8 

7.07 

216.0 

GC11S 

3401.57 

0.00 

10.0 

9.6 

6.65 

51.0 

GC12S 

3401.57 

0.00 

10.0 

22.5 

7.07 

184.0 

GC13S 

2926.08 

0.00 

8.0 

23.6 

7.38 

170.0 

GC14S 

2816.35 

0.00 

9.0 

14.6 

7.05 

119.0 

GC16S 

2767.58 

0.00 

9.0 

21.9 

7.21 

162.0 

GC18S 

2657.86 

6.13 

7.0 

22.1 

7.17 

167.0 

GC20S 

2657.86 

6.23 

8.0 

14.3 

6.42 

87.0 

GC24S 

2657.86 

6.20 

8.0 
Tyndall  Gorge 

20.4 

6.93 

146.0 

MA02I 

3084.58 

5.66 

8.0 

27.4 

6.36 

29.0 

MA02IR 

3084.58 

5.66 

8.0 

27.4 

6.36 

24.0 

MA060 

3084.58 

5.61 

12.0 

23.4 

6.56 

31.0 

MA08I 

3023.62 

5.81 

11.0 

27.0 

6.46 

39.0 

MA120 

3023.62 

5.77 

15.0 

26.6 

6.77 

38.0 

MA16S 

2804.16 

Loch  Vale 

12.0 

6.60 

62.0 

MC02I 

3474.72 

6.66 

1.5 

53.3 

6.46 

46.0 

MC04L 

3474.72 

6.66 

9.0 

30.9 

6.46 

28.0 

MC06O 

3474.72 

6.48 

8.0 

31.3 

5.98 

36.0 

MC12S 

3267.46 

6.55 

9.0 

34.1 

6.38 

27.0 

MC14S 

3169.92 

5.92 

6.0 

29.5 

6.43 

32.0 

MC14SR 

3169.92 

5.92 

6.0 

29.5 

6.47 

29.0 

MC18I 

3328.42 

6.21 

2.0 

36.3 

6.04 

17.0 

MC20L 

3328.42 

0.00 

— 

— 

6.53 

22.0 

MC220 

3328.42 

7.05 

11.0 

31.9 

6.42 

32.0 

MC24I 

3316.22 

6.24 

12.0 

35.7 

6.44 

24.0 

MC280 

3316.22 

6.04 

13.0 

34.5 

— 

32.0 

MC30I 

3108.96 

6.51 

13.0 

38.8 

6.46 

40.0 

MC32L 

3108.96 

6.81 

19.0 

37.1 

6.74 

43.0 

MC340 

3108.96 

5.93 

15.0 

33.4 

6.66 

47.0 

MC36S 

2987.04 

5.94 

16.0 

10.0 

6.57 

43.0 
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Table  C 

-1.  (continued) 

Sample 
ID 

Eleva 

F  PHb 

Temp 

Condd 

L  pHG 

Alkf 

G1 

acier  Creek 

MD02I 

3535.68 

6.10 

8.0 

28.0 

5.66 

10.0 

MD060 

3535.68 

6.15 

12.0 

18.5 

5.89 

11.0 

MD08I 

3230.88 

6.22 

11.0 

26.7 

— 

— 

MD10I 

3230.88 

6.21 

6.0 

20.8 

6.09 

30.0 

MD140 

3230.88 

6.34 

11.0 

24.4 

6.56 

52.0 

MD24I 

3474.72 

6.17 

3.0 

33.7 

6.52 

33.0 

MD26L 

3474.72 

6.82 

12.0 

26.6 

6.66 

38.0 

MD280 

3468.62 

6.49 

12.0 

28.2 

6.66 

35.0 

MD32L 

3413.76 

6.89 

14.0 

29.3 

6.65 

43.0 

MD340 

3413.76 

6.90 

13.0 

23.9 

6.52 

64.0 

MD38S 

3145.54 

6.47 

14.0 

25.5 

6.63 

53.0 

MD40S 

3035.81 

6.44 

14.0 

26.5 

6.52 

54.0 

MD420 

3035.81 

6.51 

16.0 

25.7 

6.79 

47.0 

MD44L 

3035.81 

6.48 

18.0 

26.5 

6.66 

50.0 

MD460 

3035.81 

6.61 

15.0 

26.3 

6.75 

45.0 

MD48S 

2974.85 

5.97 

16.0 

9.9 

6.73 

65.0 

MD50S 

2865.12 

5.88 

12.0 

10.1 

6.73 

65.0 

MD52S 

2755.39 

0.00 

0.0 

10.1 

0.00 

53.0 

Upper 

Colorado  River 

V01I 

3511.30 

0.00 

11.0 

10.5 

6.91 

85.0 

V02S 

3230.88 

0.00 

0.0 

76.0 

7.93 

682.0 

V030 

3499.10 

0.00 

0.0 

30.2 

7.38 

280.0 

V04S 

3255.26 

0.00 

0.0 

29.6 

7.62 

242.0 

V06S 

3182.11 

0.00 

0.0 

53.0 

7.61 

427.0 

V08S 

3169.92 

0.00 

0.0 

56.4 

7.69 

607.0 

VIOS 

3108.96 

0.00 

0.0 

46.6 

7.58 

375.0 

VIOSR 

3108.96 

0.00 

0.0 

0.0 

7.61 

366.0 

V12S 

3115.06 

6.95 

10.0 

72.9 

8.03 

580.0 

VMS 

3121.15 

6.71 

9.0 

64.8 

7.56 

374.0 

V16S 

3121.15 

6.70 

10.0 

52.0 

7.44 

361.0 

V18I 

3486.91 

6.82 

9.0 

18.7 

7.29 

146.0 

V20L 

3486.91 

9.00 

11.0 

41.6 

7.72 

342.0 

V220 

3486.91 

6.85 

110.0 

41.6 

7.65 

341.0 

V24S 

3230.88 

7.19 

7.0 

45.8 

7.58 

357.0 

V25S 

3121.15 

0.00 

0.0 

0.0 

7.69 

367.0 

V28S 

3121.15 

6.70 

8.0 

64.0 

7.76 

480.0 

V28SR 

3121.15 

0.00 

0.0 

0.0 

7.76 

480.0 

V30S 

3133.34 

0.00 

0.0 

0.0 

7.33 

202.0 

V32S 

3133.34 

0.00 

0.0 

26.1 

6.56 

50.0 

V34S 

3133.34 

0.00 

0.0 

0.0 

0.00 

0.0 

V34SR 

3133.34 

0.00 

0.0 

0.0 

0.00 

0.0 

V36S 

3133.34 

0.00 

0.0 

48.9 

7.78 

443.0 

V38S 

3108.96 

6.40 

9.0 

22.9 

7.20 

189.0 
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Table  Ol.  (concl 

uded) 

Sample 

a 

b 

r 

H 

p 

,  f 

ID 

Eleva 

F  pHD 

Temp 

Conda 

L  pHe 

AlkT 

i 

East  Inlet 

MW02I 

3316.22 

5.78 

7.0 

13.0 

0.00 

80.0 

MW060 

3316.22 

8.13 

11.0 

12.0 

6.92 

80.0 

MW08I 

3169.92 

0.00 

12.0 

11.4 

6.79 

72.0 

MW120 

3169.92 

0.00 

16.0 

10.9 

6.84 

93.0 

MW15I 

3145.54 

0.00 

8.0 

8.8 

7.00 

94.0 

MW16I 

3145.54 

0.00 

15.0 

10.2 

6.80 

68.0 

MW200 

3145.54 

0.00 

16.0 

9.8 

7.10 

76.0 

MW22I 

3108.96 

0.00 

13.0 

9.3 

6.73 

60.0 

MW260 

3108.96 

0.00 

15.5 

10.0 

6.84 

83.0 

MW30S 

3084.58 

0.00 

0.0 

11.6 

6.96 

77.0 

MW32I 

3023.62 

0.00 

0.0 

12.6 

6.85 

73.0 

MW360 

3011.42 

0.00 

17.0 

10.2 

6.81 

91.0 

MW38S 

2865.12 

0.00 

0.0 

12.3 

6.76 

80.0 

MW40S 

2889.50 

0.00 

0.0 

15.4 

6.45 

121.0 

MW40SR 

2889.50 

0.00 

0.0 

15.4 

6.86 

120.0 

MW42S 

2865.12 

0.00 

0.0 

14.4 

6.64 

94.0 

MW44S 

2791.97 

0.00 

0.0 

13.4 

7.10 

96.0 

*Elev  = 

FpH  = 

dTemp  = 

elevation 

in  me 

ters. 

field  pH  in  SU. 

temperature  in 

°c. 

Cond  = 
ft-  PH  = 

conductivi 

ty  in 

umohs/cm, 

lab  pH  in 

SU. 

Alk  =  alkalinity  in  ueq/1 . 
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REGION  3 

Regional  Director 
U.S.  Fish  and  Wildlife  Service 
Federal  Building,  Fort  Snellinj 
Twin  Cities,  Minnesota  55111 


REGION  4 

Regional  Director 
U.S.  Fish  and  Wildlife  Service 
Richard  B.  Russell  Building 
75  Spring  Street,  S.W. 
Atlanta,  Georgia  30303 


REGION  5 

Regional  Director 

U.S.  Fish  and  Wildlife  Service 

One  Gateway  Center 

Newton  Corner,  Massachusetts  02  1  58 


REGION  6 

Regional  Director 

U.S.  Fish  and  Wildlife  Service 

P.O.  Box  25486 

Denver  Federal  Center 

Denver,  Colorado  80225 


REGION  7 

Regional  Director 
U.S.  Fish  and  Wildlife  Service 
1011  b.  Tudor  Road 
Anchorage,  Alaska  99503 


As  the  Nation's  principal  conservation  agency,  the  Department  of 
the  Interior  has  responsibility  for  most  of  our  nationally  owned 
public  lands  and  natural  resources.  This  includes  fostering  the 
wisest  use  of  our  land  and  water  resources,  protecting  our  fish  and 
wildlife,  preserving  the  environmental  and  cultural  values  of  our 
national  parks  and  historical  places,  and  providing  for  the  enjoy- 
ment of  life  through  outdoor  recreation.  The  Department  assesses 
our  energy  and  mineral  resources  and  works  to  assure  that  their 
development  is  in  the  best  interests  of  all  our  people.  The  Depart- 
ment also  has  a  major  responsibility  for  American  Indian  reservation 
communities  and  for  people  who  live  in  island  territories  under  U.S. 
administration. 
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